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Abstract 

We discuss the importance of high-p^ hadron and jet measurements in nucleus- 
nucleus collisions at the LHC. 
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1. INTRODUCTION 



U.A. Wiedemann 

This report summarizes the current understanding of how the production of high-pr partons in 
nucleus-nucleus collisions at the LHC can be used as a "hard probe", i.e. as a diagnostic tool for the 
produced QCD matter either in thermaUzed (quark gluon plasma) or in other non-equilibrated but dense 
forms. 

The production of high-pT partons (observed as high-pT hadrons or jets) involves a "hard" per- 
turbative scale Q 3> Aqcd- This report mainly considers the case when this scale is harder than any 
momentum scale characterizing the medium produced in the nucleus-nucleus collision. Momentum 
scales proposed to characterize the medium (such as the initial temperature T or the saturation momen- 
tum Qs) may be perturbatively large themselves, in which case Q ^ T,Qs and T, Qg ^ Aqcd- This 
restriction aims at insuring that hard partonic production processes are not part of the "bulk matter": they 
occur in the primary partonic colhsions on temporal and spatial scales At ^ 1/Q, Ar ^ 1/Q which are 
sufficiently small to be unaffected by the properties of the produced matter. This makes them promising 
candidates of processes whose primary partonic production process is unaffected by the presence of a 
medium, while the development of the final (and possibly initial) state parton shower leaving (enter- 
ing) the hard partonic subprocess is sensitive to the medium. If colhnear factorization is applicable in 
nucleus-nucleus collisions, then inclusive cross sections of high-p^ partons measured in proton-proton 
colhsions or calculated in perturbative QCD can be used as benchmark against which one can search for 
the actual signals and properties of the hot and dense matter created in nucleus-nucleus collisions at the 
LHC. 

Section 2. discusses benchmark calculations for jet spectra and identified high-pr hadronic spectra 
in nucleus-nucleus collisions at LHC, calculated in the framework of collinear factorized QCD. The 
question to what extent collinear factorization can be expected to hold in nucleus-nucleus collisions, and 
how its validity can be tested experimentally, is discussed in chapter [1] of this workshop report. 

Section 3. addresses the main theoretical arguments for strong final state medium-modifications 
in jet production at the LHC. A jet is the hadronized remnant of a final state parton shower related to 
a produced highly virtual parton. Sections 3.1 and 3.2 focus mainly on the modification of this final 
state parton shower due to multiple parton scattering in a spatially extended dense medium. The current 
understanding of the additional medium-induced radiative contributions and the transverse momentum 
broadening of the parton shower is discussed. Section 3.3 compares radiative and coUisional modifica- 
tions of the parton shower. Section 3.4 discusses calculations of the main observables in which these 
medium-modifications are expected to show up. Related results obtained in the formalism of medium- 
enhanced higher twist expansion are summarized in section 3.5. Finally, the section 3.6 on other poten- 
tially large medium-modifications discusses medium effects which may become important at moderate 
scales, Q up to ~ 10 GeV. 

Once medium-modifications of hard probes are determined, the question arises to what extent 
the properties of the hot and dense matter produced in nucleus-nucleus collisions differ from those of 
normal cold nuclear matter. To this end. Section 3.14 summarizes what is known about the medium- 
modifications of hard probes in cold nuclear matter. 

The remainder of this report summarizes the experimental situation. Section 4. gives a short 
overview of the data available for Au-i-Au collisions at ^/s = 200 GeV from the Relativistic Heavy Ion 
Colhder RHIC at Brookhaven. Several measurements at RHIC indicate that strong final state medium- 
modifications of the hadron production in central Au-Au colhsions persist at RHIC up to the highest 
transverse momenta explored {pT < 15 GeV). This further supports the theoretical expectations of 
sttong final state medium-modifications in nucleus-nucleus collisions at LHC where a much wider range 
in transverse energy is experimentally accessible. Finally, Section 5. discusses the current status of how 
the LHC experiments ALICE, ATLAS and CMS will measure jets and their medium-modifications in 



the high-multiplicity background of a heavy ion collision. 
In summary: 

• Jets and high-pT hadrons are the most abundant hard probes produced in nucleus-nucleus collisions 
at LHC. Within one month of running at design luminosity and in the absence of strong medium- 
modifications, jet spectra up to at least Ex = 200 GeV and leading hadron spectra up to pT = 100 
GeV are accessible. — > Section 2.. 

• Both theoretical arguments ( — > Section 3.) and data at lower center of mass energy ( — > Sec- 
tion 4.) suggest large medium-modifications of hadronic high-pr spectra in nucleus-nucleus colli- 
sions at LHC. 

• A variety of theoretical approaches (comparison with benchmark calculations — > Section 2. ) 
and experimental techniques (comparison to benchmark measurements, dependence of medium- 
modifications on nuclear geometry — > Section 3.46 , etc.) are available to quantify the medium- 
dependence of jet production at LHC. The mutual consistency of these different approaches is a 
prerequisite for any characterization of the produced hot and dense medium from the medium- 
dependence of hard probes (see also Ref. [2, 3]). We emphasize that data from p-A collisions at 
the LHC are an important part of this program [1]. 

2. BENCHMARK CROSS SECTIONS 

2.1 Jet and Dijet Rates in Nucleus-Nucleus Collisions 

A. Accardi, N. Armesto, I. P. Lokhtin 

Jet studies will play a central role as a proposed signature of the formation of QGP in AB colh- 
sions. Energy loss of energetic partons inside a medium where color charges are present, the so-called 
jet quenching [74], has been suggested to behave very differently in cold nuclear matter and in QGP. It 
has been postulated as a tool to probe the properties of this new state of matter [162, 280, 240, 263]. 

On the other hand, jet calculations at NLO have been successfully confronted with experimental 
data in hadron-hadron collisions [32]. Monte Carlo codes have become available: among them, we will 
use that of [138, 139, 140] adapted to include isospin effects and modifications of nucleon pdf inside 
nuclei, see the section on jet and dijet rates in pA collisions [1] for more information. Here we will 
present the results of 'initial' state effects, i.e. no energy loss of any kind will be taken into account. 
These results can be considered as the reference, hopefully to be tested in pA, whose failure should 
indicate the presence of new physics. As in pA collisions, we will work in the LHC lab frame, which for 
symmetric AB collisions coincides with the center-of-mass one, and the accuracy of our computations, 
limited by CPU time, is the same as in the pA case: 

• For the transverse energy distributions, 2 % for the lowest and 15 % for the highest Ex-him. 

• For the pseudorapidity distributions, 3 %. 

• For the dijet distributions of the angle between the two hardest jets, 20 % for the least populated 
and 3 % for the most populated bins. 

All the energies will be given per nucleon and, in order to compare with the pp case, all cross sections 
will be presented per nucleon-nucleon pair, i.e. divided by AB. 

Unless explicitly stated and as in the pA case, we will use as nucleon pdf MRST98 central gluon 
[214] modified inside nuclei using the EKS98 parameterizations [125, 126], a factorization scale equal 
to the renormalization scale /i = /Xi? = iir = Et/'2., with Et the total transverse energy of all the jets 
in the generated event, and for jet reconstruction we will employ the /cy-clustering algorithm [102, 121] 
with D = 1. The kinematic regions we are going to consider are the same as in the pA case: 

• |r7i| < 2.5, with r]i the pseudorapidity of the jet; this corresponds to the acceptance of the central 
part of the CMS detector. 

• Exi > 20 GeV in the pseudorapidity distributions, with Eti the transverse energy of the jet; this 
will ensure the validity of perturbative QCD. 
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Fig. 1: Isospin and nuclear pdf dependence of jet cross sections (pp results: solid lines; PbPb results without modification of 
nucleon pdf inside nuclei: dashed lines; PbPb results with EKS98 modifieation of nucleon pdf inside nuclei: dotted lines) versus 
transverse energy of the jet (for < 2.5, upper plot) and pseudorapidity of the jet (for En > 20 GeV, middle plot), and dijet 
cross sections (lower plot) versus angle between the two hardest jets for Sti > 20 GeV, i5T2 > 15 GeV and |?7i|, j?72| < 2.5, 
for collisions at 5.5 TeV. Unless otherwise stated default options are used, see text. 



• Eti > 20 GeV and Et2 > 15 GeV for the (/)-distributions, with Eti (£'t2) the transverse energy 
of the hardest (next-to-hardest) jet entering the CMS acceptance, and (p the angle between these 
two jets. 

For more information, we refer to the chapter [1] of this report. The centrality dependence is not studied 
in either contribution. 

The words of caution about our results which were given in the pA Section are even more relevant 
in AB collisions, as our ignorance on soft multiparticle production in this case is even larger than in pA 
colUsions. For example, the number of particles produced at midrapidity in a central PbPb collision at 
the LHC may vary as much as a factor 3 [58, 101] among different models which, in principle, are able 
to reproduce the available experimental data on multiplicities at SPS, RHIC and TeVatron. Therefore, 
these issues of the underlying event [166, 131] and multiple hard parton scattering [10, 11, 48, 5, 6] 
demand extensive Monte Carlo studies including full detector simulation. PreUminary analysis, based 
on the developed sliding window-type jet finding algorithm (which subtracts the large background from 
the underlying event) and full GEANT-based simulation of the CMS calorimetry, shows that even in the 
worst case of central PbPb collisions with maximal estimated charged particle density at mid-rapidity 
dN^/dy\y=o = 8000, jets can be reconstructed with almost 100 % efficiency, low noise and satisfactory 
energy and spatial resolution starting from Eti '-^ 100 GeV (see the Section on Jet Detection at CMS). 
In the case of more reaUstic, lower multiphcities, the minimal threshold for adequate jet reconstruction 
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Fig. 2: Jet cross sections versus transverse energy of the jet (for \rii\ < 2.5, plots on the left) and pseudorapidity of the jet (for 

En > 20 GeV, plots in the middle), and dijet cross sections versus angle between the two hardest jets for Eti > 20 GeV, 
Et2 > 15 GeV and |?7i|, |?72| < 2.5 (plots on the right), for PbPb collisions at 5.5 TeV (upper plots) and ArAr collisions at 6.3 
TeV (lower plots). Default options are used, see text. 



could even decrease. 

As in the pA case, see the previously mentioned section on pA collisions, the influence of discon- 
nected collisions on jet production in AB collisions may be studied using simple estimates on the number 
(n) of nucleon-nucleon collisions involved in the production of jets with Exi greater than a given Eto, 
which can be obtained in the Glauber model [100, 97] in the optical approximation: (n)(6, E'to) = 
ABTABib)a(ETo)/cTABib, Eto), with b the impact parameter, TAsip) = J (fsTA{s)TB{b - s) the 
convolution of the nuclear profile functions of projectile and target normalized to unity, a{ETo) the 
cross section for production of jets with Exi greater than Eto in PP collisions, and aAB{b,ETo) = 
1 - [1 - TABib)(T{ETo)]'^^ ■ Taking cr(^To) = 294, 0.463 and 0.0012 /ib as representative values in 
PbPb collisions at 5.5 TeV for Eto = 20, 100 and 200 GeV respectively (see results in Fig. 2 below), the 
number of nucleon-nucleon collisions involved turns out to be respectively 3.5, 1.0 and 1.0 for minimum 
bias collisions (i.e. integrating numerator and denominator in apA{b, Eto) between 6 = and oo), while 
for central collisions (integrating between 6 = and 1 fm) the numbers are 8.9, 1.0 and 1.0 respectively. 
So, in AB collisions at LHC energies the contribution of multiple hard scattering coming from different 
nucleon-nucleon collisions seems to be negligible for transverse energies of the jets greater than ~ 100 
GeV, while for Eti smaller than ~ 50 GeV this effect might need to be taken into account more carefully 
in our computations. 



Collision 


Ecm per nucleon (TeV) 




Number of jets/events per month per //b/(AB) 


ArAr 


6.3 




1.6 ■ 10« 


ArAr 


6.3 


3 • 10'-^'^ 


4.8 • 10*^ 


PbPb 


5.5 


5 • 10^*^ 


2.2 • 10^ 



Table 1: Luminosities C in units of cm~^s~^ and £ x 10^ s in units of fih/AB for different collisions at the LHC. 
The numbers of expected jets and dijets in a certain kinematic range are obtained by multiplying the latter column 
by the cross sections given in Figs. 1, 2 (jets) and 2 (dijets). 



2.11 Uncertainties 

Uncertainties on the renormalization/factorization scale, on the jet reconstruction algorithm and on nu- 
cleon pdf, have been discussed in the mentioned Section on pA collisions and show very similar features 

in the AB case, so we will discuss them no longer. Here we will focus, see Fig. 1, on isospin effects 
(obtained from the comparison of pp and PbPb without any modifications of nucleon pdf inside nuclei at 
the same energy per nucleon, 5.5 TeV) and on the effect of modifications of nucleon pdf inside nuclei, 
estimated by using EKS98 [125, 126] nuclear corrections. 

On the transverse momentum distributions isospin effects are negligible, while effects of EKS98 
result in a ~ 3 % increase. On the pseudorapidity distributions, isospin effects apparently tend to fill a 
small dip at ~ present in the pp distribution, while EKS98 results in some increase, but nevertheless 
effects never go beyond 5 % and are not very significant when statistical errors are considered. On the 
dijet angular distributions, isospin effects are negligible while EKS98 produces an increase of order 10 % 
at maximum. 

2.12 Results 

The expected LHC luminosities in different collisions are collected into Table 1, and also shown are 
£ X (1 month) in units of ijLh/{AB). Using this Table and the cross sections for inclusive one-jet 
production and dijet production in the Figures, it is possible to extract the number of expected jets (Figs. 1 
and 2) or dijets (Fig. 2) in different ranges of the kinematic variables. For example, examining the solid 
line in Fig. 2 (upper-left) one can expect, within the pseudorapidity region we have considered, the 
following number of jets per month in PbPb colUsions at 5.5 TeV with a luminosity of 5 • 10^^ cm~^s~^: 
2.2 • 10'' jets with Eti ~ 50 GeV (corresponding to a cross section of 1 /ib/(AB)), and 2.2 ■ 10^ jets with 
Eti ~ 250 GeV (corresponding to a cross section of 10^^ ;ub/(AB)). 

A detailed study of jet quenching [74, 162, 280, 240, 263] and of associated characteristics as 
jet profiles, which should be sensitive to radiation from the jet [73], should be feasible with samples of 
~ 10^ jets. Looking at the results given in Fig. 2, it becomes evident that, from a theoretical point of 
view, the study of such samples should be possible up to a transverse energy Exi ~ 275 GeV with a run 
of 1 month at the considered luminosity: indeed, from Table 1, 10^ jets for PbPb would correspond to a 
cross section of 4.5 • 10^^ /ib/(AB), which in Fig. 2 (upper-left) cuts the curve at Eti ~ 275 GeV. 

The centraMty dependence of the observables has not been examined due to our poor knowledge 
of the centrality behavior of the modification of nucleon pdf inside nuclei; if this behavior becomes clear 
in future experiments at eA colliders [59, 117, 8], such study would become very useful [204]. In any 
case, a variation of nuclear sizes should allow a systematic study of the dependence of jet spectra on the 
size and energy density of the produced plasma. 

2.2 Benchmark Particle Cross Sections 

Ivan Vitev 

Hadron production in leading order pQCD is reviewed. The shape of the single inclusive particle 
spectra is well described for > 2 — 3 GeV at center of mass energies from 20 GeV to 2 TeV. The 



phenomenological K-factor is found to decrease systematically with ^/s. For ultra-relativistic heavy ion 
reactions the calculation is augmented with the effects of initial multiple parton scattering and final state 
radiative energy loss. Baseline CERN-LHC predictions for hadron production inp + p and suppression 
in central Pb + Pb reactions at ^/s = 5.5 TeV are given in comparison to the corresponding results at 
BNL-RHIC and CERN-SPS energies. 

The purpose of this section is to present a lowest order (LO) analysis of inclusive hadron produc- 
tion up to the Tevatron energies and discuss hadron differential cross sections and composition at the 
LHC. 



2.21 Hadroproduction in Factorized pQCD 

The standard factorized pQCD hadron production formalism expresses the differential hadron cross 
section in N + N h + X as a convolution of the measured parton distribution functions (PDFs) 
fa/NiXai Q%) for the interacting partons {a = a, 6), with the fragmentation function (FFs) Df^/^iz, Ql) 
for the leading scattered parton c into a hadron of flavor h and the parton-parton differential cross sections 
for the elementary sub-process da^"'^^^'''^ / di: 



1 1 1 



da f f f 

Eh = Knlo^ J '^^'^ J J ■l'<^/p^^"-'^a)fb/p{xb,Q, 



abed o T ■ ■ 

^ •^-a mm '*'omm 



^ ^^{ab-^cd) 

xDh/cizc,Ql) — -f Sis + u + i). (1) 

A list of the lowest order partonic cross sections can be found in [230]. In Eq. {\) Xa,Xb are the initial 
momentum fractions carried by the interacting partons and Zc = Ph/Pc is the momentum fraction of the 
observed hadron. lo is a phenomenological factor that is meant to account for next-to-leading order 
(NLO) corrections. It is ^/s and scale dependent and takes typical values ~ 1 — 4. One usually finds 
that Eq. (1) over-predicts the curvature of the inclusive hadron spectra \dp,j,da'^\ at transverse momenta 
Pt < 4 GeV. This can be partly corrected by the introduction of a small intrinsic (or primordial) kx- 
smearing of partons, transversely to the collision axis, and generalized parton distributions fa{x, kT,Q^) 
motivated by the pQCD initial state radiation. For the corresponding modification of the kinematics in 
(1) in addition to the / d'^k^ J d'^k^{- ■ ■) integrations see [230]. The generaUzed parton distributions are 
often approximated as 

fa{x,kT,Q ) fa{x,Q )g{kT), gikx) = ^^^2^ ' 
where the width {k^) of the Gaussian enters as a phenomenological parameter. 



Perturbative QCD fits to data [112, 251, 51, 46, 37, 80, 89, 4] use different coupled choices for 
Knlo and (k^) and the extracted values are thus not directly comparable. However, similar agreement 
between data and theory at the level of spectral shapes and the ^/s dependence of the corrective factors 
discussed above is found. In [293] the factorization and fragmentation scales were set to Qpdf = Pt/2 
and Qff = Prf^Zc and no K^lo factors were employed. The extracted decreases from 2.7 GeV^ 
at ^/s ~ 50 GeV to 0.75 GeV^ at ~ 2 TeV. Alternatively, in [128] no primordial /cy-smearing was 
used and the scales in the calculation were fixed to be Qpdf = Qff = Pt- The deduced Knlo 
decreases from ~ 6 at ^/s ~ 50 GeV to ~ 1.5 at ^/s ~ 2 TeV. 

In the fits shown in Fig. 3 we have used the GRV98 LO PDFs [148] and the BKKLO FFs [84]. Pro- 
ton-i-antiproton fragmentation has been parameterized as in [275], inspired from PYTHIA [248] results. 
A fixed {k'^)pp = 1.8 GeV^ has been employed, leading to a K]\fLo parameter that naturally exhibits 
a smaller variation with ^/s. A ±25% error band about the Knlo value, fixed by the requirement to 
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Fig. 3: Extracted Knlq from comparison of LO pQCD calculation to data [51, 46, 37, 80, 89, 4] at and about 
mid-rapidity in the range 2 < px < 25 GeV. A systematic decrease of lo with y/s is observed and illustrated 
in the bottom right panel. The projected 50% uncertainty at -y/s = 5.5 — 8.8 TeV is also shown. 



match the moderate- and high-pr behavior of the data, is also shown. The fragmentation and factoriza- 
tion scales were fixed as in [128]. In the lower right panel the systematic decrease of the next-to-leading 
order K-factor is presented. Two fits to Kj^lq have been used: linear K^lq = 2.7924 — 0.0999 In s 



and quadratic Knlq = 3.8444 — 0.3234 Ins + 0.0107 In^ s in Ins. For center of mass energies up to 
1 TeV the two parameterization differ by less than 15% but this difference is seen to grow to 30%-50% 
at ^ = 5 - 10 TeV. 




Fig. 4: The predicted LO differential cross section daP^ /dycPpT for inclusive neutral pion and charged hadron 
production at midrapidity y — in p + p (p + p) reactions is shown for y/s — 17, 200, and 5500 GeV. The ratio of 
neutral pions to inclusive charged hadrons versus pt is given in the right panel. 



In Fig. 4 the predicted transverse momentum distribution of neutral pions and inclusive charged 
hadrons is shown, corresponding to the quadratic in In s fit to K^lq for energies typical of SPS, RHIC, 
and the LHC. The significant hardening of the spectra with y/s has two important consequences for p+A 
and A + A collisions: a notably reduced sensitivity to initial state kinematic effects (smaller Cronin) and 
larger variation of the manifested final-state multi-parton scattering (energy loss) with px [263]. We 
have also investigated the effect of isospin asymmetry between p + p and p + p reactions in vr*^ and 
/i^ + /i^ production and found it to be small. More quantitatively, at y/s = 5.5 TeV the fractional 
difference \daPP - daPP\/daPP varies from 2.5% at px = 5 GeV to 4.8% at pr = 150 GeV. This is 
insignificant as compared to the projected 50% uncertainty that comes from the extrapolation of Knlo 
in LO calculations (see Fig. 3) or the choice of scale in NLO calculations. A recent study showed no 
deviation from DGLAP evolution at = 10 GeV^ down to x = 10"^ in N + N reactions [129]. The 
nuclear amplification effect oc A^^^ ~ 10 for a large nucleus is still insufficient to enable measurements 
of high initial gluon density QCD at RHIC, but will play an important role at the LHC. 

2.22 Perturbative QCD Hadron Composition 

The predicted hadron composition mp + p {p + p) reactions is plotted in the right panel of Fig. 4. The 
proton+kaon fraction is seen to increase systematically with pT (xt = 'ipx / \fs) and is reflected in the 
decreasing 7r'^/0.5(/i^ + h^). At RHIC and LHC energies this ratio becomes ~ 0.5 at px — 15 GeV 
and px ~ 75 GeV, respectively. At transverse momenta ~ 2 — 4 GeV the contribution of baryons and 
kaons to /i+ + /i~ is < 20%. This is significantly smaller compared to data on + reactions, with 



the discrepancy being amplified in central A + A. Possible explanations include: enhanced baryon pro- 
duction via topological gluon configurations (junctions) and its interplay with jet quenching [271, 154] 
in A -\- A [261, 262], hydrodynamic transverse flow [254], uncertainty of the fragmentation functions 
Dp/cizcQ"^) into protons and antiprotons [294], and quark recombination driven by unorthodox (ex- 
tracted) parton distributions inside nuclei [167]. The approaches in Refs. [261, 262, 254] also address 
the centraUty dependence of the baryon/meson ratios in heavy ion collisions at RHIC. In [262] in has 
been shown that similar nuclear enhancement is expected in A, A production (as compared to kaons). 
The combined low-pr baryon enhancement and the growth of the non-pionic hadron fraction in the per- 
turbative regime may lead to an approximately constant pion to charged hadron ratio in the full measured 
Pt region at RHIC at a/s^vat = ^00 GeV. We propose that the LHC may play a critical role in resolving 
the mystery of enhanced baryon production in A + A through the significantly larger experimentally 
accessible pT range. Effects associated with baryon transport and transverse flow are not expected to 
extend beyond px = W — 15 GeV and may result in a detectable minimum of the baryon/meson ratio 
versus pr before a secondary subsequent rise. On the other hand, fragmentation functions (possibly en- 
hanced at large Zc relative to current parameterizations) are expected to exhibit a much more monotonic 
behavior. 

3. FINAL STATE EFFECTS IN DENSE AND HOT MATTER 

Bjorken argument In August 1982 J. D. Bjorken published a preprint [85] on "Energy Loss of Energetic 
Partons in Quark-Gluon Plasma: Possible Extinction of High p± Jets in Hadron-Hadron ColUsions", in 
which he discussed that high energy quarks and gluons propagating through quark-gluon plasma (QGP) 
suffer differential energy loss, and where he further pointed out that as an interesting signature events 
may be observed in which the hard collisions may occur such that one jet is escaping without absorption 
and the other is fully absorbed. 

The arguments in this work have been based on elastic scattering of high momentum partons from 

quanta in the QGP, with a resulting ("ionization") loss —dE/dz ~ a^V^. with e the energy density of 
the QGP. The loss turns out to be less than the string tension of 0(1 GeV/fm) [259] . 

However, as in QED, bremsstrahlung is another important source of energy loss [154]. Due to 
multiple (inelastic) scatterings and induced gluon radiation high momentum jets and leading large p± 
hadrons become depleted, quenched [153] or may even become extinct. In [68] it has been shown 
that a genuine pQCD process (Fig. 5) is responsible for the dominant loss: after the gluon is radiated 
off the energetic parton it suffers multiple scatterings in the medium. Indeed, further studies by [70, 
69, 288, 289, 72, 290, 291, 292, 284, 285, 286, 287, 159, 158] support this observation. For reviews, 
see [74, 164, 191]. 

3.1 Radiative Energy Loss and Medium-Induced Gluon Radiation 

3.11 Qualitative Arguments 

R. Baier, UA. Wiedemann 

After its production in a hard collision, the energetic parton radiates a gluon which both traverse 
a finite size L medium. Due to its non-abelian nature and its interaction with the medium, this gluon 
follows a zig-zag path (Fig. 5), with a mean free path A > 1/fi, which is the range of screened multiple 
gluon interactions. We estimate the medium-induced gluon radiation in two different limits, requiring 
that the gluon is emitted from the hard parton if it picks up sufficient transverse momentum to decohere 
from the partonic projectile. 

In the multiple soft scattering Umit, the average phase ip accumulated by the gluon due to multiple 
scattering is 

/ \ qL ujc 

= ( Az ) ~ —L = —. (3) 



Ej parton 




Here, the medium dependence is controlled by the transport coefficient 

q ~ n'^/X ^pj (fq^ ql da/(fq± , (4) 

where p is the density of the medium (a nucleus, or partons) and a the cross section of the gluon-medium 
interaction. According to (3), gluons are emitted from a hard parton traversing a finite path length L in 
the medium, if the phase (p > 1. Thus, their energy has to be smaller than the "characteristic gluon 
frequency" 

C0c=lqL\ (5) 

For an estimate of the shape of the energy distribution, consider the number Nco^^ of scattering centers 
which add coherently in the gluon phase (3), ~ Ncoh A*^- Based on expressions for the coherence 
time of the emitted gluon, tcoh — p- — and Nco^. = = ■\J'j^^ one estimates for the gluon 
energy spectrum per unit path length 

dl 1 d/lscatt a, , , 

'^^T^-ITt — ^ ^ J - 7 — "^asX -. (6) 
dujaz iVcoh dwaz tcoh V co 

Again, this 1/y^-energy dependence of the medium-induced non-abelian gluon energy spectrum is 
expected for sufficiently small uj < cUc- The average energy loss of the parton (in the limit Eparton —>■ oo) 
due to gluon radiation with a spectrum ^ is then determined by the characteristic gluon energy lJc as 
follows, 

LodI 

= / dco c^asCOc , (7) 

It shows a characteristic quadratic dependence on the in-medium pathlength. The medium-induced 
BDMPS gluon spectrum (valid for finite size L >> A and for soft gluon energies (jJqb = < 
UJ « Eparton oo) with the characteristic behavior: ^ ^ as , a; < is suppressed by 

^/Ncoh for UJ > ujgb with respect to the incoherent Gunion-Bertsch spectrum. For comparison with 
QED the Landau-Pomeranchuk-Migdal [195, 196, 217] (LPM) suppressed photon spectrum behaves as 

Opacity expansion: we turn now to the limiting case in which the radiation pattern results from an 
incoherent superposition of very few uqE single hard scattering processes positioned within path length 
L. Consider a hard partonic projectile which picks up a single transverse momentum p by interacting 
with a single hard scatterer. An additional gluon of energy uj decoheres from the projectile wave function 



2uj 

of the parton and the position of the scatterer. The relevant phase is 



if its typical formation time tcoh = is smaller than the typical distance L between the production point 



which indicates a suppression of gluons with energy a; larger than the characteristic gluon energy 



Uc = ^fi''L. (9) 

The gluon energy spectrum per unit path length can be estimated in terms of the coherence time fcoh and 
of the average number no L of scattering centers contributing incoherently 

dl^=^ as 

— ~ (noL) ^ =i (no-L)as — . (10) 

auj az icoh ^ 

This is the typical 1 / oj-dependence of the non-abelian gluon radiation spectrum in the absence of LPM- 
type destructive interference effects. It will result again in a quadratic L-dependence of the average 
energy loss [159]. 



3.12 Quantitative Results 
C.A. Salgado, U.A. Wiedemann 

There are several calculations of the inclusive energy distribution of medium-induced gluon radi- 
ation from Feynman multiple scattering diagrams [70, 69, 288, 289, 72, 290, 291, 292, 284, 285, 286, 
287, 159, 158]. They can be obtained as particular limiting cases of the following compact expres- 
sion [286, 287] 

= 7?r^:2— 2 2Re / dyi / dyi / du / dkr e u ^ 2 Jy, 



dy du iy=o L Jyi 2 



^2 n{Oa (r) 



(11) 



Here, denotes the transverse momentum of the emitted gluon. The limit kx = I^t] < on the 
transverse phase space allows to discuss gluon emission into a finite opening angle O, x = sin 6. For 
the full angular integrated quantity, x = 1- The radiation of hard quarks or gluons differs by the Casimir 
factor Cr = Cf or Ca, respectively. 

The two-dimensional transverse coordinates u, y and r emerge in the derivation of (11) as dis- 
tances between the positions of projectile components in the amplitude and complex conjugate ampli- 
tude. The longitudinal coordinates yi, yi integrate over the ordered longitudinal gluon emission points in 
ampUtude and complex conjugate ampUtude, which emerge in time-ordered perturbation theory. These 
internal integration variables play no role in the following discussion. They are explained in more detail 
in Ref. [286]. 

The properties of the medium enter Eq. (11) in terms of the product of the time-dependent density 
n(^) of scattering centers times the strength of a single elastic scattering <j{r). This dipole cross section 
a{r) is given in terms of the elastic high-energy cross section |a(q)p of a single scatterer, 

= ^ / (^^I'^^'^^l' (l-e^^-) . (12) 

The full expression (11) has been studied in two limiting cases: 

1. Multiple soft scattering limit 

For arbitrary many soft scattering centers, the projectile performs a Brownian motion in transverse 
momentum. This dynamical limiting case can be studied in the saddle point approximation of the 
path-integral (11), using[288, 290] 

n(0a(r)~ig(0r^ (13) 



Here, q{^) is the transport coefficient[69] which characterizes the medium-induced transverse mo- 
mentum squared {q^)med transferred to the projectile per unit path length A (for details and numer- 
ical estimates, see section 3.14). For a static medium, the transport coefficient is time-independent, 

q = . (14) 

A 

In the approximation (13), the path integral in (11) is equivalent to that of a harmonic oscillator. 

Technical details of how to evaluate (11) are given in Ref. [191]. In the multiple soft scattering 
approximation, the gluon energy distribution (11) depends on the characteristic gluon energy ooc 
and a dimensionless parameter R, 

R = UcL. (15) 

The spectrum is shown in Fig. 6. For the case of medium showing one-dimensional Bjorken 

expansion, R can be related to the initially produced gluon density [240]. In the limit i? — > oo, the 
full spectrum (11) reduces to a compact analytical expression first derived in [69], 



dl agCR 

lim uj— = in 

R-^oo aw TT 




2 r^c -2 I '^r. 



cosh ^1 sin' 

2oj V 2a; 




(16) 



In the limit of large and small gluon energies, this expression coincides with the qualitative expec- 
tations: it shows a characteristic l/-v/cJ-energy dependence for small u which is suppressed above 
the characteristic gluon frequency Uc- 



dl 2asCR ( for '^<Wc, 

lim w— ~ < 1 ^ (17) 

R-^oo du 71" I, 1^ (^j for u>Uc. 

2. Opacity Expansion 

In the opacity expansion [286, 159, 158], the path integral 



ICir{yi),yi;riyi),yi\u;) = Jvrexp (^i^r^ - ^niO^^ (r)^ (18) 

in (11) is expanded in powers of the dipole cross section. To first order, the entire medium- 
dependence comes from the interaction of the hard parton with a single static scattering center, 
multiplied by the number uqL = L/X of scattering centers along the path. Modeling the single 
scatterer by a Yukawa potential with Debye screening mass /j, one finds [241] 

dl''=' ^asCR. p , r - sin{r) 

u — ; = 2 (noL) dr 

du TT Jo 

xf^- , ' V (19) 

This result is also obtained to leading order in opacity from the reaction operator approach (for 
details, see section 3.13). The energy distribution (19) depends via the phase factor 7 = ^ on the 
characteristic gluon energy Uc in (9), and on the dimensionless quantity R = UcL. The energy 
distribution ( 1 9) is plotted in Fig.7. In the limit R ^ 00, the characteristic 1 / w-energy dependence 
of the estimate (10) is recovered for sufficiently large gluon energies u > Uc, 

dl^=^ asCR /-^ 1 r-sin{r) 

iim u; — ; = 2 [no L) j dr 7^ 

R-^00 dw TT Jo r + j 
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Fig. 6: The medium-induced gluon energy distribution tj in Fig. 7: The medium-induced gluon energy distribution a;J^ 
the multiple soft scattering approximation for different values for a hard quark in the = 1 opacity expansion, calculated 
of the kinematic constraint R = uJcL. for different values of the kinematic constraint .R. 



In both limiting cases, the multiple soft and the single hard scattering limit, the gluon energy distributions 
show similar dependencies on the gluon energy and the dimensionless "kinematic constraint" R = u>cL 
(R = u)c L). In the opacity expansion, one additional model parameter enters since one specifies both the 
average momentiun transfer jj, per scattering as well as the average number no L of scattering centers in- 
volved. One can estabUsh, however, a numerical relation between transport coefficient, Debye screening 
mass and opacity, for which both approximations lead to comparable results [241]. 

3. 13 Energy Loss in the Reaction Operator Approach 
I. Vitev 

In this section we review the finite opacity GLV approach [158, 159] to the computation of the in- 
duced gluon radiative energy loss in dense nuclear matter. This calculational framework is well equipped 
for practical apphcations [263, 261, 157] and underUes the numerical results presented in sections 3.43 
and 3.46. This section also complements the discussion and numerical results on the opacity expansion 
in section 3.12. 

We first discuss some of the important physical constraints in the case of nucleus-nucleus collisions 
that lead to the development of the Reaction Operator approach: 

• A prerequisite for the consistency of all current theoretical approaches to non-abelian jet energy 
loss is the path (or time) ordering of the exchanged gluons between the propagating jet-i-induced 
gluon system the dense nuclear matter. This approximation holds as long as the range of the 
typical scattering in the medium, R = is much smaller than the mean free path Xg, in 
which case diagrams with crossed gluon exchanges are suppressed by ~ e"''^" [70, 156] (see 
also section 3.11). This condition by itself puts a theoretical constraint on the applicability of the 
large number of scattering centers limit for the case of heavy ion reactions where the typical size 
of the medium (L) ~ 5 fm. The conclusion is reinforced by the fact that the hot and dense quark- 
gluon plasma, that is expected to be created in energetic A + A reactions, expands in the final state, 
which leads to additional growth of the mean free path. 

• The inherent dynamical nature of heavy ion reactions and the final state expansion of the system 
require careful treatment of the interference phases along the eikonal line of jet-i-gluon propagation 
that are the basis of the LPM destructive interference effect [195, 196, 217] in QCD. Symmetry 
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Fig. 8: Left panel: diagrammatic representation of tlie action of the direct insertion operator Di (single hit) at position Zi on a 
jet+gluon state described by an amplitude A. The generated kinematic modification, color factors or color rotation, symmetry 
factors and phases from the energy difference before and after the momentum exchange are explicitly shown. Right panel: 
diagrammatic representation of the action of the direct insertion operator Vi (double hit) at position Zi. Figure is adapted from 
Ref. [158]. 

arguments for the exchange gluons do not apply since their collective properties, e.g. the Debye 
screening mass, and correspondingly the transport coefficient q = fi^/Xg are strongly position 
dependent. Most of the contribution to the interference phases is accumulated in the early stages 
of jet propagation. The explicit solution in the GLV reaction operator approach for the medium 
induced bremsstrahlung spectrum keeps the exact position and momentum information in the jet 
and gluon phases and propagators, see Eq. (24). This particular formulation of the problem of 
multiple radiative parton scattering can therefore answer the important question about the rate at 
which the transition to the asymptotic large number of interactions limit might takes place [158, 
159]. 

• For phenomenological applications [263, 261, 157] the Reaction Operator solution for the radiative 

spectrum uodN^ /duo and the mean energy loss (A£') does implement finite kinematic bounds (e.g. 
H < kx < ijJpi ~ jU < w < and q?j. < s/A ) [158, 159], where the effective parton mass 
is determined by the medium properties [93, 88]. In the infrared and coUinear regions plays 
the a role similar to the mass of a heavy quark [120, 118, 119]. The finiteness of the available 
phase space is particularly important at RHIC energies, but this also holds true at LHC energies 
for pt < 20 — 50 GeV. The analytic solutions discussed in the previous sections and in the sections 
that follow relax the kinematic constraints. It is difficult to a posteriori adequately account for the 
overestimate of the available phase space even on average. Instead, the integration Umits have to 
be imposed directly in the full solution Eq. (24). 

• The dominance of the lowest order terms in the opacity expansion series, which has been demon- 
strated in [158, 159], is not unique to the gluon bremsstrahlung problem. For example, a perturba- 
tive expansion of nuclear enhanced power corrections in the number of quark-nucleon scatterings 
in DIS on nuclei was recently computed and resummed [237]. The full solution is well approxi- 
mated by the first few terms of the series for a wide parameter range. 

A powerful way to address multiple interactions of systems traversing abelian and non-abelian 
media is to decompose the complex multi-parton dynamics in a product of basic operator insertions that 
represent the interaction with a single scattering center [158, 159, 237, 236, 163], subject to the assump- 
tion that the mean free path of the propagating system significantly exceeds the range of the scattering 
potential as discussed above. Let Aij...i^_-^ (x, k, c) be the amplitude of the propagating jet-i-gluon system 
that has already undergone n — 1 scatterings where x = k'^/p'^ « w/E' is the gluon momentum fraction. 



k is its transverse momentum and c - its color matrix. When the composite system passes by a scattering 
center it can miss, which is which leaves its amphtude unchanged (1). If the system exchanges a single 
momentum with the scattering, left panel of Fig. 8, there corresponding modification of the color and 
kinematics of its amplitude is given by the direct insertion operator D and reads: 



-H) 



1 \N^{Ai-,...i^_^) 



(21) 



where B„ = H — C„ = k/k^ — (k — q„)/(k — qn)^) is the so-called Bertsch-Gunion amplitude for 
producing a gluon with transverse momentum k in an isolated single collision with scattering center n. 
The momentum transfer to the jet is q^. The notation ojn = (k — q„)^/2a; is for a gluon with energy 
u {ujQ = k^/2a;) and a„ is the color matrix in the du dimensional representation of the jet with color 
Casimir Cr. Ny = J2n7=i ^im,2 counts the number of virtual interactions in A^J^...i^_-^^. Te;(.Aj^...i^_j) 
is the elastic color factor associated with all n — 1 momentum transfers from the medium to the jet line. 
Similarly, for the case of two momentum transfers given by the virtual insertion operator V, right panel 
of Fig. 8, 



/ 1 \N„ Ca 
-(-2) 



Cr + Ca 



•^ii-'-in-i (2^) k, c) e 



i{uJo—uj„)z„ 



n-1 



(^w^i- ■ -in-lip^ 1^ Qn)[C)On]) 



(22) 



To build one power of the elastic scattering cross section, or equivalently one power of opacity 
X = fi = L / \g, two gluon exchanges at a fixed position Zn are needed. Therefore, it is easy to see that 
the basic operator unit that represents one additional scattering with the medium - the GLV Reaction 
Operator - has the form: 

Rn = bibn + Vn + V^ , (23) 

where -D„, Vn are defined in Eqs.(21,22). The full solution [158] for the medium induced gluon radiation 
off jets produced in a hard collisions inside the nuclear medium of length L and to all orders in the 
correlations between the multiple scattering centers is computed via the iterative action of the Reaction 
Operator on an initial condition given by the vacuum bremsstrahlung and averaging over the momentum 
transfers and the positions of the scattering center, respectively. 
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tion times, and Azk = Zk — Zk-i are the separations of subsequent scattering centers. The momen- 
tum transfers q^ are distributed according to a normahzed elastic scattering cross section |wj(qi)p = 
a~j^dU(.i/d^cii = 11^ /iric^ + /x^)^ for the color-screened Yukawa type and the radiative spectrum can be 
evaluated from (24) for any initial nuclear geometry with an arbitrary subsequent dynamical evolution 




Fig. 9: The radiative energy loss of a quark jet with energy Ejet = 5, 50, 500 GeV (at SPS, RHIC, LHC) is plotted as a 
function of the opacity L/Xg. for a static medium (Ag = 1 fm, jj. = 0.5 GeV). Solid curves show the first order in opacity 
results. The dashed curves show results up to second order in opacity, and two third order results are shown by solid triangles 
for SPS energies. Figure is adapted from Ref. [158]. 

of the matter density. It is this stage of the calculation, Eq. (24), at which the finite kinematic constraints 
have to be imposed for the remaining q, k and x integrals. 

As argued above, the expUcit all order solution for the double differential radiative spectrum, 
Eq. (24), provides an unambiguous way to study the convergence of the opacity series and the relative 
importance of its terms. It can be seen from Fig. 9 that at large jet energies the lowest order correlation 
between the jet production point one of the scatterings that follow gives the dominant contribution to 
the non-abelian energy loss. It also gives the quadratic dependence of AE on the size of the plasma, 
AE oc for static media [159]. For realistic plasmas higher order opacity corrections may become 
important only for large number of scatterings n > 5 and small jet energies £^ ~ 5 — 10 GeV. 

Despite the dominance of the first order in the opacity expansion [159], to improve the numerical 
accuracy for small parton energies we include corrections up to third order in The left panel of 
Fig. 10 shows the radiation intensity dl/dx with an infrared cut-off at small x = uj/E given by the 
plasmon mass ^. The dynamical expansion of the bulk soft matter is assumed to be of Bjorken type. 
The medium induced radiative energy loss is proportional to the density of the scattering centers in 
the medium and for the cases of 1+lD and 1+3D expansions it has been shown that a useful to drive 
the calculation by dN^/dy [157, 161] since the gluon rapidity density can be related to the hadron 
multiplicities and the number of participants in A + A collisions. The right panel of Fig. 10 shows 
the probability distribution of fractional energy loss e = Yln^n/E, numerically computed as in [162] 
from the gluon number distribution dN{x,E)/dx, in the Poisson approximation [162, 76, 240, 56] of 
independent gluon emission 

P(e, E) = f2 Pn{e, E) , P„+i(e, E) = C "° dx^ ^^^^^!i^P„(e - x^, E) . (25) 

n=Q n -h i Ja-g ax 

The number of radiated gluons {N^{E)) is finite and small which leads to an explicit finite no-radiation 
contribution Po{e,E) = e-<^'(-^)>5(e). 
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Fig. 10: Left panel: the radiative spectrum dl /dx computed numerically up to third order in opacity from Eq. (24) in a 
dynamical Bjorken expansion scenario with a gluon plasma rapidity density dN^ /dy = 2000. Right panel: the fractional 
energy loss probability distribution P(e = X^n^"/^) °f fractional radiative energy loss computed as in [162]. The no- 
radiation 5-function contribution is not shown. 



Analytic examples in the one scattering center limit 

In the case of 1+lD Bjorken longitudinal expansion with initial plasma density po and forma- 
tion time To, i.e. p{t) = po {tq/t), it is possible to obtain a closed form analytic formula [157] un- 
der the strong asymptotic no kinematic bounds assumption. For a hard jet penetrating the quark-gluon 
plasma the LPM effect originates from the formation physics function defined in [157] as /(x,r) = 
Jq°° du [1 — cos ( uZ{x, t))] /[u{l + u)]. With Z{x,t) = (r — tq)p?{t)/2xE being the local forma- 
tion physics parameter, two simple analytic limits apply: for x ^ Xc = Lfi'^{L)/2E, in which case the 
formation length is large compared to the size of the medium, the small Z{x,t) limit applies, leading to 
f{Z) « 7rZ/2. The interference pattern along the gluon path becomes important and accounts for the 
the non-trivial dependence of the energy loss on L. When x <C Xc, i.e. the formation length is small 
compared to the plasma thickness, one gets f{Z) w log Z. The bremsstrahlung intensity distribution 
reads: 
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where Cr = 4/3 (3) for quark (gluon) jets. In Eq. (26) At is the transverse size of the medium, e.g. 
At = ttR"^ for central nucleus-nucleus collisions. The mean energy loss (to first order in x) integrates to 

^ 9Cij7ra3 1 dN9 /2E 3 \ 

A^(i) = — — — L ln^ + - + -- - • (27) 

4 At dy \ p?L vr y 

We emphasize the Unear rather than quadratic dependence of the energy loss on the size of the medium [157] 
in the Bjorken expansion case. The logarithmic enhancement with energy comes from the Xc < x < 1 
region [159]. In the case of sufficiently large jet energies {E — > oo) this term dominates. 

In the reaction operator approach, medium induced radiative energy loss in transversely expand- 
ing plasma is discussed in [162]. To derive an analytic expression taking transverse flow into account, 
we consider an asymmetric expanding sharp elliptic density profile the surface of which is defined by 
x^{Rx + VxT~'^ + y'^{Ry + VyT)^"^ = 1. The area of this elliptic transverse profile increases with time, 
r, as At{t) = tt{Rx + VxT){Ry + VyT). A short calculation for the oc InE' term in the opacity series 



leads to 



4 Hxiiy ay — ay ji'^L 

where ax = Vx/Rx^dy = VyjRy. This expression is a central result for transversely expanding media 
and provides a simple analytic generalization that interpolates between pure Bjorken 1+lD expansion 
for small a^^yT, and 3+lD expansion at large ax^yT. In the special case of pure Bjorken (longitudinal) 
expansion with Vx = Vy = ^ Eq. (28) reduces to Eq. (27) with At = irRxRy. We also note that for a 
jet originating near the center of the medium anA fully penetrating the plasma the enhanced escape time 
due to expansion r = — vt) compensates for the 1/(1 + vtt/R) dilution factor. Therefore, in 

this isotropic case, the extra dilution due to transverse expansion has in fact little or no effect of the total 
energy loss A£^|^(6 = fm) k, A£^|^(6 = fm), modulo logarithmic factors which become sizable 
only for large vt- An important consequence of our finding is that the inclusive azimuthally averaged jet 
quenching pattern in central collisions is approximately independent of the transverse flow. 



3.14 Estimates for Cold Nuclear Matter Transport Coefficients 
F. Arleo 

The modification of high-p-r hadro- and jet production due to multiple medium-induced interac- 
tions depends on the spatial extension of the medium, and on the probability and strength of the multiple 
scatterings which the hard partons suffer. To characterize medium-modifications of high-pr jets pro- 
duced in nucleus-nucleus collisions, and to relate them to the properties of hot and dense QCD matter 
produced in the collision region, knowledge about the multiple scattering strength of cold nuclear mat- 
ter is a baseline of obvious importance. Several parameterizations, suited for different processes, have 
been proposed to characterize this strength of multiple scattering which a hard parton undergoes while 
propagating through cold nuclear matter. Here we summarize the information currently available from 
theoretical predictions, as well as from data analysis of processes in which incoming or outgoing quarks 
propagate through nuclear matter, and we comment on the relation between different parameterizations. 

There are two parameters often used for the characterization of the strength of multiple scattering. 



BDMPS transport coefficient q 

In the approach of Baier, Dokshitzer, Mueller, Peigne, and Schiff (BDMPS) [69, 70, 72], this 
transport coefficient of the medium is given by q = p'^/Xmfp- Here p is the typical transverse mo- 
mentum exchanged with the medium and Amfp denotes the parton mean free path in the medium. 
In the following, q will refer to the transport coefficient of a propagating gluon. 
The transport coefficient g is a measure of the scattering strength of the QCD medium. It is related 
to the local density of color charges. In the BDMPS framework, q is related to the elastic scattering 
cross section cr of a parton on a scattering center in the medium, see Eq. (4). For cold nuclear 
matter, it is given by the gluon density of the nucleon (see Appendix B of [69]). This allows for a 
simple expression of the gluon transport coefficient in terms of the gluon distribution xG{x) and 
the nuclear density p, 

'i = ^^^P^G{x,Q'). (29) 

Taking Og — 1/2, p 0:^ 0.16 fm~^ and xG{x) ~ 1, Baier et al. have estimated the value for 
cold nuclear matter, q = 0.045 GeV^/fm. Since also the saturation momentum Qg of gluons for 
central gluon-nucleus (radius Ra) collisions at small x is given in terms of the gluon distribution 
function, there is a linear relation between Qg and q, [77] 

Ql = 2RAq. (30) 



• LQS scale parameter A 
In the perturbative QCD approach developed by Luo, Qiu, and Sterman (LQS) [211], A denotes 
the strength of twist-4 matrix elements. These determine the strength of double parton scattering, 
see section 3.51. 

To determine a numerical estimate, LQS calculated the momentum imbalance of di-jets in photo- 
production on nuclei. This is defined as 

A (4) = {pI,+pI,) X sinA$, (31) 

where p±. denote the transverse momenta of the partons and A<I> the angle between the two corre- 
sponding jets. Their analysis assumes that the nuclear enhancement seen in the data is due to the 
rescattering of one of the produced partons (either quark or gluon). Assuming only the rescattering 
of the quark when x is not too small, they can rewrite the momentum imbalance as 

A(4) =CR7r2a,AMV3. (32) 

Comparing with the measurements (see Table 2) reported at Fermilab by the E683 collaboration, 

A{k^) = 2 X 0.2161n A - 2 x 0.216^4^/^ , (33) 

they extract A^ ~ 0.1 GeV^. Assuming moreover the non-perturbative scale A to be greater than 
Aqc"/:), they conclude A^ = 0.05 — 0.1 GcV^. The original LQS estimate is based on rescattering 
of partons in the final state. X.F. Guo [151] gave a different estimate based on the {k^) of Drell- 
Yan pairs produced in h-A collisions. This observable shows a nuclear enhancement due to the 
multiple scatterings of the incoming quark entering the nucleus, 

A{k^) = {k^)hA-{k^)hN. (34) 

In a calculation [151] to leading order in Ug but taking into account nuclear enhanced power 
corrections, this quantity was shown to be proportional to the four-parton correlation function in 
nuclei, Tq{x, A), given in the Luo, Qiu, and Sterman model by 

Tg{x,A) = X'A'/^f^{x,A), (35) 

which depends hnearly on the length ~ A^^^ covered by the hard incoming quark, and where A 
is the unknown non-perturbative scale. Assuming only one quark flavor to contribute to the DY 
process, the LO (k^) broadening (34) reads 

A(4) = (^)a^^V3. (36) 

Experimentally, the {k"^) broadening has been measured by the NAIO and E772 collaborations in 
pion and proton induced reactions on nuclei respectively, see Table 2. The value 

A^ = 0.01 GeV^ (37) 

has been extracted from a comparison of Eq. (36) with these data. This is a factor at least 5 smaller 
than the original LQS estimate. It may come from the strong interference beyond the leading order 
for the initial-state interactions [235]. 

The following discussion focuses on these two parameters. Also, we shall express numerical 
estimates in terms of the mean energy loss per unit length, —dE/dz. Other parameterizations in the 
literature can be related to them. For example, the product no C is used in the works of Zakharov [290] 
and Wiedemann [286, 287] on parton energy loss. It denotes the product of the density no of charges 
in the medium times their scattering strength C, and can be expressed in terms of the BDMPS transport 
coefficient, q = 2 uqC. 
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Data 
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Vasiliev et al. (E866) 
[260] 


Drell-Yan 
xi dependence 


E866 


p(800 GeV)- Be, Fe, W 




Johnson et al. (E772) 
[175] 


Drell-Yan 
xi and M dependence 


Bill 
E866 


K800 GeV)- D, C, Ca, Fe, W 
p(800 GeV)- Be, Fe, W 


incorr 


Arleo [55] 


Drell-Yan 
xi dependence 


NA3 


7r-(150 and 280 GeV)- p, Pt 




Guo [151] 


Drell-Yan 
kr broadening 


NAIO 
E772 


TT" (140 and 280 GeV) - D, W 
p (800 GeV) - C, Ca, Fe, W 




Gyulassy, Vitev 
[263] 


Hadroproduction 
Cronin enhancement 


bJOO 
E605 


p (400 and oOO GeV) - Be, W 


going q 


Lug, Qiu, Sterman 
[211] 


Photoproduction 
dijet momentum 
imbalance 


E683 


7 - D, Be, C, Al, Cu, Sn, Pb 
= 21 GeV 


3 

O 


Wang, Wang [280] 


DIS fragmentation 
functions 


HERMES 


e- D, N, Kr = 7.2 GeV 



Table 2: Summary of various data analysis to extract the amount of quark energy loss in nuclear matter 

Relation between BDMPS transport coefficient, LQS scale parameter and mean energy loss The 
BDMPS transport coefficient q and the LQS scale parameter A can be connected assuming the kx broad- 
ening and the dijet momentum imbalance to be directly comparable [69]. In the BDMPS framework, the 
broadening of an incoming parton (with color Cr) is given by the transport coefficient and the length it 
has traveled though the medium, 

( BDMPS = ^ <?-Z>- (38) 

In the LQS approach to jet broadening [21 1], it is given by Eq. (32). Comparing (38) and (32), a simple 
expression between q and A is found [69] 

41/3 

q = C^TT^a.-^A^ (39) 

= ^CATT^a, X (^1^ X A^ (40) 

where L = 3/4 ro A^^^ with rg ~ 1.2 fm. The transport coefficient can also be related to the mean 
energy loss per unit length —clE/dz. This parameter depends linearly on the length L of the medium 
and is therefore proportional to the atomic mass number A^^^ [72] 

dE\ 1 



.asCnqL (41) 
= \alCACRTT^ A^/"" \^ (42) 



for an outgoing parton, while 



(43) 
(44) 



in the case of partons approaching the medium. 

For a numerical comparison of estimates for —dE/dz with q and A, we use the mean energy loss 
—dE/dz in a L = 5 fm nucleus, and obtain the following relations 

1 f-^\ = f-^\ [GeV/fm] = 1.39g [GeVVH = 22.8 [GeV^] , (45) 
^ V dz J \ dz J -^^ 

where ~ 0.5 and Cr = 4/3. 





Q 


A^ 


i-dE/dz),^ 


i-dE/dz),,, 


theory 


Brodsky, Hoyer [99] 


< 0.72 


< 0.022 


< 0.5 


< 0.5 


Baier etal. [69] 


0.045 


0.0029 


0.063 


0.19 


data analysis 


Vasiliev et al. 
(E866) [260] 


< 0.24 


< 0.014 


< 0.33 


< 0.99 


Johnson et al. 
(E772) [175] 


2.0 ±0.3 ±0.4 


0.12 ±0.02 ±0.02 


2.7 ±0.4 ±0.5 


8.2 ± 1.1 ± 1.5 


Arleo [55] 


0.14 ±0.11 


0.009 ± 0.007 


0.20 ±0.15 


0.60 ± 0.45 


Guo [151] 


0.16 


0.01 


0.23 


0.69 


Gyulassy, Vitev 
[263] 


0.1 


0.0061 


0.14 


0.42 


Luo, Qiu, Sterman 
[211] 


0.82-1.6 


0.05-0.1 


1.14-2.28 


3.4-6.8 


Wang, Wang [280] 


0.12 


0.0073 


0.17 


0.5 



Table 3: Compilation of the different estimates for the magnitude of quark energy loss, given either in terms of 
q, A, or the mean energy loss per unit length for an incoming {—dE/dz)-^^^ and outgoing quark {—dE/dz)^^^. 
The correspondence between the variables has been tempted and is explained in the text. In bold are the original 
estimates given by the various groups. 

Comparison of estimates for multiple scattering in cold nuclear matter 

The results of different numerical estimates for the scattering properties of cold nuclear matter are 
tabulated in Table 3 and summarized in Fig. 3.14. In the last subsection, the origin of the estimates of 
Baier et al., Luo, Qiu and Sterman, as well as X.F. Guo have been reviewed already. Here, we explain 
how the other estimates were obtained and we give arguments for the discrepancies between different 
estimates. 
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Fig. 1 1 : Compilation of the different estimates for the magnitude of an incoming quark mean energy loss per unit length, 
{—dE/dz)^^, in a L = 5 fm nucleus (see text and Table 3). 



Theoretical arguments [99, 69] The arguments leading Baier et al. to relate the energy loss of cold 
nuclear matter to the gluon distribution in the nucleon are discussed above. Without comparison to 
experimental data, S. Brodsky and P. Hoyer had suggested an upper bound for parton energy loss [99]. 
It is based on the argument that the minimum longitudinal momentum transfer to the hard parton due 
to gluon radiation is set by the uncertainty principle, Apz x L > 1, where L is the distance between 
two scattering centers. Let AE be the energy carried away by the emitted gluon and kx its transverse 
momentum. The longitudinal momentum transfer then reads Apz ~ 2AE, leading to 

X L > 1. (46) 



2AE 

While this expression comes directly from the Heisenberg principle, a similar expression has been ex- 
plicitly derived by Brodsky and Hoyer in a simpler QED model. The maximal radiative energy loss for 
partons (note that it should apply to both quarks and gluons) now amounts to 



-f < f (47) 

dz 2 

where k'^ can be related to the typical transverse momentum which partons acquire in the medium. 
Brodsky and Hoyer mainly emphasize that this radiative energy loss is not proportional to the energy of 
the scattering particle. They point out that a previous analysis of Drell-Yan data [143] based on (x E 
violates the uncertainty principle at large E. With the estimate k^ = 0.1 GeV^ for partons traversing 
cold nuclear matter, Brodsky and Hoyer arrive at — ^ < 0.5 GeV/fm, taking into account a similar 
energy loss as (47) due to elastic scattering. Clearly, this upper bound depends on the choice of kj, and 
is violated if k'^ turns out to be larger. Thus, the analysis of Brodsky and Hoyer does not constrain 
the scattering properties of cold nuclear matter, but it constrains the energy dependence of — ^ in the 
ultra-relativistic limiting case. 

Estimates for outgoing quarks [211, 280] The original estimate for the LQS scale parameter 
A in cold nuclear matter is based on the transverse momentum broadening of di-jets measured in pho- 
toproduction on nuclei. Another analysis of the energy loss of outgoing quarks in cold nuclear matter 
was given by E. Wang and X.-N. Wang [280]. In their approach, the multiple scattering of the produced 
quarks escaping the nucleus modifies the fragmentation functions in nuclei, D{z, Q'^,A). The strength 
of this modification (here denoted C) is again related to nuclear enhanced twist-4 parton correlation 
functions. 

The HERMES collaboration measured hadron production in e-A collisions on D, N, and Kr targets 
{^/s = 7.2 GeV) as a function of the virtual photon energy u and the momentum fraction carried by 



the produced hadron, z. These measurements give a direct access to the nuclear dependence of the 
fragmentation functions. Comparing with these preliminary data, E. Wang and X.-N. Wang found a 
good agreement provided that 

Ca^ ~ 0.00065 GeV^. (48) 
They translate this quantity into a mean energy loss of 

— =0.5 GeV/fm (49) 

for a L = 5 fm nucleus. 

Estimates for incoming quarks [151, 260, 175, 174, 55] Several works have attempted to 

parametrize multiple scattering effects of cold nuclear matter from Drell-Yan measurements in hadron- 
nucleus collisions. The estimate of X.F. Guo based on the nuclear enhanced transverse momentum 
broadening of Drell-Yan pairs was reviewed. Three other groups estimated the parton energy loss from 
the xi -dependence of Drell-Yan data. 

The data analysis of M.A. Vasiliev et al. (E866) [260] is based on Drell-Yan measurements in 
800 GeV proton induced reactions on Beryllium, Iron, and Tungsten targets. The data cover a wide 
range in the momentum fraction of the projectile parton, xi, integrated over the invariant mass interval 
4 < M < 8.4 GeV. The data analysis assumes [260] that the multiple scatterings of the incoming 
(anti)quark in the nucleus shift the momentum fraction of the hard parton on average by a 



s 

where k parametrizes the strength of the energy loss. The leading-order Drell-Yan cross section 

aP^{xi) ~ fP{xi + Axi) X /^(x2) X a (51) 

is then fitted to the data with k considered as a free parameter. The effects of nuclear shadowing 
f^{x2) / f^{x2) in (51) have been taking into account using the EKS98 parameterization [126]. 

The amount of quark energy loss was found to be negligible in these data sets, with a one-standard 
deviation upper Umit k < 0.10 GeV^. Assuming the length covered by the incoming parton to be given 
by the nuclear radius, L = 3/4 ro A^^^, one can relate easily the n parameter to the BDMPS transport 
coefficient. Using (50) in (43), it is given by 

' (52) 



mprg 

with ro ci! 1.2 fm. The upper Umit extracted by this group amounts to 

q < 0.237 GeVVfm. (53) 



M. Johnson et al. (E772) [175, 174] performed a different analysis of the E772 and E866 Drell- 
Yan measurements. They raised the point that E772 Drell-Yan data have also been used in the EKS98 
analysis. In principle, the small upper bound Eq. (53) could originate from an erroneous attribution of 
the sizable nuclear dependence to shadowing effects. The authors [175, 174] then attempt to disentangle 
the effects of quark energy loss and nuclear shadowing on the basis of a light-cone formulation of the 
Drell-Yan process which allows to calculate shadowing corrections. For the E772 and E866 data sets, 
they extract 

^-^^ = 2.73 ± 0.37 ± 0.5 GeV/fm (54) 



where the errors are statistics and systematics respectively. This value is an order of magnitude larger 
than any other estimate for cold nuclear matter and depends entirely on the vaUdity of a theoretical 
Ught-cone calculation of nuclear shadowing, whose uncertainties are difficult to evaluate. 

As illustrated by the two analysis above, the poorly known shadowing corrections in the Drell- 
Yan process render the extraction of quark energy loss from Drell-Yan data difficult. F. Arleo [55] 
emphasized that Drell-Yan production in pion induced reactions at lower beam energy is much less 
sensitive to shadowing effects [55] mainly for two reasons. First, the pion beam favors the fusion of 
valence quarks for which shadowing corrections are well fixed from DIS measurements only. Moreover, 
the low beam energy probes a target momentum fraction range X2 ^ 0.1 where shadowing is known to 
be small. 

The Drell-Yan cross section has been computed in the QCD-improved parton model to leading- 
order. In these calculations, the energy loss e is modeled by the BDMPS probabihty distribution [76] 
P{e) which was determined from the medium-induced BDMPS gluon energy distribution. The cross 
section reads 

^(xi) ~ (xi + e/E^) X /^(X2) xax P{e) , (55) 

In the fit of Eq. (55) to NA3 data, the transport coefficient was considered a free parameter. It has been 
found that 

q = 0.144 ± 0.108 GeVVfm (56) 
which corresponds to a mean energy loss per unit length 

^-^^ = 0.20 ±0.15 GeV/fm (57) 

in a large (L = 5 fm) nucleus. 

An alternative way to determine the nuclear matter transport coefficient is to analyze the multiple 
scattering and associated transverse momentum broadening of incoming partons. The Cronin effect 

observed inp + A reactions relative to the Glauber-scaled p + p result [112, 251, 51] has been analyzed 
by Gyulassy and Vitev [263] in the framework of multiple initial state scatterings of partons in cold 
nuclei. Parton broadening due to random elastic scatterings is computable from [163]. The possibility of 
hard fluctuations along the projectile path leads to a power law tail of the kx distribution that enhances 
{Ak^) beyond the naive Gaussian random walk result iJ?L/\. For a high energy parton with transverse 
momentum pT produced in a p -F A reaction Q^^x ^ Pt- ^he Cronin effect is modeled by using 

{kT)pA - {kT)pp ~ ln(l + cp^/^^^) . (58) 

Calculations are consistent with the energy, ^/s = 27.4, 38.8 GeV, and pT dependence observed in 
p + W/p + Be reactions with parameters of the nuclear medium set as follows: c//i^ = 0.18 /GcV^ and 
yU^/A = 0.05 GeV^/fm. The corresponding transport coefficient and initial state mean energy loss per 
unit length using Eq. (10) are: 

q ~ 0.1 GeVVfm (59) 

and 

~ 0.14 GeV/fm (60) 

in 

for a quark jet approaching the nucleus. The uncertainties in the above estimates are correlated with the 
uncertainties in the fragmentation functions. 




3.15 The Transport Coefficient qfor a Hot and Expanding Medium 

R. Baier, U.A. Wiedemann 

Various data on DIS electron-nucleus and hadron-nucleus collisions indicate that the multiple 
scattering properties of cold nuclear matter can be described by a cold nuclear matter transport coefficient 
^Inuciear matter < 0.5 — 1 GcV^/fm. This information is compiled in the previous section 3.14. For hot 
equilibrated matter, the estimated dependence of q on the energy density e is shown in Fig. 12. For 
example, for a quark gluon plasma, the number density is translated into e as p{T) ~ T'^ ~ e^^^. A 
"smooth" increase of q with increasing e is observed, such that 

^IqGP > ^Inuclear matter • (61) 

The QCD phase transition near e ~ 1 GeV/fm^ [176] leaves no structure in the e-dependence of q. In 
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Fig. 12: Transport coefficient as a function of energy density for different media: cold, massless hot pion gas (dotted) and 
(ideal) QGP (solid curve). Figure taken from [77] 

contrast, plotting q versus temperature, one would find a sharp increase of the transport coefficient at the 
critical temperature. 

What matters in practice for jet quenching in a heavy ion collision is for how long the transport 
coefficient takes values which are significantly above the cold nuclear matter reference value. 

In nucleus-nucleus collisions at collider energies, the produced hard partons propagate through a 
rapidly expanding medium. The density of scattering centers and thus the transport coefficient q^r) is 
expected to reach a maximal value q^ around the plasma formation time tq, and then decreases with time 
r rapidly due to the strong longitudinal and - to a lesser extent - transverse expansion, 

Qir)=qd{^y- (62) 

Here, q = characterizes the static medium discussed above. The value a = I corresponds to a one- 
dimensional, boost-invariant longitudinal expansion and approximates the findings of hydrodynamical 
simulations. The formation time tq of the medium may be set by the inverse of the saturation scale 
Psat [127] and is 0.2 fm/c at RHIC and « 0.1 fm/c at LHC. Since the time difference between the 
formation of the hard parton and the formation of the medium bulk is irrelevant for the evaluation of the 
radiation spectrum (11), one can replace in (1 1) the production time tq of the parton by 0. 




Fig. 13: LHS: The medium-induced gluon energy distribution for expanding collision regions (62) with expansion parameter 
a = 0, 0.5, 1.0 and 1.5. The value of the transport coefficient qo is taken at initial time to. RHS: The same gluon radiation 
spectrum with parameters rescaled according to (63). Figure taken from [241]. 

For a dynamically evolving medium of the type (62), the path-integral (18) in (11) is the path 
integral of a 2-dimensional harmonic oscillator with time dependent (imaginary) frequency $7^ (r) = |^ 
and mass u) [71]. In this way, one can calculate the medium-induced gluon energy distribution (11) for 
a dynamically expanding medium [71]. The result is shown in Fig.l3. The radiation spectrum uj^ 
satisfies a simple scaling law which relates the radiation spectrum of a dynamically expanding collision 
region to an equivalent static scenario. The linearly weighed Une integral [240] 

_ 2 /"To+i 

q = J dr (t -tq) q{T) (63) 

~ q{L) for TO ^ 0, (64) 

z — a 

defines the transport coefficient of the equivalent static scenario. The linear weight in (63) implies that 
scattering centers which are further separated from the production point of the hard parton are more 
effective in leading to partonic energy loss. In contrast to earlier believe that parton energy loss is most 
sensitive to the hottest and densest initial stage of the collision, this impUes for a dynamical expansion 
following Bjorken scaling [a = 1 in Eq. (62)] that all time scales contribute equally to the average 
transport coefficient. This makes partonic energy loss a valuable tool for the measurement of the quark- 
gluon plasma Ufetime. 

3.16 Angular Dependence of Radiative Energy Loss 
R. Baier, U.A. Wiedemann 

Hard jets, when produced in a heavy ion collision, will be measured in a typical calorimeter 
experiment within an angular cone of opening angle ^?cone- Here we summarize what is known about the 
angular dependence of medium-induced radiative energy loss. 

Average energy loss The mean energy loss due to gluons, induced by the medium and radiated outside 
the cone, has been investigated in [73] and more recently in [287, 75, 241]. The calculations are based 



on the integrated energy loss outside an angular cone of opening angle 9, 



Fig. 14 shows numerical results for this angular dependence obtained from evaluating Eq. (11) in the 
multiple soft scattering approximation with different values for the BDMPS transport coefficient q and 
the in-medium pathlength L. Here, ^(O) does not decrease monotonously with increasing but has 
a maximum at finite jet opening angle. The reason is that the radiative energy loss outside a cone angle 
Q receives additional contributions from the Brownian k^-broadening of the standard DGLAP vacuum 
radiation, 

(ky + qrY 

This redistribution in transverse momentum space does not affect the total k^- -integrated yield = 
0), but shows up as soon as a finite cone size is chosen. Thus, strictly speaking, the total ky-integrated 
radiative energy loss ^(0 = 0) is not the upper bound for the radiative energy loss outside a finite jet 
cone angle ^(0). 
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Fig. 14: The fraction of the total radiative energy loss AE/E emitted outside a jet cone of fixed angle Q. This calculation is 
for a jet of total energy E = 100 GeV and R = L'^ , ujc ^ q L'^ . 

To sum up, there is a simple physical reason for the non-monotonic behavior of AE{6conc) as a 
function of the jet cone, namely the redistribution of the vacuum component of gluon radiation in trans- 
verse phase space. However, the size of the effect remains unclear. First, one may expect that the region 
of small 9cone is dominated by higher order QCD contributions which are not yet taken into account. 
Second, the effect shown in Fig. 14 is seen in the multiple soft scattering approximation of Eq. (11). 
However, in the single hard scattering approximation, A£'(0conc) decreases monotonously [241]. Irre- 
spective of these differences for small opening angle, it is worth emphasizing that both approximations 
agree quantitatively for ^conc > 10° [241]. 




Universality of angular dependence The integrated mean loss A£'(6'conc)> normalized to AE = 
AE{dconc = 0), is defined by 

i?(^cone) = ^^^^ (67) 



with R{9cone = 0) = 1.. In the BDMPS Umiting case, the total energy loss AE depends only on the 
characteristic gluon energy, see Eq. (7). In this limit, it was observed that the ratio R{9cone) depends on 
a single dimensionless variable, which includes all the medium dependent parameters, namely [73] 

cone ) , (68) 

where 

c\L) = ^q{L/2f. (69) 

The typical dependence on ^cone is shown in Fig. 15. When comparing hot and cold QCD matter we 
recall that for fixed in-medium pathlength L, 

COLD- 

(70) 

This is a consequence of the temperature dependence of the BDMPS transport coefficient q. R{6 cone) 
is also universal in the sense that it is the same for an energetic quark and for an energetic gluon jet. In 
a recent study[241], the ratio R{9cone) was calculated in the presence of kinematic constraints on the 
transverse phase space. Small deviations from the universal behavior (68) were observed, but for all 
practical purposes these are neghgible. 




Fig. 15: The mean energy loss radiated outside an opening cone Scone normalized to the total average energy loss. LHS: in 
the EDMS calculation [75], the ratio is universal, depending on a single dimensionless variable only. RHS: if the kinematic 
constraints in transverse phase space are taken into account, deviations from the universal behavior remain small [241]. 



3.2 Multiple Gluon Emission and Quenching Weights 

R. Baier, U.A. Wiedemann 

Irrespective of the number of additionally radiated gluons, what matters for the medium modi- 
fication of hadronic observables is how much additional energy AE is radiated off a hard parton. In 
this section, we first discuss the so called quenching weight which is the probabiUty distribution P{AE) 
of the additional medium-induced energy loss. For independent gluon emission, this probability is the 



normalized sum of the emission probabilities for an arbitrary number of n gluons which carry away a 
total energy A£':[76] 

" dl{uji) 
auji 



'■~SJ ^ It. n 



n=0 



8\^E-Yj^iY'^'^'^ ■ (71) 

In general, the quenching weight (71) has a discrete and a continuous part, [240] 

P(A£;) = po 5(A£;) + p(A£;) . (72) 

The discrete weight po emerges as a consequence of a finite mean free path. It determines the probability 
that no additional gluon is emitted due to in-medium scattering and hence no medium-induced energy 
loss occurs. 

In order to determine the discrete and continuous part of (72), it is convenient to rewrite Eq. (71) 
as a Laplace transformation [76] 

P(A£;) = I^^^V{i^)e^^^, (73) 



'P{u) = exp 



Jo duj 



(74) 



Here, the contour C runs along the imaginary axis with Rev = 0. 

For the further discussion, it is useful to treat the medium-induced gluon energy distribution Loj^ 
in Eq. (11) explicitly as the medium modification of a "vacuum" distribution [241] 

duj duj duj 

From the Laplace transform (73), one finds the total probabiUty 

/•oo 

pitot) = dE P{AE - E) P(""^) {E) . (76) 
Jo 

This probability {A.E) is normalized to unity and it is positive definite. In contrast, the medium- 
induced modification of this probability, P{AE), is a generahzed probability. It can take negative values 
for some range in AE, as long as its normalization is unity, 

/ dEP{E) =po+ dEp{E) = 1 . (77) 
Jo Jo 

We now discuss separately the properties of the discrete contribution po and the continuous one p{E). 



Discrete part of the quenching weight The discrete part of the quenching weight can be expressed in 
terms of the total gluon multiplicity, 

Po = lim V{u) = exp [-N{oj = 0)] , (78) 

u—>oo 

where the multiplicity N{u) of gluons with energy larger than u emerges by partially integrating the 
exponent of (74), 

For the limiting case of infinite in-medium pathlength, the total multiplicity N{lo) diverges and the dis- 
crete part vanishes. In general, however, po is finite. A typical dependence of po on model parameters 
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Fig. 16: The discrete part po of the quenching weight calculated in the multiple soft scattering limit as a function of R. Figure 
taken from [241]. 

is shown in Fig. 16 for the radiation spectrum calculated in the multiple soft scattering hmit. A quali- 
tatively similar behavior is found in the opacity expansion. Remarkably, po can exceed unity for some 
parameter range, since the medium modification utj^ to the radiation spectrum (75) can be negative. 
The value po > 1 then compensates a predominantly negative continuous part p{AE) and satisfies the 
normaUzation (77). It indicates a phase space region at very small transverse momentum, into which less 
gluons are emitted in the medium than in the vacuum. This effect is more pronounced for gluons than 
for quarks. 

Continuous part of the quenching weight The continuous part p{AE) of the probability distribution 
(72) is shown in Fig. 17 calculated in the multiple soft scattering limit. In the opacity expansion, it 
looks qualitatively similar. With increasing density of scattering centers (i.e. increasing R = ^qL^) the 
probability of loosing a significant energy fraction AE increases. Also, since the interaction between 
partonic projectile and medium are larger for a hard gluon than a hard quark, the energy loss is larger 
for gluons. This can be seen in Fig 17 from the larger width of p{AE) for the gluonic case. Finally, 
as expected from the normahzation condition (77), the continuous part p{AE) shows predominantly 
negative contributions for the parameter range for which the discrete weight po exceeds unity. 

In the multiple soft scattering limit and for infinite pathlength, the quenching weight was found 
to be fit very well by a two-parameter log-normal distribution [56]. Analytically, an estimate of the 
quenching weight can be obtained [76] in the limit R —>■ oo from the small-a; approximation uj^cx^ 
in the multiple scattering limit 



This approximation is known to capture [240] the rough shape of the probability distribution for large 
system size, but it has an unphy sical large e-tail with infinite first moment J dee -Pbdms^ (^) • Remarkably, 
Equation (80) provides a semi-quantitative understanding of the degree of partonic energy loss shown in 
Fig. 17. In particular, comparing for i? ^ oo the maxima of the curves in Fig. 17 for quarks and gluons, 
one finds a displacement by a factor « 5. This agrees well with the square of the relative Casimir factors 
(Ca/Cf)'^ by which the maximum of c-Fbdms'I^) changes. 
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Fig. 17: The continuous part of tlie quencliing weight (72), calculated in the multiple soft scattering limit for a hard quark 
(upper row) or hard gluon (lower row). Figure taken from [241]. 



3.3 Collisional versus Radiative Energy Loss 

/. Lokhtin, A.M. Snigirev 

The collisional energy loss due to elastic scattering with high-momentum transfer have been orig- 
inally estimated by Bjorken in [85], and recalculated later in [223, 258] taking also into account the loss 
with low-momentum transfer dominated by the interactions with plasma collective modes. The method 
for quantum field-theoretic calculating energy loss in the low exchange momentum region of the col- 
lisions (screening effect in the plasma) have been developed by Braaten and Pisarski [93, 94, 133]. It 
allows one to calculate the hard thermal loop (HTL) corrections to the propagator of the exchanged gluon 
in the Qq — > Qq and the Qg Qg processes. 

The average collisional energy loss per mean free path A can be written as [85] 

3TE/2 

dz A j ' J dt 2k ^Tu\ J dt 

Here, A = 1// d^kp{k)a, and p{k) defines the thermal density of partonic scatterers (the sum over 
which is implicit). The factor ^ denotes the energy transfer per scattering times the flux factor of the 
incident participating partons [74]. The dominant contribution to the differential cross section da/dt for 
scattering of a parton with energy E off the "thermal" partons with energy 3T E at temperature T is 



taken to be the LO perturbative scattering cross section [85] 

da ^ 27raj{t) 



(82) 



dt t2 • 

The integrated parton scattering cross section a is regularized by the Debye screening mass squared 

There are marked differences between coUisional and radiative energy loss. For collisional energy 
loss, the scattering centers act incoherently. The value AEcoU is independent of in-medium path length, 
and it depends only logarithmically on the initial parton energy. It is determined mainly by the medium 
temperature [85] 

<x a^,T^ In E/T. (83) 

The dependence of the total collisional energy loss on in-medium pathlength can be weaker than Unear 
for an expanding medium and it is linear for a static one. 

The angular-integrated radiative energy loss of a high energy projectile parton is known to dom- 
inate over the collisional energy loss by up to an order of magnitude [154, 272, 226]. However, the 
angular dependence of the lost (i.e. redistributed) energy is very different for both mechanisms. With in- 
creasing parton energy, the maximum of the angular distribution of bremsstrahlung gluons shifts towards 
the direction of the parent parton. This means that measuring the jet energy as a sum of the energies of 
final hadrons moving inside an angular cone with a given finite size 6*0 will allow the bulk of the gluon 
radiation to belong to the jet and thus the major fraction of the initial parton energy to be reconstructed. 
Therefore, the medium-induced radiation mainly softens the particle energy distributions inside the jet, 
and increases the multipUcity of secondary particles. Only to a lesser degree does it affect the total 
jet energy. It is important to notice that the coherent Landau-Pomeranchuk-Migdal radiation induces a 
strong dependence of the radiative energy loss of a jet (but not a leading parton) on the jet angular cone 
size [203, 71, 286, 287, 159]. 

On the other hand, collisional energy loss turns out to be practically independent of jet cone 
size, because the bulk of "thermal" particles knocked out of the dense matter by elastic scatterings fly 
away in almost transverse direction relative to the jet axis [203]. In fact, in relativistic kinematics, 
E mo = 3 T, in the rest system of the target with effective mass mo we get the following estimate 
for the transverse pj,* and longitudinal p^' momenta of the incident and "thermal" particles: pip ~ ^/t, 
p\ ~ t/2mo; ~ —Vi, p\^ cri E — t/2mQ. Scattering angle di of the incident parton vanishes in the 
relativistic limit, tan^j = p\^/p\ ^ Vt/E 0. The scattering angle 9t of a struck "thermal" particle 
with respect to the initial direction of the fast parton can be estimated as tan^t = p^/p^ — 2mo/\/t. 
The minimal and maximal values of tan^j are tan^^""^ 2± 2mo/fiD and tan^^"*" 2± 2mo/\/03mo£' 
respectively. It is straightforward to evaluate the average (tan 6t) as [203] 

3TE/2 

/2mo\ 6T f da 1 

Neglecting a weak as{t) dependence we obtain (tan6't) ~ |tan^™"^ ~ 4mo/3iJ,D- Substituting hd, 
we arrive at < 6t >~ 60*^ for T > 200 MeV. This value exceeds typical cone sizes 6o ~ 10°— 30*^ used to 
experimentally define hadronic jets. This means that the major part of "thermal" particles will fly outside 
the cone of the jet and thus cause the "jet energy loss". This study indicates that radiative energy loss 
indeed dominates the medium-dependence of jets for small cone opening angles ^o- However, collisional 
energy loss may have a significant contribution to jet quenching for larger cone opening angles 6o, see 
Fig. 18. 

Here, we have considered only massless partons propagating through a dense QCD-matter. Al- 
though a full description of the coherent gluon radiation from a massive color charge is still lacking, finite 
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Fig. 18: Angular dependence of the coUisional energy loss for a 100 GeV quark-initiated jet, according to Ref.[203]. In 
comparison to the radiative energy loss shown in Fig. 14, the contribution is relatively small at small angles, but cannot be 
neglected at large angles. 

quark mass effects are expected to lead to a relative suppression of medium-induced radiation of heavy 
(especially b) quarks [120]. In this case the influence of collisional energy loss on experimental observ- 
ables of "heavy quark quenching" (such as high-mass dileptons and secondary charmonium production) 
can be comparable with the effect of medium-induced radiation [203], see section 3.44. 

3.4 Observable Consequences of Radiative Energy Loss 

3.41 ET-Distributions and pT-Spectra: Quenching Weights 

R. Baier, U.A. Wiedemann 

Assume that a hard parton looses an additional energy fraction AE while escaping the collision 
region. The medium-dependence of the corresponding inclusive transverse momentum spectra can be 
characterized in terms of the quenching factor Q [76] 



where P{AE) is the quenching weight given in Eq. (71). Here, the last line is obtained by assuming a 
powerlaw fall-off of the ;jT-spectrum. The effective power n depends in general on energy and px- It 
is n ~ 7 for the kinematic range relevant for RHIC, and it is smaller for LHC. Alternatively, instead of 
the quenching factor (85), the medium modification of hadronic transverse momentum spectra is often 
characterized by a shift factor S{pt), 





(85) 



dp^ dp^ 



(86) 



which is related to Q{pt) by 




(87) 




Fig. 19: The quenching factor (85) calculated in the multiple soft scattering limit. Upper row: calculation in the — > cx)-limit 
but with a variable sharp cut-off on the infrared part of the gluon energy distribution. Lower row: the same calculation is 
insensitive to infrared contributions if the finite kinematic constraint R = uicL < (x is included. Figure taken from [241]. 

Most importantly, since the hadronic spectrum shows a strong powerlaw decrease, what matters for the 
suppression is not the average energy loss (A£^) but the least energy loss with which a hard parton is 
likely to get away. One concludes that S{pt) < {^E) and depends on transverse momentum [76]. 

Fig. 19 shows a calculation of the quenching factor (85) in the multiple soft scattering limit. A 
qualitatively similar result is obtained in the opacity expansion. In general, quenching weights increase 
monotonically with px since the medium-induced gluon radiation is independent of the total projectile 
energy for sufficiently high energies. At very low transverse momenta, the calculation based on (11) is 
not reUable and the interpretation of the medium modification of hadronic spectra in nucleus-nucleus 
collisions will require additional input (e.g. modifications due to the Cronin effect). Fig. 19 suggests, 
however, that hadronic spectra at transverse momenta pt > 10 GeV, can be suppressed significantly due 
to partonic final state rescattering. 

Finally, Fig. 19 allows to comment on the sensitivity of the perturbative calculation of coj^ on 
uncontrolled non-perturbative soft physics. The gluon energy distribution in (11) allows in principle for 
the emission of arbitrarily soft gluons. It is clear, however, that the calculation cannot be reliable in 
this soft regime. To quantify the sensitivity of the calculation to the low momentum region, Baier et 
al.[76] introduced a sharp cut-off on the i? — >^ oo gluon energy distribution which was varied between 
"^cut = and ojcut = 500 MeV. However, phase space constraints (i.e. finite R) deplete the gluon 
radiation spectrum in the soft region, see Fig. 6. As seen in Fig. 19, this decreases significantly the 
sensitivity of quenching factors to the uncontrolled infrared properties of the radiation spectrum [241]. 



3.42 Medium-Modified Fragmentation Functions 
Carlos Salgado 

In proton-proton collisions, the inclusive production of a hadron h of high enough px can be de- 
scribed by the factorization (LO) formula (1). Both the parton distribution functions and the fragmenta- 
tion functions entering this expression are obtained from global fits to experimental data. The procedure 
is well know: an initial condition containing all the non-perturbative information is evolved, by DGLAP 
equations, to scales and /jij^ respectively and then fitted, in a recursive procedure, to available data. A 
third scale, the renormaUzation scale |Lt|j, is contained in the perturbative cross section da*-'^^ through 
the running of 0! 9(/i|j). Eq. (1) leads to a fair description of the shape of high-px hadronic spectra 
while the normalization has to be adjusted by an energy-dependent ivT-factor (see e.g. [128]). Also to 
NLO, the disagreement between theory and experiment lies essentially in an albeit reduced normalization 
factor [63]. However, the theoretical K-factor a^^^ ja^*^ shows, some pT-dependence. 

For proton-nucleus or nucleus-nucleus collisions, Eq. (1) is also expected to work, though, due to 
the enhanced power corrections in the nucleus, the range of validity would be for larger than in proton- 
proton. Apart from geometrical factors, the generalization to pA or AA colUsions needs of nuclear PDF 
and medium-modified FF. This medium that modifies the parton fragmentation could be the nucleus it- 
self {cold nuclear matter) in both pA and AA collisions and/or eventually the produced high-dense state 
{hot and dense medium) in AA collisions. The nuclear PDF have been studied in several approaches and 
global fits similar to the ones for the proton are available (see the pA section in this Yellow Report for 
the state-of-the-art in the field). The case for the FF is less clear, ideally one should perform a new global 
fit for these medium-FF using modified evolution equations. These new evolution equations would take 
care of the evolution, in the medium, of a highly virtual parton to the final hadrons. Whether some- 
thing like this could be obtain for (factorized) leading twist FF, and how the evolution equations would 
be modified by the medium, is still very unclear. Finite temperature modifications to DGLAP evolu- 
tion have been calculated in [231] in the framework of a thermal field theory. The modified spUtting 
functions depend in this case on the temperature of the medium. However, multiple scattering effects 
as induced radiation or interference (LPM) are not included. There has been an attempt of constructing 
medium-modified fragmentation functions from a twist expansion in Refs. [152][278] [280][232]. In 
these references, the medium-modifications of the FF are given by one additional colUsion of the parton 
with the medium (the first term in an opacity expansion). These terms are of higher-twist nature. On the 
other hand, the medium induced radiation that could eventually lead to modified leading-twist evolution 
equations including multiple scattering effects has been computed in several approximations [70]-[159]. 
In summary, a full leading-twist DGLAP-like evolution containing all the relevant features of the prob- 
lems is still missing. Most of the present approaches using leading twist FF rely on a model proposed in 
Ref. [273]. The medium-modified fragmentation functions are, in this model, given by 

I?£f(^,Ml)= r ^ d,PE{e)-^D^^^{-^,ixl). (88) 
7o i — e i — e 

The picture is the following. The high energetic quark or gluon loses some fraction e of its energy 
when traveling through the medium and then fragments in the vacuum with the normal vacuum FF, 
Dk^yXzjX — e^i^p), with the corresponding shifted momentum fraction. Any modification of the virtu- 
ality dependence of the FF by the medium is neglected. Also, the hadronized remnants of the medium- 
induced soft radiation are neglected in the definition of (88). However, these remnants are expected to be 
soft, and their inclusion would thus amount to an additional contribution to d''^^^^ {z ^ [P'p] for z < 0.1 
say. The only ingredient needed in Eq. (88) is the probability distribution Pe{^) for a parton of energy E 
to lose a fraction e of this energy. These quenching weights are normally computed in the independent 
gluon emission approximation (71) with P{/S.E) = PE{e = /S.E/E)/E. 




Fig. 20: Medium-modified fragmentation functions for media of different densities (upper two panels). Multiplying these FF 
by (lower two panels) the position of the maximum gives the relevant z values in the integration of Eq. (1). The vacuum 
fragmentation functions are from Ref. [183]. 



Another approximation that has been taken for the quenching weight is just 

PE{e) = b[e-^. (89) 

It has been argued that using Eq. (89) produces a much stronger effect due to the rapid pT-dependence 
of the production spectrum (or equivalently to the rapid z-dependence of the FF). This is clearly the case 
for the multiple soft scattering approximation for which the quenching weight has a sharp maximum. In 
the hard scattering opacity expansion, the longer tails of the distributions makes the difference between 
using Eqs. (71) and (89) smaller. 

Medium-modified quark to pion fragmentation functions for different media are plotted in Fig. 20 
[240] [241]. They are calculated from Eq. (88) using the multiple scattering approximation for the 
quenching weights and the LO KKP [182] parametrization of (5^)- For this calculation, the 

virtuality Q of Dhiq{z, Q"^) is identified with the (transverse) initial energy Eq of the parton. This is 
justified since Eq and Q are of the same order, and Di^/q{z, Q^) has a weak logarithmic Q-dependence 
while medium-induced effects change as a function of e = ^ ~ ^(^)- ^ collision region ex- 
panding according to Bjorken scaling, the transport coefficient can be related to the initial gluon rapidity 
density [69, 157]. 



1 L2 dN9 



That's what is done in Fig. 20. Interestingly, Eq. (90) indicates how partonic energy loss changes with the 
particle multiplicity in nucleus-nucleus collisions. This allows to extrapolate parton energy loss effects 
from RHIC to LHC energies [240]. 



In principle, the medium modified fragmentation function should be convoluted with the hard par- 
tonic cross section and parton distribution functions in order to determine the medium modified hadronic 
spectrum. For illustration, however, one may exploit that hadronic cross sections weigh D^^^\z, Q^) 

by the partonic cross section dai/dp\ ~ and thus effectively test z'^^^^'^^'^D'^^f^ {z, Q^) [128]. 

The value n = 6 which characterizes [128] the power law for typical values at RHIC (-y/i = 200 GeV 
and pt ~ 10 GeV). Then, the position of the maximum Zmax of z^Z)^™^'^^ (z, Q^) corresponds to the 
most hkely energy fraction Zmax-E'q of the leading hadron. And the suppression around its maximum 
translates into a corresponding relative suppression of this contribution to the high-pT hadronic spec- 
trum at Pt ~ ZmaxEq. In general, the suppression of hadronic spectra extracted in this way is in rough 
agreement with calculations of the quenching factor (85). 

Medium-modified fragmentation functions have been applied to lepton-nucleus and nucleus-nucleus 
collisions in this framework. The idea is to compare the strength of the effect in both systems. In this way, 
information about the relative densities of the media could be obtained. In this section we summarize the 
results obtained up to now for both lepton-nucleus and AA collisions. 
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Fig. 21: Ratios of medium-induced FF over those for vacuum 
computed in the twist expansion from Ref. [280]. 

Fig. 22: Theoretical calculations with (solid lines) and without 
(dashed lines) finite formation times are computed [57] in the 
BDMPS formalism with q = 0.72GeV/fm'^. This corresponds 
to AE/L ~ 0.6 GeV/fm. (See Ref. [57] for details.) 
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In the case of cold nuclear matter, experiments of lepton-nucleus scattering measure the inclusive 
particle cross section for different nuclei. In the kinematic regime of present experiments, the valence 
quarks give the dominant contribution, so, the ratio of cross sections give direct information about the 
ratio of FF: 

' ' dah^{z,v)/dudz " Dh{z,Q^,Dy ^ 

where is the energy of the virtual photon. HERIVIES experiment has measured the ratios (91) for dif- 
ferent nuclei [35]. These data has been studied in [280] in the the twist-expansion previously mentioned 
and in [57] using the BDMPS [70] [69] gluon radiation spectrum for the energy loss. See Figs. 21 and 22 
for the comparison with data. Though the two approximations are rather different, it is interesting that 
they both result in a similar average energy loss dE/dL ~ 0.5 GeV/fm. 



i 



0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 - 
, 



7T°, Au+Au 200 GeV (PHENIX). 

Multiple soft scat. 

Single hard scat. 








Fig. 23: Suppression of 7r° production at high-pT measured by PHENIX collaboration [25]. Theoretical suppression is esti- 
mated [241] from ratios of z^Di^'''^{z, Q^) at the maximum (see Fig. 20). Solid lines are for multiple soft scattering approxi- 
mation and dashed for single hard scattering in the opacity expansion of Ref. [286]. 



In the case of nucleus-nucleus collisions, the produced medium is expected to be the main source 
of energy loss for particles produced in the central rapidity region. The new data from RHIC has been 
used to fix the amount of energy loss needed to reproduce the observed suppression of particles produced 
at large pt- The size of the effect is compatible with this jet quenching explanation. The fact that the 
transport coefficient is proportional to the density of the medium allows to estimate the effect for the 
LHC. This can be read out from the lower two panels in Fig. 20. For instance, the suppression for 10 
GeV quarks in a medium of 350 gluons per unit rapidity is similar to the one for 50 GeV quarks in a 
medium of five time larger density. In Fig. 23 the suppression using medium-modified fragmentation 
functions is compared with the experimental data from PHENIX [25]. 



3.43 Nuclear Modification Factors 
I. Vitev 

Dynamical nuclear effects mp + A and A + A reactions are detectable through the nuclear modi- 
fication ratio 

dfjf ^ , A da'fP 



Rba{pt) 



I- 



dyd'^PT dyd^pT 
dN^'^ib) ^TAA{b) daPP 



in p + A 



in A + A 



(92) 



dyd'^PT dyd'^pT 

where A and T^Aib) = / c^^rr4(r)rB(r — b) in terms of nuclear thickness functions r^(r) = 
/ dz pa{t^tz) are the corresponding Glauber scaling factors [147] of daPP. We note that in Rba{pt) 
the uncertainty associated with the Kjq factors, discussed in the previous sections, drops out. The ref- 
erence calculations that follow include shadowing/antishadowing/ElVIC-effect (here referred to as "shad- 
owing"), the Cronin effect, and the non-abelian energy loss of jets. The scale dependent nuclear PDFs 
read: foc/A{x,Q'^) = Sa/A{x,Q^) {Z/Af^/p{x,Q^) + N/Af^/n{x,Q'^)), where we take the isospin 
effects on average and the EKS'98 parameterization [126] of the shadowing functions S^/a{x, Q"^)- Ini- 
tial state multiple elastic scatterings have been discussed in [II, 235, 163]. From [163] the transverse 



momentum distribution of partons that have undergone an average x = L/X incoherent interactions in 
the medium can be evaluated exactly for any initial flux dN^{p): 




Gel d^QiJ d^p 



1 d(Tel^ dN'' 



(p - qi Qn) . 



(93) 



Numerical estimates of (93) show that for thin media with a few semi-hard scatterings the induced trans- 
verse momentum broadening exhibits a weak logarithmic enhancement with and is proportional to 
L oc A^/^. The transverse momentum transfer per unit length in cold nuclear matter is found to be 
/x^/Ag ~ 0.05 GeV^/fm [263] from comparison to low energy p + A data [112, 251, 51]. The left 
top and bottom panels of Fig. 24 show the predicted Cronin-i-shadowing effect in p + Pb collisions at 
= 8.8 TeV and central Pb + Pb at ^/s = 5.5 TeV without final state medium induced energy loss. 
The 4% (10%) enhancement of Rba slpx — 40 GeV comes from antishadowing and in not related no 
multiple initial state scattering. The observed difference between 7r° and 0.5(/i+ -I- ^~) reflects the differ- 
ent Sa{x, Q"^) for quarks and gluons. Cronin effect at the LHC results in slowing down of the decrease 
of Rba at small x as seen in the pt ^ limit. In contrast, at RHIC one finds ~ 30% enhancement in 
central d + Au reactions at y/s = 200 GeV and ~ 60% effect in central Au + Au relative to the binary 
collision scaled p + p result [263, 264]. At CERN-SPS energies of ^/s = 17 GeV the results are most 
striking, with values reaching 250% in d + Au and 400% in central Au + Au at — 4 GeV. For a 
summary of results on midrapidity Cronin effect at the LHC see [1, 12]. 

The manifestation of multiple initial state scattering and nuclear shadowing at forward and back- 
ward rapidities y = ±3 in p + Pb at the LHC (for CMS r? < 2.5) and d + ^-u at RHIC (for BRAHMS 
r) < 3) has also been studied in the framework of a fixed (or slowly varying) initial parton interaction 
strength [264]. At LHC energies at y = +3 (in the direction of the proton beam) the effect of the se- 
quential projectile interactions is again small (due to the much flatter rapidity and transverse momentum 
distributions) and is overwhelmed by shadowing, which is found to be a factor of 2-3 times larger than 
the y = result at small pr ~ few GeV and vanishes (Rba = 1) at pr — 50 GeV. As previously 
emphasized, initial state gluon showering can significantly change the low-px behavior of the hadronic 
spectra at the LHC beyond the current shadowing parameterization. At RHIC in d + Au reactions at 
= 200 GeV the nuclear modification ratio is qualitatively different. While near nucleus beam (back- 
ward y = —3) rapidity Rba ^0.9 — 1 at forward rapidities y = -1-3 the nuclear modification factor 
exhibits a much more dramatic px dependence. At small transverse momenta pT ^ 1 GeV hadron pro- 
duction is suppressed relative to the binary collision scaled p + p result, Rba < 0.8. The maximum 
Cronin enhancement Rb'a — ^-"^ (30%) is essentially the same as at midrapidity but slightly shifted 
to larger pT- We emphasize that both the suppression and enhancement regions are an integral part of 
the Cronin effect [112, 251, 51] that is understood in terms of probability conservation and momentum 
redistribution resulting from multiple initial state scattering [263, 293, 147, 11, 235, 163, 1, 12]. At 
forward (in the direction of the deuteron beam) rapidities a calculation as in [263, 264] demonstrates a 
broader Cronin enhancement region with Rba — 25% at = 5 GeV. This is understood in terms of 
the significantly steeper fall-off of the hadron spectra away from midrapidity that enhances the effect of 
the otherwise similar transverse momentum kicks. While the discussed moderate px interval [264] lies 
at the very edge of BRAHMS acceptance (at y = +3) the same quahtative picture holds aty = +2. 

The full solution for the medium induced gluon radiation off jets produced in a hard collision 
inside the nuclear medium of length L and computed to all orders in the correlations between the multiple 
scattering centers via the GLV reaction operator approach [158] has been discussed in section 3.13. At 
large jet energies the lowest order correlation between the jet production point one of the scatterings 
that follow has been shown to dominate and lead to a quadratic mean energy loss dependence on the 
size of the plasma, AE oc for static media [159]. To improve the numerical accuracy for small 
parton energies we include corrections to third order in opacity [263] from Eq. (24). The dynamical 
expansion of the bulk soft matter is assumed to be of Bjorken type. In the Poisson approximation of 




Pt [GeV] Pt [GeV] 

Fig. 24: The antishadowing and Cronin effects in p + Pb and central Pb + Pb without energy loss at the LHC 
(y/s = 5.5 and 8.8 TeV) are shown in the left top and bottom panels. The right panel demonstrates the dominance 
of final state radiative energy loss effects at the LHC with a much stronger dependence compared to RHIC. The 
possible restoration of the participant scaling through hydrodynamic-like feedback at pT ^ is also shown [263]. 



independent gluon emission [76, 162, 280, 240] the probability distribution P{e,E) of the fractional 
energy loss e = J2i ^i/E can be obtained iteratively from the single inclusive gluon radiation spectrum 
dN{x, E)/dx [162] as in Eq. (25). If a fast parton looses eE of its initial energy prior to hadronization 
its momentum fraction is modified to z* = ph/pd^ — ^) = ^c/(l — e)- The observable suppressed 
hadron differential cross section can be computed from Eq. (1) with the substitution 

deP{e,p,)^Df,/M,Q^) ■ (94) 

The nuclear modification factor Raa{pt) at the LHC is shown on the right panel of Fig. 24 and is 
completely dominated by final state interactions (see left panel). It shows a significantly stronger px 
dependence as compared to RHIC, where jet quenching was predicted to be approximately constant over 
the full measured moderate- to high-transverse momentum range [263] - the result of an interplay of 
shadowing, Cronin effect, and radiative energy loss. The variation of Raa at the LHC is a factor of 5: 
from 10-20 fold suppression at pT = 10 GeV to only a factor 2-3 suppression at px = 100 GeV. The 
reason for such a prominent variation is the hardening of the particle transverse momentum spectra (see 
Fig. 4) and the insufficient balancing action of multiple initial state scattering. In fact, the prediction from 
Fig. 24 is that the suppression in central Pb + Pb at \/sjvAf = 5.5 TeV at px ~ 40 GeV is comparable 
to the factor of 4-5 suppression currently observed at RHIC. 

The extrapolation of the LHC quenching calculations to small px results into suppression 
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Fig. 25: Comparison of the prediced Cronin effect 
from [263, 264] to the measured small enhancement of sin- 
gle inclusive neutral pion production in d + Au at \/s^jy ~ 
200 GeV. Data is from PHENIX [29]. Bottom panel: a test 
of a suggested interpretation of high pT-hadron suppres- 
sion as a result of initial state wavefunction modification, 
RdAu ~ VRauAu. Data from BRAHMS, PHOBOS and 
STAR [62, 67, 16] exclude such possibility. Figure adapted 
from [265]. 



Fig. 26: Predicted [263] suppression for tt (top panel) and 
+ (bottom panel) in Au + Au compared to PHENIX and 
STAR data [25, 14]. Similar quenching is found by BRAHMS 
and PHOBOS [66, 62]. Figure adapted from [265]. 



below participant scaling. Ivlore careful examination of the mean energy loss of partons, in particular 
for gluons radiating in nuclear matter at LHC densities, reveals sizable regions of phase space with 
Ai? > E. This indicates complete absorption of jets in nuclear matter. There is experimental evidence 
that this regime of txtrtrae, final state densities may have been achieved at RHIC [218, 170, 194]. In this 
case Eq. (94) has to be corrected to include the feedback of the radiated gluons into the system. This 
hydrodynamic-like feedback is expected to recover the Npart scaling in the soft pT region [263] - also 
illustrated on the right panel of Fig. 24. The effective initial gluon densities derived from the rapidity 
densities used in Fig. 24 are pg{RHIC) = 30 - SO/fm^ and pg{LHC) = 130 - 275/fm3. These are 
one to two orders of magnitude larger than the density of cold nuclear matter and are suggestive of a 
deconfined QCD state - the quark-gluon plasma. Interestingly, a recent study of non-equilibrium parton 
transport in central Au + Au and Ph + Pb at ^/sj^j^ = 200 GeV and ^/s = 5.5 TeV has found initial 
parton densities corresponding to the lower bound of the intervals quoted above [110]. 

At the time of the completion of the CERN Yellow Report experimental data on hadroproduction at 
RHIC ^/s = 200 AGeV in d+Au and Au + Au reactions became available for comparison to theoretical 
predictions. In the top panel of Fig. 25 the Cronin enhancement, resulting from initial state parton 
broadening [263, 236] is seen to compare qualitatively to the shape of the PHENIX vr*^ measurement [29] 
in minimum bias d + Au. Larger enhancement of h'^ + h~ production, consistent with results form 
low energy p + A data, is also shown [62, 67, 16]. The bottom panel rules out the scenario for the 
initial wavefunction origin of moderate and high-pr hadron suppression, see Fig. 26, since in this case 



RdAu ~ VRauAu- Fig- 25 compares the predicted [263] approximately constant suppression of 7r° 
and h+ + h- in ^/s = 200 AGeV Au + Au collisions at RHIC to PHENIX and STAR data [25, 
14]. The overall quenching magnitude and its centrality dependence are set by {{L)/AT)dN^ /dy oc 
A^p^rt, dN9/dy = 1150. We again emphasize that the shape of RauAu is a result of the interplay of all 
three nuclear effects: Cronin, shadowing, and jet quenching. The full numerical calculation takes into 
account the dynamical Bjorken expansion of the medium, finite kinematic bounds, higher order opacity 
corrections and approximates multiple gluon emission by a Poisson distribution [263]. The remarkably 
good agreement between the predicted nuclear modification factors and the experimental measurements 
give confidence in projecting the anticipated nuclear effects over a much wider dynamical moderate- and 
high-pT range at the LHC. 

Conclusions: In central ^4-1-74 reactions the nuclear modification factor Raa{pt) at the LHC 

is shown to be completely dominated by final state multi-parton interactions [263]. For comparison, 
at RHIC Cronin effect and nuclear shadowing also play an important role, leading to an approximately 
constant suppression ratio. At the SPS initial state multiple elastic scatterings dominate, resulting in 
a net enhancement of hadron production. At forward {y = +3) rapidities in d + Au at RHIC the 
Cronin enhancement region is predicted to be broader in comparison to the y = case. In contrast in 
p + Pb at the LHC nuclear shadowing dominates but in order to detect a sizable reduction relative to 
the binary collision scaled p + p cross section measurements at close to proton rapidity (ymax = 9-2 for 
-y/s = 8.8 TeV) are needed. 

The predicted decreasing Raa with pT at the LHC, if confirmed, may have important experimental 
consequences. Comparative large-E'y measurements of the difference in the full structure of the jet cone 
inp + p and A + A reactions may prove difficult for weak signals and large backgrounds. We emphasize 
that one of the easiest and most unambiguous approaches for detecting the non-abelian jet energy loss 
and performing jet-tomographic analysis of the properties of the hot and dense matter created in ultra- 
relativistic heavy ion reactions is through the suppression pattern of leading hadrons. Therefore these 
measurements should enter as an important part of the experimental programs at the LHC. 

3.44 Heavy Quark Energy Loss Observables 
R. Vogt 

Heavy quarks (HQ) are good probes of the QCD medium produced in heavy ion colfisions. They 
are produced perturbatively in the initial hard nucleon-nucleon collisions at timescales on the order of 
l/niQ. Their production during other stages of the evolution of the system is unlikely, except perhaps in 
the pre-equilibrium phase of the plasma, because mq 3> T. (See Refs. [244, 245, 225, 197] for some 
estimates of thermal charm production.) Thus the initially produced heavy quarks experience the full 
collision history. 

While the heavy quarks are in the medium, they can undergo energy loss by two means: elastic 
collisions with light partons in the system (coUisional) and gluon bremsstrahlung (radiative). We will 
briefly review some of the predicted results for —dE / dx of heavy quarks for both coUisional and radiative 
loss. We then show the predicted effect on the charm and bottom contributions to the dilepton continuum 
for both ALICE and CMS using a fixed value of —dE/dx. 

The coUisional energy loss of heavy quarks through processes such as Qg — > Qg and Qq Qq 
depends logarithmically on the heavy quark momentum, —dE/dx oc ln(gmax/Q'min)- Treatments of the 
coUisional loss vary with the values assumed or calculated for the cutoffs. These cutoffs are sensitive to 
the energy of the heavy quark and the temperature and strong coupling constant in the medium. Thus 
the quoted value of the energy loss is usually for a certain energy and temperature. The calculation was 
first done by Bjorken [85] who found -dE/dx « 0.2 GeV/fm for a 20 GeV quark at T = 250 MeV. 
Further work refined the calculations of the cutoffs [256, 257, 223], with similar results. Braaten and 
Thoma calculated the coUisional loss in the limits E ^ 'rn^/T and E S> run/T in the hard thermal 



loop approximation, removing the cutoff ambiguities. They obtained —dE/dx w 0.3 GeV/fm for a 20 
GeV charm quark and 0.15 GeV/fm for a 20 GeV bottom quark at T = 250 MeV [95]. 

Other models of heavy quark energy loss were presented in the context of J /ifj suppression: Could 
aproduced ccpak stay together in the medium long enough to form a J/V'? Svetitsky [252] calculated the 
effects of diffusion and drag on the cc pair in the Boltzmann approach and found a strong effect. The drag 
stopped the cc pair after traveling about 1 fm but Brownian diffusion drove the c and c apart quickly. The 
diffusion effect increased at later times. Essentially he predicted that the heavy quarks would be stopped 
and then go with the flow. His later calculations of D meson breakup and rehadronization [253] while 
moving through plasma droplets reached a similar conclusion. Koike and Matsui calculated the energy 
loss of a color dipole moving through a plasma using kinetic theory and found —dE/dx ~ 0.4 — 1.0 
GeV/fm for a 10 GeV QQ [185]. 

Thus the coUisional loss was predicted to be rather small, less than 1 GeV/fm for reasonable 
assumptions of the temperature. The loss increases with the energy and temperature. Using the hard 
thermal loop approach, Mustafa et al. found —dE/dx 1 — 2 GeV/fm for a 20 GeV quark at T = 500 
MeV [226]. 
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Fig. 27: The ratio of quenching factors Qh{p±) /Ql{p±) for charm and light quarks in hot matter with (J — 0.2 GeV (L = 5 
fm upper panel, L — 2fm lower panel). Solid lines correspond to unrestricted gluon radiation, while the dashed lines are based 
on the calculation with the cut lu > 0.5 GeV on gluon energies. From Ref. [120]. 

Radiative energy loss of light quarks has been extensively studied and is discussed in detail in 
the remainder of this chapter. The first application of radiative loss to heavy quarks was perhaps by 
Mustafa et al. [226]. They included the effects of only a single scattering/gluon emission, Qq — > Qqg or 
Qg — > Qgg. In this case, the loss grows as the square of the logarithm ln((7max/9min)> one power more 
than the coUisional loss, but is of the same order in the strong coupling constant [95]. Thus the radiative 
loss is guaranteed to be higher than the coUisional in this approximation. The heavy quark mass enters 
their expressions only in the definition of qma.x so that the mass dependence of the energy loss is rather 
weak. They found, for a 20 GeV quark at T = 500 MeV, -dE/dx « 12 GeV/fm for charm and 10 
GeV/fm for bottom. 

These large values suggested that energy loss could be quite important for heavy quarks. If true, 
there would be a strong effect on the QQ contribution to the dilepton continuum. Shuryak [246] was the 
first to consider this possibility for AA collisions. He assumed that low mass QQ pairs would be stopped 
in the medium, suppressing the dilepton contribution from these decays substantially. However, the 
stopped heavy quarks should at least expand with the medium rather than coming to rest, as discussed 
by Svetitsky [252]. Lin et al. then calculated the effects of energy loss at RHIC, including thermal 

"His drag coefficient A{p'^) is related to the energy loss per unit length through A{p'^) = {—dE/dx) jp^ . 




fluctuations, for a constant —dE/dx = 0.5 — 2 GeV/fm [200]. These results showed that the heavy 
quark contributions to the dilepton continuum would be reduced albeit not completely suppressed. In 
any case, the energy loss does not affect the total cross section. The heavy quarks are thus piled up at 
low pt and at midrapidity if stopped completely. The work by Lin et al. was extended to the LHC for 
—dE/dx = 1 GeV/fm [202]. These results are shown here. 

Other calculations of effects on the dilepton continuum have focused on higher mass lepton pairs. 
Lokhtin and Snigirev have calculated the effect of colUsional and radiative energy loss on the correlated 
bb contribution to the dilepton continuum in the CMS acceptance [206] and find a large depletion in the 
mass range 20 < M < 50 GeV. If the loss is large, the Drell-Yan and thermal dileptons could emerge 
from under the reduced hb decay contribution at large masses. Gallmeister et al. [141] have recently 
considered the amount of energy loss that the RHIC charm data can support. 

Before presenting the results of the model calculations of Ref. [202], we note that Dokshitzer and 
Kharzeev recently pointed out that soft gluon radiation from heavy quarks is suppressed at angles smaller 
than ^0 = ''t^q/E [120]. Thus bremsstrahlung is suppressed for heavy quarks relative to light quarks 
by the factor (1 + 9^/9'^)^'^, the 'dead cone' phenomenon. The radiative energy loss of heavy quarks 
could then be quite small. In fact, PHENIX sees little to no energy loss in their charm data [20]. These 
calculations and their implications are discussed in detail in Chapter [2] of this report. We show how the 
ratio of quenching factors for heavy to light quarks, Qh{pt)/Ql{pt), depends on L, the path length 
through the medium, in Fig. 27. These results are indicative of how the D/tt ratio might be modified in 
a heavy ion colUsion. 

We now turn to an illustration of how a constant 1 GeV/fm energy loss might affect the heavy quark 
contribution to the dilepton continuum in 5.5 TeV Pb+Pb collisions at the LHC [202]. It is based on the 
picture that a heavy quark with a transverse path. It, and mean-free path. A, undergoes on average S = 
Zt/A scatterings. The main model assumption is that the actual number of scatterings, n, is generated 
from the Poisson distribution, P{n, S) = e~^S"'/n\. In the rest frame of the medium, the heavy quark 
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Fig. 28: The dilepton invariant mass distributions in the ALICE acceptance. The e+e~ (a), p+At" (b) and efj, (c) channels are 
shown. The dashed and dotted curves are the DD and summed single B and BB decays respectively without energy loss. 
The solid and dot-dashed curves are the corresponding results with —dE/dx = 1 GeV/fm. The Drell-Yan rate is given by the 
dot-dot-dashed curve in (a) and (b). From Ref. [202]. 

then experiences momentum loss Ap = nX dE / dx so that its final momentum is p'x = Pt — Ap. The 
heavy quark will thermalize if is smaller than the average transverse momentum of thermalized heavy 
quarks with a temperature T. These thermalized quarks are given a random thermal momentum in the 



rest frame of the fluid and are then transformed back to the center-of-mass frame of the collision. The 
calculation assumes —dE/dx = 1 GeV/fm, A = 1 fm and T = 150 MeV. Even a small energy loss will 
suppress high pT and large invariant mass QQ pairs as long as \dEldx\ > {pt)/Ra ~ 0.4 GeV/fm in 
Pb+Pb collisions at 5.5 TeV where {pr) is the average transverse momentum of the heavy quark and Ra 
is the nuclear radius. 

The results for DD and BB decay contributions to the dilepton continuum in ALICE are shown 
in Fig. 28 for the e+e~, jU+/Lt~ and e/j, channels. The pseudorapidity cuts are |r7| < 0.9 for electrons 
and 2.5 < ry < 4 for muons. A momentum cut of > 1 GeV is used for both lepton types. Full 
azimuthal coverage is assumed for both the central barrel and forward muon arm. Note that the BB 
contribution also includes opposite sign lepton pairs from chain decays of individual B and B mesons. 
In the calculation, 540 cc pairs and 7 bb pairs were created in central Pb+Pb collisions using the MRS D-' 
[213] parton densities [202]. More recent parton densities such as MRST HO [214] give smaller charm 
cross sections, ~ 6 mb, instead of the 17 mb obtained with MRS D-. This would reduce the charm rate 
relative to the bottom rate by nearly a factor of five. The charm and bottom yields for low mass pairs are 
similar without loss but energy loss suppresses the charm yield much more strongly than bottom. The 
moderate loss assumed here still predicts a larger BB contribution to the e+e~ and continua than 

the Drell-Yan yield. 
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Fig. 29: The dilepton invariant mass distributions in tiie CMS acceptance. The dashed and dotted curves are the DD and 
summed single B and BB decays respectively without energy loss. The solid and dot-dashed curves are the corresponding 
results with -dE/dx = 1 GeV/fm. From Ref. [202]. 

The CMS muon acceptance is in the range \rj\ < 2.4 with a lepton pr cut of 3 GeV. After these 
simple cuts are applied, the results are shown in Fig. 29 for both DD and BB decays. Whereas for 
M < 15 GeV, the DD decays would dominate those of BB before the cuts, the measured BB decays 
are everywhere larger than those from charm mesons both before and after energy loss. The generally 
larger momentum of muons from B decays and the rather high momentum cut result in larger acceptance 
for BB decays. No DD decay pairs with M < 5 GeV survive the momentum cut. A factor of 50 loss 
in rate at M ~ 10 GeV is found before energy loss. A loss in rate by a factor of 100 is obtained when 
energy loss is included. The corresponding acceptance from BB decays is significantly larger, with a 
loss in rate of a factor of 8 before energy loss and ^ 15 with energy loss. Interestingly, the leptons in 
the decay chain of a single B meson are energetic enough for both to pass the momentum cut, causing 
the peak at M ~ 2 — 3 GeV. These results suggest that rather than providing an indirect measurement 
of the charm cross section, as postulated in Ref. [144], the dilepton continuum above the T family could 



instead measure the bb production cross section indirectly. A comparison with the spectrum from pp 
interactions at the same energy would then suggest the amount of energy loss, —dE/dx, of the medium. 
For a calculation of the effects of the dead cone on higher mass dileptons in CMS, see Ref. [208]. 

We have so far shown results for the dilepton continuum. However, the PHENIX measurement 
was of single leptons [20]. The single leptons are not as sensitive to the magnitude of the energy loss as 
the dilepton continuum [202]. 

Single leptons can be categorized as those from thermalized heavy quarks and those from heavy 
quarks energetic enough to escape after energy loss. The former mainly reflects the effective thermaliza- 
tion temperature while the latter can provide us with information on the energy loss. Single leptons with 
energies greater than 1 — 2 GeV are mainly from energetic heavy quarks and thus are more sensitive to 
the energy loss. Before energy loss, the single leptons from D decays are larger than those from b hadron 
decays for pT < 2.5 GeV. After energy loss, the b hadron decays dominate the spectra over all pr- 

We show the effect of energy loss on single electrons and muons within the ALICE acceptance in 
Fig. 30. A comparison of the pT distributions of single muons in the CMS acceptance from the decays 
of D and B mesons can also provide a measure of the b cross section, shown in Fig. 31. The muon pT 
distribution is clearly dominated by B decays. 
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Fig. 30: The pr spectrum of single electrons (a) and muons Fig. 31: The pr spectrum of single muons from charm and bot- 
(b) from charm and bottom decays within the ALICE accep- tom decays within the CMS acceptance. The dashed and dotted 
tance. The dashed and dotted curves are the D and B meson curves are the D and B meson decays respectively without en- 
decays respectively without energy loss. The solid and dot- ergy loss. The solid and dot-dashed curves are the correspond- 
dashed curves are the corresponding results with —dE/dx = 1 ing results with —dE/dx = 1 GeV/fm. From Ref. [202]. 
GeV/fm. From Ref. [202]. 



3.45 Medium-Modified Jet Shapes and Jet Multiplicities 
C.A. Salgado, U.A. Wiedemann 

To discuss the medium-dependence of jet shape observables, one can start from the probability 
-ftot(e) 0) that a fraction e = ^ of the total jet energy E is emitted outside the cone angle 0. Under the 
assumption that gluon emission follows an independent Poisson process (see section 3.2), this probability 
is given by 

The expression (95) takes into account the angular energy distribution of the parton fragmentation 
process in the vacuum, as well as its medium-modification ^^^^p = /@ d@ Since both 
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Fig. 32: The gluon multiplicity distribution (101) inside a cone size R = Gc, measured as a function of kt with respect to the 
jet axis. Removing gluons with energy smaller than Ecut from the distribution (dashed and dotted lines) does not affect the 
high-fct tails. Figure taken from Ref. [242]. 



contributions are additive, the total probability (95) can be written as a convolution of the vacuum and 
the medium-induced probability, 



^'tot(e, 0) = / '^ei i^vac(ei, e) P^ed(e - ei, 6) . 



(96) 



We define the vacuum contribution Pvac(e> ©) in terms of the jet shape p{r). This jet shape is measured 
in elementary {pp, pp or e+e^) collisions as the average fraction of calorimeter cell Et in a jet subcone 
of radius r = y^(Ar/)2 + (A^)^, 

= TT^Ett^. (97) 

We use the Fermilab DO parametrization [7] of p{r) which is based on jet shapes measured for average 
transverse energies of « 50 — 150 GeV. We work in the dijet center of mass where the jet width in 
pseudorapidity Ar/ and azimuth A^> is related to the gluon emission angle @ of our calculation as r = 0. 
In general, Pvac(ei 6) is a probability distribution of some width in e whose first moment determines the 
jet shape. In the presence of medium-effects, however, the vacuum part emits only a fraction ^'^^ 
of the total energy, and thus we have 

(e>vac(e) = j deePvac(e,e) = [l-/5(r = 6)] ^ . (98) 

Since we have no experimental data about the width of Pvac(e, ©), we choose 

/ E - AE \ 
P..c{e, e) = 5(e ^ [1 - p{r)]] |,=e • (99) 

The medium-modified jet shape Pmcd('^) = 1 — (e)tot(0) is then defined in terms of the average jet 
energy fraction radiated outside an angle 0. One finds 



Pmed(r) = 1 - (e)tot(e) = 1 - y deePtotie,&) 



= p{r) - 



AE{G) AE{e = 0) 
E E 



(l-p(r)) 



(100) 



where ^(6) = J dee P„ied{^: ©)• For realistic parameters, one finds that this jet shape is modified by 
a few percent only [242]. 

While the jet energy distribution is little affected by the medium, the multiplicity distribution 
inside the jet cone is expected to change significantly. This is seen from the medium-induced additional 
number of gluons with transverse momentum k± = |k|, produced within a subcone of opening angle 6c, 



In Fig. 32, this distribution is compared to the shape of the corresponding perturbative component, 
oc Y[ log{E sin 9c/ k±). The total partonic jet multiplicity is the sum of both components. For 
realistic values of medium density and in-medium pathlength, medium effects are seen to increase this 
multiplicity significantly (by a factor > 2 — 5) in particular in the high-Zc^ tails. Also, the shape and 
width of the distribution (142) changes sensitively with the scattering properties of the medium. More- 
over, since gluons must have a minimal energy u > k±/ sin @c to be emitted inside the jet cone, this 
high-A;^ tail is unaffected by "background" cuts on the soft part of the spectrum, see Fig. 32. This sug- 
gests that the measurement of the transverse momentum distribution of hadrons with respect to the jet 
axis is very sensitive to the transverse momentum broadening of the underlying parton shower and should 
be detectable above background. 

3.46 Jet Quenching and High-px Azimuthal Asymmetry 
LP. Lokhtin, A.M. Snigirev and I. Vitev 

The azimuthal anisotropy of high-pr particle production in non-central heavy ion collisions is 
among the most promising observables of partonic energy loss in an azimuthally non-symmetric volume 
of quark-gluon plasma. We discuss the implications of nuclear geometry for the models of partonic 
energy loss in the context of recent RHIC data and consequences for observation of jet quenching at the 



In order to interpret data on nuclear collisions from current experiments at the Relativistic Heavy 
Ion Collider (RHIC) and future experiments at the Large Hadron Collider (LHC), it is necessary to 
have knowledge of the initial conditions. There are large uncertainties in the estimates of the initial 
produced gluon density, Pg{To) ~ 15 — 50/fm'^ in central Au + Au at ^/s = 130,200 AGeV and 
Pgijo) ~ 100 — 400/fm^ in central Ph + Ph reactions at y/s = 5500 AGeV, since widely different 
models (e.g. see [270, 127]) seem to be roughly consistent with data [64]. It is, therefore, essential to 
check the energy dependence of the density of the produced quark-gluon plasma (QGP) with observables 
complementary to the particle multiphcity dN'^^/dy and transverse energy dEx/dy per unit rapidity. 
High-pT observables are ideally suited for this task because they provide an estimate [271] of the energy 
loss, AE, of fast partons, resulting from medium induced non-abelian radiation along their path, as 
first discussed in [154, 272] in the context of relativistic heavy ion reactions. The approximate linear 
dependence of AE' on pg is the key that enables high-pr observables to convey information about the 
initial conditions. However, AE also depends non-linearly on the size, L, of the medium [72, 159] 
and therefore differential observables which have well controlled geometric dependences are also highly 
desirable. 

A new way to probe AE in variable geometries was recently proposed in Refs. [277, 157]. 
The idea is to exploit the spatial azimuthal asymmetry of non-central nuclear collisions. The depen- 
dence of AE on the path length i^((/>) naturally results in a pattern of azimuthal asymmetry of high-pT 
hadrons which can be measured via the differential elliptic flow parameter (second Fourier coefficient), 
V2{pt) [267, 234, 228]. Before we show the sensitivity of the high-pr V2{pt > 2 GeV) to different 
initial conditions we briefly discuss the various model calculations for the "elUptic flow" coefficient V2: 




(101) 



LHC. 



1. The elliptic flow parameter V2 was first introduced in the context of relativistic hydrodynam- 
ics [228] and reflects the fact that due to the macroscopic sizes of large nuclei many aspects of 
A + A collisions are driven by nuclear geometry In non-central collisions the interaction region 
has a characteristic "almond-shaped" form as shown in Fig. 33. 




Fig. 33: The nuclear overlap region in non-central A + A collisions shows the importance of reaction geometry. 
Model calculation described here convert the spatial anisotropy illustrated above into momentum anisotropy of 
measured hadrons. 

Hydro calculations convert the elUpticity of the reaction volume into momentum space azimuthal 
asymmetry 

through the higher pressure gradient along the small axis. The elhptic flow is thus perfectly cor- 
related to the reaction plane and can be used for its determination [229]. Hydrodynamic sim- 
ulations [187, 254] typically describe well data from relativistic nucleus-nucleus collisions at 
\fsj^j^ = 200 GeV up to — 1.5 — 2 GeV and it is not unlikely that at LHC energies of 
V^NN = 5.5 TeV the region of validity of those calculations may extend to pt — 5 GeV. 

2. Initial conditions can also be mapped onto final state observable distributions by solving covariant 
Boltzmann transport equations as in cascade models (partonic, hadronic, and multi-phase). Elliptic 
flow in this approach is generated via multiple elastic scatterings. Calculations are sensitive to the 
choice of initial conditions [219] and are currently limited by statistics to ~ 6 GeV. It is 
interesting to note that they can match the high-pr behavior of the V2 but require extremely large 
initial gluon rapidity densities dN^ /dy ~ 16000 [219] and/or string melting [201]. 

3. Memory of the initial parton density, reaction geometry, and the consequent dynamical evolution 
is also retained by large transverse momentum partons (and fragmented hadrons) through their jet 
quenching pattern [277, 157]. While this approach is discussed in more detail below, it is important 
to emphasize here that at the single inclusive jet (or hadron) level the resulting high-p^ azimuthal 
asymmetry is also perfectly coupled to the reaction plane. It has been suggested that in the limit 
of very large energy loss the momentum asymmetry is driven by jet production from the boundary 
of the interaction volume [247]. 

4. Recently, a classical computation of the elUptic flow at transverse momenta > in the frame- 
work of gluon saturation models has been performed [255]. It was found that the azimuthal asym- 
metry is generated already at proper time r = 0, i.e. it is built in the coherent initial conditions. 



The resulting elliptic flow coefficient was found to vanish quickly V2{kT) oc /c^^(i?~^ — Ry"^) 
above Qs (~ 1 GeV for RHIC and ~ 1.4 — 2 for LHC energies) which is not supported by the 
current data. 

5. An approach that does not associate azimuthal asymmetry with the reaction plane has also been 
presented [190]. Both high-pT and low-pT V2 emerge as a back-to-back jet correlation bias (with 
arbitrary direction relative to the reaction geometry for every pT bin). For large transverse mo- 
menta V2 OC InpT/n suggest an easily detectable factor of 3 increase in going from = 5 GeV 
to pt = 100 GeV at LHC. The -integrated V2 oc \/Qs at LHC exhibits ~ 50% reduction rel- 
ative to RHIC. (It can also be deduced that f 2 is larger at the SPS in comparison to RHIC in this 
model.) 

The methods for f 2 analysis can be broadly divided in two categories: two-particle methods dis- 
cussed, e.g., in [269] and multi-particle methods [91, 92]. In two-particle methods the error on the 
determined V2 from non-flow (non-geometric) correlations is 0{1/ {V2M)), where M is the measured 
multiplicity. With multi-particle methods this error goes down typically to 0(1/(^2-/^^)), i.e., smaller 
by a factor of order M. Although it is not possible to completely eliminate the non-flow components 
to V2, experimental techniques based on higher oder cumulant analysis [91, 92] will be able in many 
cases to clearly distinguish between between reaction geometry generated azimuthal asymmetry and 
back-to-back jet bias. 

Parton energy loss and nuclear geometry For nucleus-nucleus collisions the co-moving plasma pro- 
duced in an A-\-B reaction at impact parameter h at formation time tq has a transverse coordinate distribu- 
tion at mid-rapidity Pg{Y, 2; = 0, tq). In studying jet production and propagation in nuclear environment 
it is not always technically possible to perform the Monte-Carlo averaging over the jet production points 
coincidentally with the simulation of parton fragmentation. It is therefore useful to separate the medium 
dependence of the mean jet energy loss as a function of the extent of the nuclear matter traversed and the 
azimuthal angle (t> relative to the reaction plane. The total energy loss is proportional to a line integral 
along the jet trajectory r(r, (f>) = r -\- v{(f>){T — tq), averaged over the distribution of the jet production 
points 

Fi,,,4,) = jSr^-^^S^^jyrr[^Y^lr + Hmr-r.)). (103) 



Ta{t) = J dz Pa{t^, z) and TAsib) = J d^r T4(r)TB(r — b) depend on the geometry. In particular, for 
a sharp uniform cylinder of radius i?cfr one readily gets T4(r) = {A/7rR'^g)9{Rcs — |r|) and Tab(O) = 
A'^/ttR'^q. We can therefore define the effective radius of the sharp cylinder equivalent to a diffuse 
Wood-Saxon geometry via F(0, <;/!))wood-Saxon = ^^(0, 0)sharp cylinder- For Au -\- Au collisions and 
a = 1 the above constraint gives Res ~ 6 fm. 

For a non- vanishing impact parameter b and jet direction v{4>), we calculate the energy loss as 

E F{0,(f)) E ^ E ' ^ ' 

where the modulation function ^) captures in the linearized approximation the h and dependence 
of the jet energy loss and also provides a rough estimate of the maximum ellipticity generated via corre- 
lations to the reaction plane. Fig. 34 shows the R{h, (p) modulation factor plotted against the azimuthal 
angle (p for impact parameters 6 = 2, 6, 10 fm. Note that R{b, (p) reflects not only the dimensions of the 
characteristic "almond-shaped" cross section of the interaction volume but also the rapidly decreasing 
initial plasma density as a function of the impact parameter. 

In order to compare to data at < 2 GeV at RHIC and < 5 GeV at LHC, one must also take 
into account the soft non-perturbative component that cannot be computed with the eikonal jet quenching 
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Fig. 34: The modulation function R{b, (f) is plotted versus for several impact parameters h = 2,6, 10 fm from 
Ref. [157]. Diffuse Wood-Saxon versus uniform sharp cylinder geometries are compared. The most drastic differ- 
ence between these geometries occurs at high impact parameters. 



formalism. The hydrodynamic elliptic flow [228] was found in [187] to have the monotonically growing 
form V2s{pt) ~ tanh(pr/(10 ± 2 GeV)) at = 200 AGeV and to be less sensitive to the initial 
conditions than the high-py jet quenching studied here. The interpolation between the low-pr relativistic 
hydrodynamics region and the high-p^ pQCD-computable region can be evaluated as in [157]. 

Fig. 35 shows the predicted pattern of high-p^ anisotropy. Note the difference between sharp 
cylinder and diffuse Wood-Saxon geometries at 6 = 7 fm approximating roughly 20-30% central events. 
While the central (6 = 0) inclusive quenching is insensitive to the density profile, non-central events 
clearly exhibit large sensitivity to the actual distribution. We conclude that 1^2 (pt > 2 GeV, h) provides 
essential complementary information about the geometry and impact parameter dependence of the initial 
conditions in ^ + A. In particular, the rate at which the v 2 coefficient decreases at high px is an indicator 
of the diffuseness of that geometry. Minimum bias STAR data at RHIC [132, 170] for pt > 6 GeV 
now seem to support the predicted [277, 157] slow decrease of V2 at large transverse momenta. Recently 
in [263] hadron suppression in Au+Au (Pb+Pb) relative to the binary scaled p-|-p result at pT — 5 GeV 
for RHIC conditions {y/sj^j^, dN^/dy) was found to be approximately equal to the quenching factor at 
LHC at a much larger transverse momentum scale pr — 50 GeV. One may thus anticipate proportionally 
large (~ 10 — 15%) azimuthal asymmetry for high pT at the LHC. 

Energy loss of jets in transversely expanding medium We next address the question of the effect of 
possibly large transverse expansion in relativistic heavy ion reactions on V2. In non-central collisions, 
the azimuthal asymmetry of the mean energy loss can be expanded in a Fourier series and characterized 
as 

Ak5)]((/>) = AE(l + 2(52(-B)cos2(/) + ---) . (105) 

It is correlated to the final measured elliptic "flow" of jets and hadrons and has been evaluated by using 
a full hydrodynamic calculation from Ref. [187]. In this case we use the parameterization eBC of [187] 
to initiaUze the system and treat gluon number as conserved current to calculate the density evolution 
needed in the fine integral Eq. (103), where it replaces the naive Bjorken (tq/t)" expansion. We average 
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over the jet formation points the density of which is given by the number of binary collisions per unit area 
as in the Woods-Saxon geometry used in Ref. [157]. We find that the azimuthal asymmetry of the energy 
loss is strongly reduced for realistic hydrodynamic flow velocities. This implies a much smaller V2 at 
high Pt than obtained in Ref. [157] where transverse expansion was not considered and poses questions 
about the observability of the effect at LHC. 

LHC-specific remarks 

There are several important aspects in which LHC and RHIC will differ significantly. We briefly 
discuss the implications of those differences for high-pr V2 measurements: 

1. Currently at RHIC at ^/s = 200 AGeV the > 2-3 GeV regime is perturbatively com- 
putable [263] (modulo uncertainties in the baryonic sector [261]). At LHC the pr region which 
is not accessible through the pQCD approach may extend to transverse momenta as high as 5- 
10 GeV. This would imply the validity of the relativistic hydrodynamics in this domain, the extent 
of which can be tested by looking for marked deviations in the growth of V2{pt), saturation, and 
turnover. 

2. Estimates of the initial gluon rapidity density at LHC vary from dN^ /dy = 2500 to dN^/dy = 
8000. This would imply very large parton energy loss, at least in some regions of phase space. 
In this case jet production for moderate transverse momenta may be limited to a small shell on 
the surface of the interaction region, leading to a constant V2{pt) purely determined by geome- 
try [247]. 

3. Since mean transverse expansion velocities at RHIC have been estimated to be on the order of vt — 
0.5 through relativistic hydrodynamics fits, it is natural to expect even larger values at LHC. This 
may lead to a significant reduction of the observed azimuthal asymmetry as discussed above. An 
important prediction of the approach put forth in [157] is that V2{pt) exhibits a slow decrease with 
increasing transverse momentum. This can be used to distinguish azimuthal anisotropy generated 
through energy loss from alternative mechanisms. 

Jet impact parameter dependence at the LHC In light of the discussion in Sec. 3.46 it is important 



to asses the feasibility of azimuthal asymmetry measurements for large-£^r jets via detailed simulations. 
The impact parameter dependence of jet rates in Pb + Pb collisions at the LHC was analyzed in [204]. 
The initial jet spectra at \/i = 5.5 TeV were generated with PYTH1A5.7 [248]. The initial distribution of 
jet pairs over impact parameter 6 of ^ + ^ collisions (without collective nuclear effects) was obtained by 
multiplying the corresponding nucleon-nucleon jet cross section, cr^^^j^, by the number of binary nucleon- 
nucleon sub-collisions [266]: 



cm 



^^) = rAA(6Xiv(v^) 



1 - ^TAA{b)a^N{V~s) 



(106) 



with the total inelastic non-diffractive nucleon-nucleon cross section cr^jy 2± 60 mb. 

The rescattering and energy loss of jets in a gluon-dominated plasma, created initially in the nu- 
clear overlap zone in Pb + Pb collisions at different impact parameters, were simulated. For details of 
this model one can refer to [204, 205]. To be specific, we treated the medium as a boost-invariant lon- 
gitudinally expanding fluid according to Bjorken's solution [86] and used the initial conditions expected 
for central Pb + Pb collisions at LHC [122, 123, 124]: formation time tq ^ 0.1 fm/c, initial temperature 
To ~ 1 GeV, gluon plasma density pg 1.95T^. For our calculations we have used the collisional 
part of the energy loss and the differential scattering cross section from [204]; the energy spectrum of 
coherent medium-induced gluon radiation was estimated using the BDMS formaUsm [72]. 

The impact parameter dependences of the initial energy density eo and the averaged over tp jet 
escape time (r^) from the dense zone are shown in Fig. 36 [204]. (r^) goes down almost linearly with 
increasing impact parameter b. On the other hand, Eq is very weakly dependent of b (6£q /eq ^ 10%) 
up to b on the order of nucleus radius Ra, and decreases rapidly only at 6 > Ra- This suggests that 
for impact parameters b < Ra, where ^ 60% of jet pairs are produced, the difference in rescattering 
intensity and energy loss is determined mainly by the different path lengths rather than the initial energy 
density. 

Fig. 37 shows dijet rates in different impact parameter bins for Ej^^ > 100 GeV and the pseudo- 
rapidity acceptance of central part of the CMS calorimeters, |77'''*| < 2.5, for three cases: (i) without 
energy loss, (ii) with collisional loss only, {in) with colUsional and radiative loss. The total impact 
parameter integrated rates are normalized to the expected number of Pb + Pb events during a two week 
LHC run, i? = 1.2 x 10^ s, assuming luminosity L = 5 x 10^^ cm~^s~^. The maximum and mean 
values of dN'^^^^^/db distribution get shifted towards the larger b, because jet quenching is much stronger 
in central collisions than in peripheral ones. Since the coherent Landau-Pomeranchuk-Migdal radiation 
induces a strong dependence of the radiative energy loss of a jet on the angular cone size [203, 71], 
the corresponding result for jets with non-zero cone size Oq is expected to be somewhere between {iii) 
{9o 0) and (ii) cases. Thus the observation of a dramatic change in the 6-dependence of dijet rates in 
heavy ion collisions as compared to what is expected from independent nucleon-nucleon reaction pattern, 
would indicate the existence of medium-induced parton rescattering. 

Of course, such kind of measurements require the adequate determination of impact parameter in 
nuclear collision with high enough accuracy. It has been shown in [113] that for the CMS experiment 
the very forward pseudorapidity region 3 < \r]\ < 5 can provide a measurement of impact parameter via 
the energy flow in the very forward (HF) CMS calorimeters with resolution 0-5 0.5 fm for central and 
semi-central Pb + Pb collisions (see details in the section on jet detection at CMS). 

Jet azimuthal anisotropy at the LHC While at RHIC the hadron azimuthal asymmetry at high-pf is 
being analyzed, at LHC energies one can hope to observe similar effects for the hadronic jet itself [205] 
due to the large inclusive cross section for hard jet production on a scale Et ~ 100 GeV. 

The anisotropy of the energy loss (A£^) goes up with increasing b, because the azimuthal asym- 
metry of the interaction volume gets stronger. However, the absolute value of the energy loss goes down 
with increasing b due to the reduced path length L (and eq at 6 > Ra, see Fig. 36). The non-uniform de- 
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Fig. 36: The impact parameter dependence of the ini- Fig. 37: The jet+jet rates for ii'y* > 100 GeV and |?7j'''| < 
tial energy density So{b)/eo{b = 0) in nuclear overlap 2.5 in different impact parameter bins: without energy 
zone (solid curve), and the average proper jet escape time loss (soUd curve), with coUisional loss (dashed curve), 
{tl) I Ra of from the dense matter (dashed curve) [204]. with coUisional and radiative loss (dotted curve) [204]. 

pendence of AE on the azimuthal angle 99 is then mapped onto the jet spectra in semi-central collisions. 
Fig. 38 from [205] shows the distribution of jets over ip for the cases with coUisional and radiative loss 
(a) and coUisional loss only (b) for 6 = 0, 6 and 10 fm. The same conditions and kinematical acceptance 
as in Fig. 37 were fulfilled. The distributions are normalized by the distributions of jets as a function of 
if in Pb + Pb collisions without energy loss. The azimuthal anisotropy becomes stronger in going from 
central to semi-central reactions, but the absolute suppression factor is reduced with increasing b. For 
jets with finite cone size one can expect the intermediate result between cases (a) and (b), because, as we 
have mentioned before, radiative loss dominates at relatively small angular sizes of the jet cone 9o{-^ 0), 
while the relative contribution of coUisional loss grows with increasing ^o- 

In non-central collisions the jet distribution over ip is approximated well by the form ^(1 + 
Bcos2ip) , where A = 0.5(iVmax + iVmin) and B = (iVmax - iVmin)/(A^max + A^min) = 2 (cos2(^). 
In the model [205] the coefficient of jet azimuthal anisotropy, = {cos 2ip^^^) event, increases almost 
Unearly with the impact parameter b and becomes maximum at 6 ~ \.2Ra. After that i?2^ drops rapidly 
with increasing b: this is the domain of impact parameter values where the effect of decreasing energy 
loss due to the reduction of the effective transverse size of the dense zone and the initial energy density 
of the medium is crucial and cannot be compensated by the stronger volume ellipticity. Anther impor- 
tant feature is that the jet azimuthal anisotropy decreases with increasing jet energy, because the energy 
dependence of medium-induced loss is rather weak (absent in the EDMS formalism and ~ In in the 
GLV formalism for the radiative part at high Et)- 

The advantage of azimuthal jet observables is that one needs to reconstruct only the direction of 
the jet, not its total energy. It can be done with high accuracy, while reconstruction of the jet energy 
is more ambiguous. However, analysis of jet production as a function of the azimuthal angle requires 
event-by-event measurement of the angular orientation of the reaction plane. The methods summarized in 
Ref. [267, 234, 227, 234] present ways for reaction plane determination. They are applicable for studying 
anisotropic particle flow in current heavy ion dedicated experiments at the SPS and RHIC, and may be 
also used at the LHC [205]. Recently a method for measuring jet azimuthal anisotropy coefficients 
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Fig. 38: The jet distribution over azimuthal angle for the cases with collisional and radiative loss (a) and coUisional 
loss only (b), E'r^^ > 100 GeV and < 2.5 [205]. The histograms (from bottom to top) correspond to the 
impact parameter values 6 = 0, 6 and 10 fm. 



without event-by-event reconstruction of the reaction plane was proposed [207]. This technique is based 
on the correlations between the azimuthal position of jet axis and the angles of hadrons not incorporated 
in the jet. The method has been generalized by taking as weights the particle momenta or the energy 
deposition in the calorimetric sectors. It was shown that the accuracy of the method improves with 
increasing multiplicity and particle (energy) flow azimuthal anisotropy, and is practically independent of 
the absolute values of azimuthal anisotropy of the jet itself. 

Conclusions The azimuthal anisotropy of high-pT hadron production in non-central heavy ion collisions 
is shown to provide a valuable experimental tool for studying both gluon bremsstrahlung in non-abelian 
media and the properties of the reaction volume such as its size, shape, initial parton (number and energy) 
rapidity densities, and their subsequent dynamical evolution. The saturation and the gradual decrease at 
large transverse momentum of the reaction geometry generated V2, predicted as a signature complemen- 
tary to jet quenching of strong radiative energy loss in a dense QCD plasma [157], seem now supported 
by preliminary data extending up to pt — 10 GeV at RHIC. 

The initial gluon densities in P6 -|- P6 reactions at i/^jviv = 5.5 TeV at the Large Hadron Collider 
are expected to be significantly higher than at RHIC, implying even stronger partonic energy loss. This 
may result in interesting novel features of jet quenching, such as moditication of the jet distribution over 
impact parameter [204] in addition to the azimuthal anisotropy of the jet spectrum. The predicted large 
cross section for hard jet production on a scale of Et ~ 100 GeV will allow for a systematic study of 
the differential nuclear geometry related aspects of jet physics at the LHC. 

3.47 Rapidity Distribution and Jets 
LP. Lokhtin, S.V. Shmatov, P.I. Zarubin 

Medium-induced parton energy loss may result in observable modifications in the rapidity distribu- 
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Fig. 39: Nomialized differential distribution of the total transverse energy dEr /dr) over pseudo-rapidity r) for 10.000 minimum 
bias Pb— Pb collisions at = 5.5A TeV for various impact parameters [113]. The two cases are included: with jet quenching 
(the top histogram) and without jet quenching (the lower histogram). 

tions of the transverse energy flow and charged multipUcity, dEx / drj, dEj, / drj, and drich/drj [243, 113]. 

Indeed, in several IVlonte Carlo simulations of ultrarelativistic heavy ion collisions, one observes 
for example the appearance of a wide bump in the pseudo-rapidity interval — 2 < < 2, which is due 
to jet quenching. Fig. 39 from [113] demonstrates the evolution of the effect in Pb— Pb collisions as a 
function of impact parameter (HIJING [270, 155] prediction for ^/s = 5.5A TeV). One can see that even 
peripheral Pb— Pb collisions show the effects of energy loss with the central enhancement still evident 
at impact parameters up to 12 fm. Since jet quenching due to final state re-interactions is effective only 
for the mid-rapidity region (where the initial energy density of minijet plasma is high enough), the very 
forward rapidity region, |r/| > 3, remains practically unchanged. A scan of collisions of different nuclear 
systems provides an additional test of jet quenching. Because smaller nuclei require a shorter transverse 
distance for the partons to traverse before escaping the system, the central enhancement due to energy 
loss decreases with system size as obvious from the comparison with and without energy loss. Although 
the effect has only been shown here for the global Et distributions dET/drj, qualitatively the same 
picture is seen when neutral or charged particle production is studied instead of Et- 

The greater the medium-induced energy loss, the more transverse energy is piled up at central r] 
values. This leads to an increase in energy density or "stopping" in the mid-rapidity region, in contradic- 
tion to the assumption of nuclear transparency. Results qualitatively similar to those shown in Fig. 39 can 
be obtained using the VENUS generator [283] or the Parton Cascade Model VNI [145, 146]. However, 
the physics of the VENUS nucleon rescattering or VNI parton rescattering modes is very different from 
that of the radiative energy loss mechanism in HIJING. This may be due to the fact that various different 



nuclear collective effects provide effective forms of "nuclear stopping". 

One can consider rapidity spectrum of jets itself or high-pr products of jet fragmentation (the latter 
case has been discussed for RHIC energy in a paper [233]). One expects an anti-correlation between 
the rapidity distribution of the hard jets and the global cIEt /d-q spectrum: in the region in which jets 
are suppressed most, the multiplicity should be the highest. Since jet quenching due to in-medium 
parton energy loss is strongest in mid-rapidity, the maximal suppression of jet rates as compared to 
what is expected from independent nucleon-nucleon interactions extrapolation can be observed at central 
rapidity, while the very forward rapidity domain remains again almost unchanged. Thus analyzing the 
correlation between rapidity distributions of global energy (particle) flow and hard jets, by scanning 
the wide rapidity region (up to |?7| < 5 under acceptance of CMS experiment [82]), might provide the 
important information about the pseudo-rapidity size of dense QCD-matter area. 

3.5 Medium Enhanced Higher Twist Effects 

3.51 Formalism of Medium Enhanced Higher Twist Effects 
R.J. Fries 

In perturbative QCD the so-called leading twist approximation is widely used to describe a large 
class of phenomena. There are quantities of non-perturbative nature which cannot be described by per- 
turbation theory, e.g. the bound states of QCD. Nevertheless it is possible to separate perturbative (short 
range) and non-perturbative (long range) physics in a scattering reaction. Factorization theorems (see 
e.g. [109]) enable us to shift non-perturbative physics into a set of well-defined, gauge-invariant (i.e. 
observable) and universal (i.e. process independent) quantities. These quantities can be expressed by 
matrix elements of parton operators between hadron states. 

It is possible to establish a hierarchy between the matrix elements in terms of an expansion in in- 
verse powers of the momentum transfer. The expansion parameter is A^/Q^, where Q is the perturbative 
hard scale and A (for massless QCD) has to be some non-perturbative scale. The leading contribution in 
this expansion is called leading twist or (in the cases relevant here) twist-2. Factorization theorems can 
strictly be proved only for certain processes and only up to a certain level of higher twist (see section 2 
ofRef.[l]). 

The leading twist contribution always consists of one hard scattering on the parton level. In the 
simple example of deep inelastic scattering, the showcase for pQCD, the hard scattering takes place 
between the virtual photon and a quark from the target. The non-perturbative part is described by a 
matrix element which encodes the process of taking one quark out of the nucleon and putting it back 
(in the complex conjugated graph). These matrix elements define the well known parton distributions 
fq ~ {qq)^ fg ~ {PP) for quarks and gluons respectively. 

In a nuclear environment, more precisely in ^4 + A collisions, the factorization theorems are still 
valid, but obviously the picture of a dominant single hard scattering process is doubtful. From the point 
of view of the twist expansion, the matrix elements which are factors in front of the expansion parameter 
A^/Q^, can be numerically larger compared to the case of the same observable in p + p collisions. This 
is clear since the matrix elements encode the non-perturbative long-range behaviour and will be sensitive 
to the size of the system. In fact parton distributions are expected to scale roughly with the mass number 
A of the nucleus, when we neglect shadowing corrections for the moment. However it can happen that 
some higher twist matrix elements scale more strongly with the nuclear size. They have to contain more 
operators of parton fields, corresponding to more partons that enter the hard scattering. On the level of 
twist-4 e.g. one has a matrix element of the form Tqg ~ (qqFF) which is a correlator of a quark and a 
gluon. When the indices of the parton fields are contracted in the right way, this matrix element scales 
with A'^/^. Generally, on the level of twist- (n + 2), there exists a set of matrix elements that scale with 
^i+n/3 xhese matrix elements are called nuclear enhanced. The reason for the additional factors of 
A^/^ is, that the different partons can come from different nucleons in the nucleus. 



These matrix elements and their nuclear enhancement explain the trivial fact that multiple scat- 
terings are important in collisions of large nuclei. Luo, Qiu and Sterman pointed out some time ago 
[210, 211, 212], that for large nuclei with ^ ^> 1 one can replace the twist expansion by an effective ex- 
pansion in }? A^l"^ IQ'^ , keeping only the nuclear enhanced nuclear effects. These correspond to multiple 
scatterings on the parton level in the nuclear collision [134]. 

On the level of twist-4 there exist calculations for jet production in lepton or photon induced 
reactions on nuclei. They deal with the transverse momentum broadening in jet production [151] and 
the cross section for dijet production [212]. The twist-4 contributions in these cases correspond to an 
additional final state interaction of the jets, more precisely a rescattering of the outgoing jet in the nuclear 
medium. No attempt was made so far to calculate these medium corrections for hadron induced jet 
production or nucleus-nucleus collisions. For jet production m-p-^-AoxA^Adi complete twist-4 
calculation would contain both initial and final state interactions of the partons. For an overview of 
calculations available in the case of p -I- ^, see Ref.[l]. 





Fig. 40: Soft hard scattering: the parton from the primary hard Fig. 41: Double hard scattering: the parton from the primary 
scattering is on the mass shell (indicated by the blob) and scat- hard scattering is on the mass shell and undergoes an interaction 
ters off a soft gluon. The soft gluon together with the parton that with a second hard parton. It has to radiate a gluon to come 
enters the primary hard scattering from the nucleus is described back to the mass shell. The two hard partons from the nucleus 
by a so called soft hard matrix element T^^ of twist-4. are described by a so called double hard matrix element T°*^ 

of twist-4. 



In the case of pure final state interactions there are two important mechanisms at the level of twist- 
4. A parton that already underwent one hard scattering and is on the mass shell afterwards can interact 
with the soft gluon field of the nucleus, see Fig. 40. The second important case is that the parton leaving 
the primary hard scattering interacts with a hard parton from the nucleus and has to radiate a gluon in 
order to fall back onto the mass shell. This medium induced radiation is shown in Fig. 41. 

The factorization formula for the cross section for the twist-4 contribution of hadron induced 
double scattering takes the form 

c7 ~ /a (8) Rab ®Tb + Ta® Hba «) Ib (107) 

where /a, /b are parton distributions for nucleus A, B respectively, describing one parton entering the 
parton cross sections Hab, Hba- Ta, Tb are nuclear enhanced twist-4 matrix elements, describing two 
partons from the respective nucleus. The generalization of double scattering (twist-4) to arbitrary nuclear 
enhanced twist is possible in situations where factorization theorems allow, see Ref.[l]. 

For quantitative estimates models for the twist-4 matrix elements have to be introduced. For soft 
hard matrix elements the effect of the soft gluon amounts to the appearance of an additional soft energy 
scale A. The dependence on the parton momentum fraction of the hard parton is taken to be the same 
as in the parton distribution of this parton. One therefore sets T^^ = X^A^/^f where / is the parton 



distribution of the hard parton normalized to one nucleon. Similarly double hard matrix elements are 
approximated by the product of the two parton distributions for both partons = C is 

normalization constant. 

3.52 Parton Energy Loss and Modified Fragmentation Functions 

E. Wang, X.N. Wang, B. Zhang 

The formalism of medium-induced higher twist effects was extended recently to the calculation 
of medium effects on fragmentation functions [152, 278]. This is of importance for jet physics at LHC 
since the energy loss of a parton cannot be observed directly. One has to resort to particle distributions 
within a jet and study the effect of parton energy loss by measuring the modification of the fragmentation 
function of the produced parton, Da^hi^^ A*^) which can be measured directly. This modification can be 
directly translated into the energy loss of the leading parton. 

Here we give an account of this approach which so far includes appUcations to eA DIS and Au+Au 
colUsions at RHIC. The main results will be seen to be consistent with calculations described in sec- 
tions 3.1. 



Parton energy loss in a nuclear medium As a first example, we consider the process eA DIS [278, 
152,295]. Here, we consider the semi-inclusive processes, e{Li) + A{p) — > e(L2) + + X, where 
Li and L2 are the four-momenta of the incoming and the outgoing leptons, and £h is the observed hadron 
momentum. The differential cross section for the semi-inclusive process can be expressed as 

'^^^'^ (108) 
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wherep= [p'^jOjOt] is the momentum per nucleon in the nucleus, q = L2 — L\ = [— (5^/2g~, g", Or] 
the momentum transfer, s = (jo + Li)^ and oem is the electromagnetic (EM) coupling constant. Lj^i^ is 
the leptonic tensor while W^u is the semi-inclusive hadronic tensor. 

In the parton model with the collinear factorization approximation, the leading-twist contribution 
to the semi-inclusive cross section can be factorized into a product of parton distributions, parton frag- 
mentation functions and the partonic cross section. Including all leading log radiative corrections, the 
lowest order contribution (0(q;°)) from a single hard 7* -|- g scattering can be written as 
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dzh 



dxf^ix, ti'j)Hj^J {x,p, q)Dg^hizh, fi^) . (109) 



Here, Hjli'J {x, p, q) is the LO hard matrix element. The momentum fraction carried by the hadron is de- 
fined as Zh = ^h/o." and = Q"^ /2p^q^ is the Bjorken variable, fx^j and ^j? are the factorization scales 
for the initial quark distributions /^(x, //j) in a nucleus and the fragmentation functions Dq^h{zh, /i^), 
respectively. 

In a nuclear medium, the propagating quark in DIS wiU experience additional scatterings with 
other partons from the nucleus. The rescatterings may induce additional gluon radiation and cause the 
leading quark to lose energy. Such induced gluon radiations will effectively give rise to additional terms 
in the evolution equation leading to the modification of the fragmentation functions in a medium. These 
are the so-called higher-twist corrections since they involve higher-twist parton matrix elements and are 
power-suppressed. We will consider those contributions that involve two-parton correlations from two 
different nucleons inside the nucleus. 

The generalized factorization is usually applied to these multiple scattering processes[210, 211, 
212]. In this approximation, the double scattering contribution to radiative correction from processes Uke 




Fig. 42: A typical diagram for quark-gluon re-scattering processes witli three possible cuts, central(C), left(L) and right(R). 



the one illustrated in Fig. 42 can be written in the following form, 



Here, T^{x, xl) twist-four parton matrix elements of the nucleus which can be expressed in terms of 
{A\-ipq{Q) 7+ F^{y2 ) F^'^{yY) il;q{y~)\A). The fractional momentum is defined as = 
z) and X = xb = Q"^ /2p^q^ is the Bjorken variable. 

Using the factorization approximation [278, 152, 210, 211, 212, 232], they can be related to the 
twist-two parton distributions of nucleons and the nucleus, 

T^g{x, XL) = —(1 - e-^i/^?^)[//(x + XL) XTf^ixT) + ft{x){xL + XT)f^{xL + Xt)U11) 
XA 

where C is a constant, xa = ^/MRa, fq{x) is the quark distribution inside a nucleus, and fg{x) is the 
gluon distribution inside a nucleon. A Gaussian distribution in the light-cone coordinates was assumed 
for the nuclear distribution, p{y^) = poexp{y~ /2i?^ ), where = \/2RaM /p^ and M is the 
nucleon mass. We should emphasize that the twist-four matrix element is proportional to 1 /xa = RaM, 
or the nuclear size [232]. 

Including the virtual corrections and the single scattering contribution, we can rewrite the semi- 
inclusive tensor in a factorized form with a nuclear modified fragmentation function, 
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where Dq^h{zh, p'^) and Dg^h{zh, fj?) are the leading-twist fragmentation functions. The modified 
splitting functions are given as 
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A-fq^gq{z,X,XL,iT) = A7q~,qg{l - Z , X , X L , ir) ■ (114) 

To further simplify the calculation, we assume xt xl <^ x. The modified parton matrix 
elements can be approximated by 

Tq^g{x,XL) « — (l-e-4M)//(x), (115) 
XA 



where C = 2CxT/f(a;T) is a coefficient which should in principle depend on and xt- Here we will 
simply take it as a constant. 

In the above matrix element, one can identify 1/xlP'^ = 2q~z{l — z)/iJp' as the formation time of 
the emitted gluons. For large formation time as compared to the nuclear size, the above matrix element 
vanishes, demonstrating a typical LPM interference effect. Additional scattering will not induce more 
gluon radiation, thus limiting the energy loss of the leading quark. 

Since the LPM interference suppresses gluon radiation whose formation time (r/ ~ Q'^/f^rP^^ 
is larger than the nuclear size MRa/P^ in our chosen frame, £y should then have a minimum value 

of 4 ~ Q'^/MRa ~ Q^/A^/^. Here M is the nucleon mass. Therefore, the leading higher-twist 
contribution proportional to agRA/'^T ~ oisR\/Q^ due to double scattering depends quadratically on 
the nuclear size Ra- 

With the assumption of the factorized form of the twist-4 nuclear parton matrices, there is only one 
free parameter C(Q^) which represents quark-gluon correlation strength inside nuclei. Once it is fixed, 
one can predict the z, energy and nuclear dependence of the medium modification of the fragmentation 
function. Shown in Fig. 21 are the nuclear modification factor of the fragmentation functions for ^'^N 
and ^^Kr targets as compared to the recent HERMES data[35, 221]. The predicted shape of the z- and 
z/-dependence agrees well [280] with the experimental data. A remarkable feature of the prediction is 
the quadratic A^/^ nuclear size dependence, which is verified for the first time by an experiment. By 
fitting the overall suppression for one nuclear target, we obtain the only parameter in our calculation, 
C(Q2) = 0.0060 GeV^ with a^iQ'^) = 0.33 at ~ 3 GeV^. 

We can quantify the modification of the fragmentation by the quark energy loss which is defined 
as the momentum fraction carried by the radiated gluon, 

{AZg){xB,n'^) = J J dz^zA-fq^gq{z,XB,XL,£T) 

Nc xaQ^ Jo z[1 - z) Jo xl 

where x^ = /2p'^ z{l — z) = a;s/z(l — z) if we choose the factorization scale as /j,^ = Q^. When 
xyi <C <C 1 we can estimate the leading quark energy loss roughly as 

(A.,)(.«,,') . a^^6^1„^. (117) 

Since xa = ^/MRa, the energy loss (Azg) thus depends quadratically on the nuclear size. 

In the rest frame of the nucleus, p"*" = ttin, = v, and = Q"^ /2p^q^ = Q"^ /2miqi'. One 
can get the averaged total energy loss as AE = v{Azg) ^ C{Q'^)al{Q'^)mNR\{CA/Nc)^hi{l/2xB)- 
With the determined value of C, {xb) ~ 0.124 in the HERMES experiment[35, 221] and the average 
distance (La) = Ra \/2 / vt for the assumed Gaussian nuclear distribution, one gets the quark energy 
loss dE/dL K, 0.5 GeV/fm inside a Au nucleus (see section 3.14 for a comparison to other cold nuclear 
matter estimates). 



Energy Loss and Jet Quenching in Hot Medium at RHIC To extend our study of modified frag- 
mentation functions to jets in heavy-ion collisions and to relate to results obtained in the opacity ex- 
pansion approach, we can assume {k"^) k, h"^ (the Debye screening mass) and a gluon density profile 
P(y) = {to/t)9{Ra — y)po for a 1 -dimensional expanding system. Since the initial jet production 
rate is independent of the final gluon density which can be related to the parton-gluon scattering cross 
section[69] [asXTG{xT) ^ iJ^cFgl, one has then 
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where r/ = 2Ez{l — z)/i\ is the gluon formation time. One can recover the form of energy loss in a 
thin plasma obtained in the opacity expansion approach[157], 
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Keeping only the dominant contribution and assuming ag « Ca^iraj/ii iCa=l for qg and 9/4 for gg 
scattering), one obtains the averaged energy loss, 
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Neglecting the logarithmic dependence on r, the averaged energy loss in a 1 -dimensional expanding 
system can be expressed as (^)irf ~ {dEo/dL){2To/RA), where dE^/dL oc pqRa is the energy loss 
in a static medium with the same gluon density as in a 1-d expanding system at time tq. Because 
of the expansion, the averaged energy loss {dE/dL)id is suppressed as compared to the static case and 
does not depend linearly on the system size. 

In order to calculate the effects of parton energy loss on the attenuation pattern of high pT partons 
in nuclear collisions, we use a simpler effective modified fragmentation function[273, 274], 
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where 2;^, Zg are the rescaled momentum fractions. The first term is the fragmentation function of the 
jet c after losing energy A£^c(Pc, 4>) due to medium induced gluon radiation. The second term is the 
feedback due to the fragmentation of the Ng{pc, 4>) = (AL/A) radiated gluons. This effective model is 
found to reproduce the pQCD result from Eq. (112) very well, but only when Az = AEc/E is set to be 
Az « 0.6{zg). Therefore the actual averaged parton energy loss should be AE/E = l.QAz with Az 
extracted from the effective model. The factor 1 .6 is mainly caused by the unitarity correction effect in 
the pQCD calculation. 

Since gluons are bosons, there should also be stimulated gluon emission and absorption by the 
propagating parton because of the presence of thermal gluons in the hot medium. Such detailed balance 
is crucial for parton thermalization and should also be important for calculating the energy loss of an 
energetic parton in a hot medium[279]. Taking into account such detailed balance in gluon emission, 
one can then get the asymptotic behavior of the effective energy loss in the opacity expansion framework 
[279], 
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where the first term is from the induced bremsstralung and the second term is due to gluon absorption in 
detailed balance which effectively reduce the total parton energy loss in the medium. 

Shown in Fig. 43 are numerical results of the ratios of the calculated radiative energy loss with and 
without stimulated emission and thermal absorption as functions of E/p for L/Xg = 3,5 and ag = 0.3. 
Shown in the inserted box are the energy gain via gluon absorption with (A^^^^) and without (AE^^) 
rescattering. For partons with very high energy the effect of the gluon absorption is small and can be 
neglected. However, the thermal absorption reduces the effective parton energy loss by about 30-10% for 
intermediate values of parton energy. This will increase the energy dependence of the effective parton 
energy loss in the intermediate energy region. One can parameterize such energy dependence as. 
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Fig. 43: The ratio of effective parton energy loss with (AE — A_E^'^'^ + AiJ^J^'^ + AE^]^^) and without (AE^''^^) absorption as 
a function of E/fi. Inserted box: energy gain via absorption with (AiJ^j^^^) and without (Ai?^°'^) rescattering. 



The threshold is the consequence of gluon absorption that competes with radiation that effectively shuts 
off the energy loss. The parameter fj, is set to be 1 GeV in the calculation. 

To calculate the modified high px spectra in yl + ^ collisions, we use a LO pQCD model [276, 

281], 
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with medium modified fragmentation funcitons D'^^i^ given by Eq. 121 and the fragmentation functions 
in free space D°y^(zc, Q^) are given by the BBK parameterization [84]. Here, Zc = Pt/ptc, U = Uc, 
a{ab — > cd) are elementary parton scattering cross sections and tA{b) is the nuclear thickness function 
normalized to f d'^btA{b) = A. We will use a hard-sphere model of nuclear distribution in this paper. 
The K !^ 1.5 — 2 factor is used to account for higher order pQCD corrections. The parton distributions 
per nucleon fa/Ai^a, Q'^,r) inside the nucleus are assumed to be factorizable into the parton distributions 
in a free nucleon given by the MRS D— ' parameterization and the impact-parameter dependent nuclear 
modification factor which will given by the new HIJING parameterization [199]. The initial transverse 
momentum distribution qaQ^t, , b) is assumed to have a Gaussian form with a width that includes 
both an intrinsic part in a nucleon and nuclear broadening. This model has been fitted to the nuclear 
modification of the pT spectra in p + ^4 collisions at up to the Fermilab energy y/s = 40 GeV [276]. 
The initial multiple scattering in nuclei can give some moderate Cronin enhancement of the high pT 
spectra. Therefore, any suppression of the high pT spectra in Au + Au collisions has to be caused by jet 
quenching. 

We assume a 1-dimensional expanding medium with a gluon density /?g(r, r) that is proportional to 
the transverse profile of participant nucleons. According to Eq. 120, we will calculate impact-parameter 
dependence of the energy loss as 
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where AL{b, r, (p) is the distance a jet, produced at f, has to travel along n at an azimuthal angle cj) relative 
to the reaction plane in a collision with impact-parameter b. Here, po is the averaged initial gluon density 
at To in a central coUision and {dE/dL)id is the average parton energy loss over a distance Ra in a 1-d 
expanding medium with an initial uniform gluon density po- The corresponding energy loss in a static 
medium with a uniform gluon density po over a distance Ra is [280] dEo/dL = (RA/2To){dE/dL)id. 
We will use the parameterization in Eq. (123) for the effective energy dependence of the parton quark 
energy loss. 

Shown in Fig. 44 are the calculated nuclear modification factors Rab{pt) = da^^/ (Ai'binary)c^fpp 
for hadron spectra (|y| < 0.5) in Au + Au collisions at ^ = 200 GeV, as compared to experimental 
data [23, 180, 18, 25]. Here, (iVbinary) = / d^hd'^rtA{r)tA{\h - r\). To fit the observed vr" suppres- 
sion (solid lines) in the most central colUsions, we have used jjl = l.b GeV, eo = 1.07 GeV/fm and 
Ao = l/(cjpo) = 0.3 fm. The hatched area (also in other figures in this paper) indicates a variation of 
eo = ±0.3 GeV/fm. The hatched boxes around Rab = 1 represent experimental errors in overall nor- 
malization. Nuclear kr broadening and parton shadowing together give a slight enhancement of hadron 
spectra at intermediate = 2 — 4 GeV/c without parton energy loss. 

The flat pt dependence of the tx^ suppression is a consequence of the strong energy dependence 
of the parton energy loss. The slight rise of Rab at < 4 GeV/c in the calculation is due to the 
detailed balance effect in the effective parton energy loss. In this region, one expects the fragmentation 
picture to gradually lose its validity and is taken over by other non-perturbative effects, especially for 
kaons and baryons. As a consequence, the {K + p)/'K ratio in central Au + Au colHsions is significant 
larger than in peripheral Au + Au ox p + p collisions. To take into account this effect, we add a nuclear 
dependent (proportional to (iVbinary)) soft component to kaon and baryon fragmentation functions so 
that [K + p)/'K 2 at ~ 3 GeV/c in the most central Au + Au collisions and approaches its 
p + p value at > 5 GeV/c. The resultant suppression for total charged hadrons (dot-dashed) and the 
centrality dependence agree well with the STAR data. One can directly relate and vr^ suppression via 
the (i^ +p)/7r ratio: i?^^ = + {K + p)/'k]aa/[^ + {K + p)/'K]pp. It is clear from the data that 

{K + p)/'K becomes the same for Au + Au and p + p colhsions at pt > 5 GeV/c. To demonstrate the 
sensitivity to the parameterized parton energy loss in the intermediate px region, we also show R\a 
0-5% centrality (dashed line) for jj, = 2.0 GeV and eo = 2.04 GeV/fm without the soft component. 

In the same LO pQCD parton model, one can also calculate di-hadron spectra, 

£^1-^2^3^^^ = ylZy d^bd^rdxadxi,d^kaTd^hTtA{r)tA{\h - r\)gA{kaT,r)gA{hTA^ - A) 
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for two back-to-back hadrons from independent fragmentation of the back-to-back jets. Let us assume 
hadron hi is a triggered hadron with pTi = p%^^- One can define a hadron-triggered FF as the back-to- 
back correlation with respect to the triggered hadron: 



similarly to the direct-photon triggered FF [273, 274] in 7-jet events. Here, zt = Pt/Pt^^ and integration 
over |yi.2| < Ay is implied. In a simple parton model, the two jets should be exactly back-to-back. The 
initial parton transverse momentum distribution in our model will give rise to a Gaussian-like angular 
distribution. In addition, we also take into account transverse momentum smearing within a jet using a 
Gaussian distribution with a width of (kx) = 0.6 GeV/c. Hadrons from the soft component are assumed 
to be uncorrelated. 




Fig. 44: Hadron suppression factors in Au + Au collisions 
as compared to data from STAR[23, 180] and PHENIX [18, 
25]. See text for a detailed explanation. 



Fig. 45: Back- to-back correlations for charged hadrons with 
Pr'*^ >PT> 2GeV/c,py «^ = 4^6GeV/cand < 0.7in 
Au + Au (lower curves) and p + p (upper curves) collisions 
as compared to the STAR[24] data. 



Shown in Fig. 45 are the calculated back-to-back correlations for charged hadrons in Au + Au 
collisions as compared to the STAR data [24]. The same energy loss that is used to calculate single hadron 
suppression and azimuthal anisotropy can also describe well the observed away-side hadron suppression 
and its centrality dependence. In the data, a background B{pt) [1+2^1 {pt) cos(2A(/))] from uncorrelated 
hadrons and azimuthal anisotropy has been subtracted. The value of V2{pt) is measured independently 
while B{pt) is determined by fitting the observed correlation in the region 0.75 < \(t)\ < 2.24 rad [24]. 

With both the single spectra and dihadron spectra, the extracted average energy loss in this model 
calculation for a 10 GeV quark in the expanding medium is {dE/dL)i(i 0.85 ± 0.24 GeV/fm, which 
is equivalent to dE^/dL 13.8 ± 3.9 GeV/fm in a static and uniform medium over a distance = 6.5 
fm. This value is about a factor of 2 larger than a previous estimate [280] because of the variation of 
gluon density along the propagation path and the more precise RHIC data considered . 

Integrating over cf), one obtains a hadron-triggered FF, D^^^^{zt,P%^^) = /J/2 d(j)D^^^^{zT, (t),p^rP^). 

Shown in Fig. 46 are the suppression factors of the hadron-triggered FF's for different values of p^'^ in 
central Au + Au collisions as compared to a STAR data point that is obtained by integrating the observed 
correlation over 7r/2 < |A0| < vr. The dashed lines illustrate the small suppression of back-to-back cor- 
relations due to the initial nuclear kr broadening md + A collisions. The strong QCD scale dependence 
on p^'^ of FF's is mostly canceled in the suppression factor. The approximately universal shape reflects 
the weak pT dependence of the hadron spectra suppression factor in Fig. 44, due to a unique energy de- 
pendence of parton energy loss. Shown in Fig. 47 are the suppresion factors for the direct-photon-tagged 
jet fragmentation function. They are very similar to the direct-triggered fragmentation function, except 
that the photon's energy is more closely related to the original jet energy. 




Fig. 46: The suppression factor for hadron-triggered frag- 
mentation functions in central (0-5%) Au + Au (d-l-Au) col- 
lisions as compared to the STAR data [24]. 



Fig. 47: The modification factor of the photon-tagged inclu- 
sive jet fragmentation function in central Au + Au collisions 
at ^ = 200 GeV for a fixed dEq/dx = 1 GeV/fm. 



3.6 Other Possible Medium-Modifications of High-p^ Hadronic Spectra 

3.61 Recombination Models at the LHC 
R.J. Fries 

Recent results from RHIC show interesting phenomena in hadron production at intermediate trans- 
verse momenta of 2 to 5 GeV/c. This is a region where perturbative QCD starts to be a valid description 
of hadron dynamics, but non-perturbative effects can still be expected to play a crucial role. The key 
observations at RHIC are the anomalous enhancement of baryon production, seen e.g. in a p/vr" ratio of 
about one and the lack of nuclear suppression for baryons between 1.5 and 4 GeV/c [28] and the different 
behavior of elhptic flow for mesons and baryons [27, 249]. 

These observations have lead to the hypothesis that hadron production at intermediate px is domi- 
nated by recombination from a hot and dense parton phase instead of fragmentation of fast partons from 
hard scatterings [136, 149, 167, 268]. In the recombination picture a quark-antiquark pair close in phase 
space can form a meson at hadronization, while three (anti)quarks can find together to be a (anti)baryon. 
The spectrum of mesons from recombination can be written as [135, 137] 

1 
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if the energy E of the meson is large compared to Aqcd and the constituent quark masses. Here <j)M{x) 
is the meson wave function in light cone coordinates, x is the momentum fraction of one of the quarks, 
CAf is a degeneracy factor and the w{(j;p) are classical phase space distributions of the partons before 
hadronization. Transverse momenta relative to the hadron momentum have been integrated out in this 
equation. A similar expression can be found for baryons. 

One can easily show that recombination is more effective than fragmentation for an exponential 
parton spectrum. On the other hand, fragmentation will win over recombination at high Pt if the parton 
spectrum follows a power law. It has been shown that the shape of the wave function plays only little 
role if the input parton spectrum is exponential [137]. It is therefore a good approximation to assume that 
the momentum is equally shared by the valence quarks (1/2 for mesons, 1/3 for baryons). Note that the 
thermal parton phase at hadronization is assumed to have effective degrees of freedom with constituent 
quarks and no dynamical gluons. 

It turns out that all spectra, the nuclear suppression factors and the anisotropic flow coefficient 
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Fig. 48: Spectra of n (left) and p (right) as a function of transverse momentum pr at midrapidity for central Pb+Pb collisions 
at ^/s — 5.5 TeV. The recombination from the thermal parton phase (long dashed line), fragmentation with energy loss from 
LO pQCD (dotted line) and the sum of both (solid line) are shown. For tt" we also give the recombination contribution for 
different values of the the radial flow vt ~ 0.65c (lower short dashed line) and vt = 0.85c (upper short dashed line). 

V2 for hadrons in Au+Au collisions at RHIC foipT > 1-5 GeV/c can be explained by the competition 
between recombination from a thermalized parton phase with temperature T = 175 MeV and radial flow 
velocity vt = 0.55c (for central collisions) and pQCD fragmentation of hard partons including energy 
loss [137]. It is worthwhile to note that it is the strong energy loss of partons in the medium that allows 
recombination to dominate for < 4 GeV/c for mesons and for pT < 6 GeV/c for baryons. 

Numerical estimates for LHC In above calculations for RHIC the parameters for the parton phase 
were determined to match existing data on hadron production. To present estimates for LHC, we fix the 
temperature of the parton phase at hadronization again at 175 MeV as predicted by lattice QCD [176]. 
The average radial flow will be increased at LHC compared to RHIC. We choose vt = 0.75c as the 
radial flow velocity in accordance with [186]. The geometric assumptions about the fireball remain the 
same as those for RHIC [137]. This is certainly a lower bound for LHC. 

The contribution from fragmentation is calculated in leading order (LO) pQCD using the parton 
spectrum given in [250] and KKP fragmentation functions [182]. The partonic energy loss is taken into 
account as in [137]. Its magnitude is fixed to match the mean nuclear suppression factor of about 0.1 for 
a 10 GeV pion at LHC estimated in [263]. 

In Fig. 48 we show the spectra for neutral pions and protons for central Pb+Pb collisions at ^/s = 
5.5 TeV. The recombination part for vr*^ is also given for two other values of the radial flow to estimate the 
theoretical uncertainty that is inherent in our ansatz for the parton phase. Larger emission volumes could 
shift the recombination curve trivially up without changing the slope while leaving the fragmentation 
contribution nearly unchanged. 

The cross over between the fragmentation domain and the recombination domain is at about 6 
GeV for pions (4 GeV at RHIC) and 8 GeV for protons (6 GeV at RHIC) using vt = 0.75c. A larger 
hadronization surface, as likely, will shift these values to even higher pT- In Fig. 49 we show the ratio 
/p from our calculation in comparison with the same quantity calculated for RHIC [137]. We note 
that the surprising baryon enhancement is shifted to even higher transverse momenta at LHC. 

3.62 Transverse Momentum Diffusion and the Broadening of the back-to-back Di-Hadron Correlation 
Function 

J. W. Qiu and I. Vitev 

Multiple parton interactions in relativistic heavy-ion reactions result in transverse momentum dif- 
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Fig. 49: The ratio p/tt" in central Pb+Pb collisions at LHC, \/s — 5.5 TeV, (solid line) and for Au+Au collisions at RHIC, 
= 200 GeV, (dashed hne) [137]. 



fusion and medium induced non-abelian energy loss of the hard probes traversing cold and hot nu- 
clear matter. The corresponding modification of the single inclusive hadron spectra carries information 
about the dynamical properties of the medium created in such reactions and constitutes the basis for 
"jet tomography" [263, 198, 162, 56, 280, 240, 264]. It has been demonstrated that the competition be- 
tween nuclear shadowing, multiple scattering and jet quenching may lead to distinctly different enhance- 
ment/suppression pattern of moderate and high-p^- hadron production inp + A and A + A collisions at 
SPS, RHIC and the LHC [263]. Additional experimental tools that can complement the single inclusive 
measurements, however, are highly desirable. A natural extension of the jet-tomographic technique, first 
quantitatively discussed in [160], is "di-jet tomography". In this case the medium response to the prop- 
agation of hard partons leads to an associated increase of di-jet acoplanarity [53, 87, 90], measured via 
the broadening of the back-to-back di-hadron correlation function [236, 265], as well as to a quenching 
of the away-side Gaussian [265, 165, 281]. These experimental observables are potentially free of the 
uncertainties in Glauber scaling of the baseline p + p result that are present in the comparison of single 
inclusive spectra. 

Particle production from a single hard scattering with momentum exchange much larger than 1/fm 
is localized in space-time. It is multiple parton scattering before or after the hard collision that is sensitive 
to the properties of the nuclear matter. By comparing the high-pr observables in p + p, p + A and 
A + A reactions, we are able to study the strong interaction dynamics of QCD in the vacuum, cold 
nuclear matter and hot dense medium of quarks and gluons, respectively. We here address the elastic 
(no-radiation) scattering of jets {-s/p^ jp^ ~ 0) in nuclear matter [236, 265, 163, 211, 212, 235] that is 
sensitive to the zeroth line integral moment, J dz z^p{z) oc {L)/\ = x, of the matter density. A closed 
form solution can be obtained via the GLV reaction operator approach [163]. Recently, we computed 
the power corrections due to the recoil of the medium [236] and related the momentum distribution of 
partons that have traversed nuclear matter to their initial distribution as follows: 

d^Nf{p+,pT), Ax 



dp^d'^PT 



For any initial jet flux the opacity series in Eq. (129) is most easily resummed in the impact parameter 
space {b~ , b^) conjugate to p^). For the case of a normalized forward monochromatic beam = 




Fig. 50: Left panel: predicted enhancement of (|kT y\) and a Far in minimum bias d + Au and central Au + Au reactions 
at RHIC from pT-diffusion [236]. Preliminary p + p data is from PHENIX [238]. Right panel: the broadening of the far-side 
di-hadron correlation function in central d + Au and Au + Au compared to scaled (xlO) STAR data [24, 16]. In the bottom 
right panel the broadening with and without suppression, approximately given by Raa, are shown. We have used the predicted 
quenching factor [263], confirmed by experimental data [25, 66, 14]. 



(P, Or, P) in the small angle scattering limit we find: 
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The medium induced px broadening and the corresponding longitudinal momentum reduction can be 
evaluated from Eq. (130): 
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In Eqs.(131,132) the factor 2 comes from 2D diffusion, ^ ~ and is the effective gluon density. 
For the 1+lD Bjorken expansion scenario At is the transverse area of the interaction region, tq is the 
initial equilibration time and dN^ /dy is the effective gluon rapidity density. We note that —Ap\\ may 
mimic small elastic energy loss if the full structure of d^N^ /dp'^d'^pr is not observed. 

As an application of the multiple initial and final state elastic scattering formalism elaborated 
in [236, 163] we consider the nuclear induced broadening of the back-to-back jet correlations associ- 
ated with hard QCD ab —>■ cd partonic subprocesses. We will limit the discussion to the Gaussian 2D 
random walk approximation, Eq. (130), to make use of its additive dispersion property. Measurements 
of intra-jet correlations find an approximately Gaussian jet cone shape. If one defines (Ijryl) to be the 



average particle transverse momentum relative to the hard scattered parent parton in the plane normal 
to the collision axis, it can be related to the width a Near of the near-side (A^ < 7r/2) di-hadron corre- 
lation function C(A(/)) = N^^^^''{^^)/N^^^^^ as follows: {\]Ty\) = (|pt|> sin(c77Vear/\/7r). It is the 
away-side {/S.(j) > vr/2) correlation function, however, that measures the di-jet acoplanarity. The total 
vacuum+nuclear induced broadening for the two partons in a plane perpendicular to the collision axis in 
p + A{A + A) reads [236]: 
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A typical range for the cold nuclear matter transport coefficient for gluons (^^/Ag)e//, iS'x^fs = 2 x 
0.1 GeV^/fm - 2 x 0.15 GeV^/fm is extracted from the analysis of low energy p + A data [263, 
264]. This can be tested via the predicted Cronin enhancement in d -|- Au collisions at RHIC s/s = 
200 AGeV [263, 264], which compares well to BRAHMS, PHENIX, PHOBOS and STAR measure- 
ments [265, 16, 62, 29, 67]. For the case of FS scattering in a 1+lD Bjorken expanding quark-gluon 
plasma the final state broadening can be evaluated from Eq. (131). The relation between (Ik^-j^l) = 

\ / {'^j') 1 parton / y O^x) '^ parton — (kji)/2, the ncar-sidc and away-side widths cpfeari (^Far 

and (IptI) 

in the hard fragmentation z = ph/pz ^ 1 limit is approximately given by: 



il^Tyl) = (IptI) cos 
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The left panel in Fig. 50 shows two measures of the predicted increase in di-jet acoplanarity 
for minimum bias d + Au and central Au + Au reactions [236]: (|krj^|) and the away-side width 
apar of the di-hadron correlation function C(A0) = N^^^'^'^ {/\4>) / N^^^'^'^ . C{A(p) is approximated 
here by near-side and far-side Gaussians for a symmetric p!f- Ri p!^ case and the vacuum widths are 
taken from PHENIX [238]. In the right panel of Fig. 50 di-hadron correlations in d + Au are shown 
to be qualitatively similar to the p + p case and in agreement with STAR measurements [24, 16]. In 
Au + Au reactions at RHIC di-jet acoplanarity is noticeably larger, but this effect alone does not lead to 
the reported disappearance of the back-to-back correlations [24]. To first approximation the coefficient 
of the away-side Gaussian (the area under C(A0), Ac/) > 7r/2), is determined by jet energy loss and 

2/3 

given by Raa oc A'^^^^ in the GLV approach [158, 159, 164]. Broadening with and without away-side 
quenching is shown the bottom right panel of Fig. 50. Combined d + Au and Au + Au experimental 
data in Fig. 50 also rule out the existence of monojets at RHIC. For further discussion on di-hadron 
correlations see [236, 165, 281]. 

The broadening of the away-side di-hadron correlation function in p + Pb and Pb + Pb reactions 
at the LHC is shown in Fig. 51. The near side width a Near at (|pr|) = 6 GeV, 8 GeV is extrapolated 
from the PHENIX [238] and STAR [24, 16] measurements. We use two baseline values for the vacuum 
radiation induced di-jet acoplanarity in p reactions to account for its possible growth with -s/s and pr 
relative to the PHENIX measurement, {\]iT y\)vac = 0.8 GeV, 1.2 GeV. In reactions the broadening 
of C(A(/)), A(j) > 7r/2 comes from transverse momentum diffusion in cold nuclear matter. The band 
reflects a range of transport coefficients //^/A = 0.1 GeV^/fm - 0.14 GeV^/fm as in Fig. 50 with ~ 30% 
increase of (jk-r j;|) relative to the vacuum case. In central Pb + Pb reactions, where the hot and dense 
quark-gluon plasma is expected to be formed, the away-side width apar grows by approximately a factor 
of two relative to the p + p case. The final state scattering strength in proportional to the gluon rapidity 
density of the medium and the band represents values in the range dN^ /dy = 2000 — 3500. In the 
right panel of Fig. 5 1 the broadening of the di-hadron correlation function in central Pb + Pb with and 
without the corresponding suppression factor is shown. We note that a direct calculation that does not 
include the hydrodynamic feedback at the LHC energy and number densities will result in suppression 
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Fig. 51: Left panel: di-jet acoplanarity quantified via the enhancement of {|kT y\) in min. bias p + Pb and central Pb + Pb 
reactions at the LHC. Two assumed baseline values for p + p vacuum broadening, {|kT y\)vac = 0.8 GeV, 1.2 GeV. Right 
panel: the broadening of the far-side di-hadron correlation function in central Pb + Pb collisions with and without suppression. 
The quenching factor Raa is taken from [263]. 

factors Raa < {Npart / / [263]. In this case the Raa has been set to {Npart/'^) /Nhin ~ 0.12 for 
Ph + Ph at ^/s = 5.5 TeV. Because of its significantly larger dynamical pT range, LHC may offer the 
best possibility to explore the relation between single inclusive hadron suppression and the broadening 
and disappearance of the back-to-back jet correlations. 

3.63 High-pT Particle Production in Saturation Models 

R. Baier, U.A. Wiedemann 

At RHIC, the production of high px hadrons in central Au-Au collisions shows substantial dif- 
ferences compared to elementary p-p collisions, see section 4.. As discussed in previous sections, the 
observed depletion/suppression may be explained due to induced multiple gluon radiation off the large 
Pt parton ("jet quenching"). Here, we consider an alternative possibility due to initial state gluon radia- 
tion effects, especially saturation effects, put forward first in [178]. 

Instead of a detailed prediction for hadron production based on saturation models [1], we concen- 
trate in the following on the problem of suppression vs. enhancement of gluon production in A-A col- 
lisions with shortly mentioning the relevant comparison in p-A scattering. This way, we do not include 
the fragmentation functions of gluons into hadrons and their possible medium dependence discussed in 
section 3.42. To proceed, we use the kx factorised formalism for calculating gluon production which is 
expected to give a qualitatively reasonable description of this process [150, 178, 78]. 

The basic factorised formula for the gluon yield at central rapidity in a collision of identical nuclei 

is 

dN A7r^asSAA{b)Nc 1 



^d3^" iV^^l ^ J d kt(t>A{y,kt)(t>A{y,PT-h)- (135) 

Here SAA{b) is the overlap area in the transverse plane between the nuclei at fixed impact parameter 



b, and y is the rapidity difference between the central rapidity and the fragmentation region. (pA{y, h) 
is the intrinsic momentum dependent nuclear gluon distribution function, related to the standard gluon 
distribution by 

(,36, 

In the following, we will also use the modified gluon distribution 

hA{kt) = k^VlMh), (137) 

which enters some calculations of the gluon yield instead of (f)A [98, 179, 36]. In general, and especially 
at low momenta, the distributions Ha and (pA are different. However, they coincide for the leading order 
perturbative distribution which at impact parameter b has the shape 

«fc„6).5i(|filfl^^. (138) 

Here, Ppartib) is the density of participants, i.e. for central coUisions /Opart(ft) oc A^/^ and Saa oc A^I'^. 

Since models for the gluon distribution are reviewed in chapter [1], we Umit the present discussion 
to shortly reviewing their main features. 

McLerran-Venugopalan gluon The McLerran-Venugopalan model [215, 216] achieves saturation by 
taking into account the Glauber-Mueller multiple scattering effects. The intrinsic gluon distribution in 
this model was calculated in [188, 172]: 

We will take the saturation momentum to be x-dependent 

'^2(x^^.) = ^^.G(x.l/x')^i!|(^, (140) 

with G{x^ kf = 1/x^) being the nucleon gluon distribution: 

Evolved gluons The MV gluon distribution does not contain any evolution in Bjorken x which is 

necessary to explore the energy dependence of the gluon spectrum. Also, small-z evolution leaves a 
distinct imprint in the /c^ -dependence of the produced gluon. It is argued [169, 224] that in the wide 
region of momenta Qs < \kt\ < ^ (typically, Qq ~ GeV), the evolved distribution behaves as 
("geometric scaling") 



(pAih) oc 



(142) 



with the anomalous dimension 7 = 0.64. However, it is important to know how the distribution behaves 
outside this scaling window since this decides about suppression vs. enhancement of large kt gluon 
production. This is illustrated by two models in [78]: 
A gluon distribution v/iihfast (F) crossover is 

^Fjkt) = ^^"^ ( ^ , (143) 



where = 27ra^A''c ^le^iM^ and ^{h) is chosen to approach 1 rapidly (like a power law) for kt ^ 
QVQo [78]. 

A slow (S) crossover can be modelled by a function with fixed anomalous dimension: 



0.64 



(144) 




Fig. 52: Cronin effect in the pt-dependence of gluon production yields for head-on A-A coUisions for = 2 GeV^. The 
solid curve is for the MV-gluon distribution normalised to the perturbative yield, the dot-dashed curve is for the evolved gluon 
distribution (143), and the dashed line is for the evolved gluon distribution (144). 



The non-linear evolution of the nuclear gluon distributions has recently been calculated numer- 
ically [36], using the Balitsky-Kovchegov (BK) evolution equation [79, 189]. The solutions for h{kt) 
plotted versus the scaled variable p = kt/Qs{x) approach a universal soliton-like shape independent of 
initial conditions if evolved sufficiently far in rapidity. This numerical solution indicates that Eq. (144) 
provides a more realistic parametrisation of the evolved gluon distribution than Eq. (143) with a fast 
crossover. 

Coming to gluon production in A-A collisions, one finds with Eq. (135) for the perturbative gluon 
distribution Eq. (138), 



dy dPp' 



S/=0 



This "reference" spectrum scales for all transverse momenta with A'coii ^ SAA{b) Ppart(^) (~ ^^^^ 
b = 0), as expected perturbatively. 

In contrast, using the saturated gluon distribution in the MV model, the gluon yield Eq. (135) is 
suppressed at small momenta compared to the perturbative one and scales with the number of partici- 
pants, A'^part(fe) = SAA{b) Ppart(^) (~ ^1 for 6 = 0). It approaches the perturbative yield (145) from 
above at large px- 

In Fig. 52, we summarise the results for the normalised ratio of central over peripheral (perturba- 
tive) yields, corresponding to the nuclear modification factor 

dyd-^p d'^b 
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Fig. 53: Ratios Raa and RpA of gluon yields in A-A (LHS) and p-A (RHS) for BK evolution, with MV as initial condition 
with QI = 0.1 GeV^ for p and 2 GeV^ for A. Lines from top to bottom correspond toy = {as Nc/n) In 1/a; = 0, 0.05, 0.1, 
0.2, 0.4, 0.6, 1, 1.4 and 2. 

here quoted at 6 = 0. For the MV model, one sees a small but clear Cronin enhancement for momenta just 
above the saturation scale. The distribution (j)^ displays a clear Cronin effect similar to the MV gluon, 
while c/)\ shows uniform suppression for the central/peripheral ratio for all momenta. This illustrates, 
indeed, that the ratio Raa is very sensitive to the way in which the distribution behaves outside the 
scaling window. 

For the numerically evolved gluon distribution h, the factor Raa is shown in Fig. 53. 

The non-linear BK evolution quickly wipes out any initial Cronin enhancement not only on the 

level of single parton distribution functions but also on the level of particle spectra. Several checks estab- 
lish that this behaviour is generic [36]. For 'realistic' initial conditions this disappearance occurs within 
half a unit of rapidity. We note that in our units the evolution from 130 GeV to 200 GeV corresponds 
to Sy ~ 0.1 for = 0.2, and thus is not sufficient to completely eliminate an initial enhancement at 
central rapidity. For forward rapidity, 6y is greater The evolution to the LHC energy corresponds to 
5y ^ I. Thus the BK evolution suggests the reduction of the Cronin effect in d-Au for forward rapidities 
at RHIC and predicts its disappearance for p-A collisions at LHC. 

This numerical study [36] strongly indicates that crossover from the scaling regime to the pertur- 
bative one is very slow and gradual, and that the Cronin effect which is present in the MV gluon is wiped 
out by the quantum evolution at high energies. Thus, (f>\ in Eq. (144) seems to provide a more realistic 
parametrisation of the evolved gluon distribution than 0^ in Eq. (143). 

On the qualitative level, however, we observe that the gluon distributions which lead to the Cronin 
effect in d-Au collisions also lead to the Cronin enhancement in the Au-Au collisions (see Fig. 53). 
And vice versa, if no Cronin effect appears in Au-Au, none is seen in d-Au collisions. Given the recent 
experimental observation of the Cronin enhancement in d-Au collisions at RHIC (see section 4.), this 
supports the view that significant final state ("quenching") effects are needed in order to account for the 
Au-Au data. 

At LHC, due to higher energies, quantum evolution according to the BK equation will suppress 
gluon production in p-A as well as in A-A collisions (see Fig. 53). 



4, EXPERIMENTAL STATUS AT RHIC 



D. d'Enterria 

We summarize the main results on hard scattering processes in Au+Au, p+p, and d+Au collisions 
at y^sjviv = 200 GeV obtained after 3 years of operation at the BNL Relativistic Heavy-Ion Collider 
(RHIC). The main observations so far at RHIC are the following: 

• The high pT yields of inclusive charged hadrons and vr^ in central Au+Au at y^s^ = 130 [18, 23, 
21] and 200 GeV [25, 14, 31, 66, 62], are suppressed by as much as a factor 4-5 compared to 
p+p and peripheral Au+Au yields scaled by T^b (or Ncou)- 

• At intermediate pr's (px ~ 2. - 4. GeV/c) in central Au+Au, at variance with mesons (7r° [25] and 
K's [15]) no suppression is seen for baryons (p, p [19, 28] and A, A [15]), yielding an "anomalous" 
baryon over meson ratio p/ir ~ 1 much larger than the "perturbative" p/n ^ 0.1 - 0.3 ratio 
observed in p+p collisions [47, 50] and in e'^e~ jet fragmentation [9]. 

• The near-side azimuthal correlations of high pT (leading) particles emitted in central and peripheral 
Au+Au reactions [24, 103] are, on the one hand, clearly reminiscent of jet-Uke parton fragmenta- 
tion as found in p+p collisions. On the other, away-side azimuthal correlations (from back-to-back 
jets) in central Au+Au collisions are found to be significantly suppressed [24]. 

• At low pt the strength of the azimuthal anisotropy parameter V2 is found to be large and consistent 
with hydrodynamical expectations for elliptic flow. Above pT ^ 2 GeV/c where the contribu- 
tion from collective behaviour is negligible, V2 has still a sizeable value with a flat (or slightly 
decreasing) behaviour as a function of px [15, 22, 27]. 

• High Pt production in "cold nuclear matter" as probed in d+Au reactions [30, 16, 67, 62] not only 
is not suppressed enhanced compared to p+p collisions, in a way very much reminiscent of the 
"Cronin enhancement" observed in p+A collisions at lower center-of-mass energies [51]. 

AH these results point to strong medium effects at work in central Au+Au collisions, and have 
triggered extensive theoretical discussions based on perturbative or "classicaF'-field QCD. Most of the 
studies on the highp^ suppression are based on the prediction [153, 154, 68] that a deconfined and dense 
medium would induce multiple gluon radiation off the scattered partons, effectively leading to a depletion 
of the high-pT hadronic fragmentation products ("jet quenching"), though alternative interpretations have 
been also put forward based on initial-state gluon saturation effects ("Color Glass Condensate", CGC) 
[178], or final-state hadronic reinteractions [142]. The different behaviour of baryons and mesons at 
moderately high px^ has been interpreted, among others, in terms of "quark recombination" (or coales- 
cence) effects in a thermaUzed partonic (QGP-like) medium [168, 149, 136], whereas the disappearance 
of the back-to-back azimuthal correlations can be explained in both QGP energy loss and CGC monojet 
scenarios. Finally, the large value of f 2 above 2 GeV/c has been addressed by jet energy loss [157], gluon 
saturation [190], and quark recombination [220] models. 

This summary report presents the pT, ^/snn, centrahty, particle-species, and rapidity dependence 
of the inclusive high pT particle production, plus the characteristics of the produced jets and collective 
elliptic flow signals extracted from the azimuthal correlations at large pT, as measured by the four exper- 
iments at RHIC (BRAHMS, PHENIX, PHOBOS and STAR) in Au+Au, p+p, and d+Au collisions. The 
whole set of experimental data puts strong constraints on the different proposed physical explanations 
for the underlying QCD medium produced in heavy-ion collisions at RHIC and at LHC energies. 



4.1 High-pT Hadron Production in p+p Collisions at ^/s = 200 GeV 

p+p inclusive cross-sections: Proton-proton collisions are the baseline "vacuum" reference to which 
one compares the Au+Au results in order to extract information about the QCD medium properties. At 

i/s = 200 GeV, there currently exist three published measurements of high pT hadron cross-sections in 
p+p{p) collisions: \]A\p+p (|r/| < 2.5, pr <1 GeV/c) [37], PHENIX p-Fp ^ 7r° i\r]\ < 0.35, 



PT < 14 GeV/c) [26], and STAR p + p ^ (|r/| < 0.5, pr < 10 GeV/c) [14]. At ^ = 130 GeV, 
an interpolation between the ISR inclusive charged hadron cross-section and UAl and FERMILAB data, 
has been also used as a reference for Au+Au at this value of ^/s. Globally the spectra can be reasonably 
well parametrized by a power-law form ^ • (1 + pt/po)~"' with the parameters reported in Table 4. We 
note that the fit parameters pQ and n are actually strongly correlated via the mean px of the collision: 
{pt) = 2po/{n - 3). 



system 


V^(GeV) 


„mm / GeV \ 

Pt ^~^) 


A / mb X 


PO(^) 


n 


p + p ^ (inel., interpolation [18]) 


130 


0.4 


330 


1.72 


12.4 


p + p^h^ (NSD, UAl) [37] 


200 


0.25 


286 


1.8 


12.14 


p + p^h^ (NSD, STAR) [14] 


200 


0.4 


286 


1.43 


10.35 


p + p ^ vfO (inel., PHENIX) [26] 


200 


1.0 


386 


1.22 


9.99 



Table 4: Parameters of the fit EcPa/dp^ = A ■ {1 + pt/po) " to the inclusive pr distributions of all existing 
P + p{P) hadron (inelastic or non-singly-diffractive) cross-section measurements at ^/s = 200 GeV. 

In general, all experimental results are consistent within each other, although it is claimed [14] 
STAR p-i-p inclusive charged yield is smaller by a factor of 0.79 it 0.18 compared to UAl p + p results 
(approximately independent of pt), the difference due in large part to differing non-singly-diffractive 
(NSD) cross section measured (35 ± 1 mb [37] in the first and 30.0 ± 3.5 mb [14] in the later). (The 
PHENIX high px tt^ cross-section is inclusive and contains, in principle, all inelastic (including diffrac- 
tive) channels.) Standard next-to-leading-order (NLO) perturbative QCD calculations describe well the 
available high px p+p data at ^/s = 200 GeV (see Fig. 54 for vr^). 




Fig. 54: High pr tt" cross-section in p+p collisions at ^/s = 200 GeV (PHENIX) compared to the results of two different NLO 
pQCD calculations: [26] (left), [81] (right). 



p+p azimuthal correlations: PHENIX [238] has studied the azimuthal correlations at high pT in pH-p 
collisions at ^/s = 200 GeV extracting several parameters characterizing the produced jets: 



• Mean jet fragmentation transverse momentum: (|j_Ly|) = 373 it 16 MeV/c, in agreement with 
previous measurements at ISR [49] and showing no significant trend with increasing ^/s. 

• Average parton transverse momentum (fitted to a constant above 1.5 GeV/c): {\k±y\) = 725 
lb 34 MeV/c. The momentum of the pair p± is related to the individual parton {\k±y\) via 

ilP^Dpair = {\k±y\). The extracted ^J {\pW)pair = 1-82 ± 0.85 GeV/c is in agreement 
with the existing systematics of dimuon, diphoton and dijet data in hadronic collisions [52]. 



4.2 High-pr Hadron Production in Au+Au Collisions 

There is a significant amount of high pr Au+Au experimental spectra (pr > 2 GeV/c) measured by the 4 
experiments at RHIC: inclusive charged hadrons at 130 [18, 23, 21] and 200 GeV [14, 31, 66, 62], neutral 
pions at 130 [18] and 200 GeV [25], protons and antiprotons at 130 [19] and 200 GeV [28], at 200 
GeV [15], and A, A at 200 GeV [15]. Moreover, all these spectra are measured for different centrality 
bins and permit to address the impact parameter dependence of high pr production. 

Details on hadron production mechanisms in AA are usually studied via their scaling behavior 
with respect to p+p collisions. On the one hand, soft processes {pT < 1 GeV/c) are expected to scale 
with Npart [83] (and they actually approximately do [17, 65, 115]). On the other hand, in the framework 
of collinear factorization, hard processes are incoherent and thus expected to scale with Ncoii- result at 
RHIC in the high px sector is the breakdown of this Ncou scaling for central Au+Au collisions. Fig. 55 
shows the comparison of the p+p vr^ spectrum to peripheral (left) and central (right) Au+Au spectra, and 
to pQCD calculations. Whereas peripheral data is consistent with a simple superposition of individual 
NN collisions, central data shows a suppression factor of 4 - 5 with respect to this expectation. 




Fig. 55: Invariant tt*^ yields measured by PHENIX in peripheral (left) and in central (right) Au+Au collisions (stars), compared 
to the Ncoii scaled p+p 7r° yields (circles) and to a NLO pQCD calculation [173] (gray line). The yellow band around the scaled 
p+p points includes in quadrature the absolute normalization errors in the p+p and Au+Au spectra as well as the uncertainties 
in Tab- Updated version of Fig. 1 of [114] with final published data [25, 26]. 

It is customary to quantify the medium effects at high px using the nuclear modification factor 
given by the ratio of the A A to the p+p invariant yields scaled by the nuclear geometry {Tab)'- 

d'NfjdydpT 

Raa[Pt) = -j;^ \ ^ ^2 ttO/^ ^ • ^^^^> 

{Tab} X d'^a^p/dydpT 



Raa{pt) measures the deviation of A A from an incoherent superposition of A^A^ collisions in terms of 
suppression {Raa <1) or enhancement {Raa >!)■ 



High pt suppression: magnitude and pt dependence Fig. 56 shows Raa{pt) for (STAR [14], 
left) and tt'^ (PHENIX [25], right) measured in peripheral (upper points) and central (lower points) 
Au+Au reactions at ^snn = 200 GeV. As seen in Fig. 55, peripheral collisions are consistent with p+p 
collisions plus Ncoii scaling as well as with standard pQCD calculations [263, 281], while central Au+Au 
are clearly suppressed by a factor ~ 4 - 5. (Although peripheral STAR charged hadron data seems to be 
slightly above Raa - 1 and PHENIX vr'^ data seems to be below, within errors both measurements are 
consistent with "collision scaling".) 
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Fig. 56: Nuclear modification factor, Raa{pt), in peripheral and central Au-(-Au reactions for charged hadrons (left) and tt'^ 
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(rigiit) measured at -^sjvjv = 200 GeV by STAR and PHENIX respectively. A comparison to theoretical curves: pQCD-I [263], 
pQCD-II [281], is also shown. 



The high pT suppression in central collisions for both and is smallest at pT = 2 GeV/c and 
increases to an approximately constant suppression factor of 1/Raa ~ 4 - 5 over pt = 5-10 GeV/c. 
Above 5 GeV/c the data are consistent within errors with "participant scaling" given by the dotted line 
at Raa ~ 0.17 in both plots (actually both STAR and PHENIX data are systematically slightly above 
this scaling). The magnitude and px dependence of Raa in the range pT = I - 10 GeV/c (corresponding 
to parton fractional momenta xi 2 = Pt/ Vs{e^^^ + e^^^) 2ptI \fs ~ 0.02 - 0.1 at midrapidity), is 
alone inconsistent with "conventional" nuclear effects like leading-twist shadowing of the nuclear parton 
distribution functions (PDFs) [126, 181]. Different pQCD-based jet quenching calculations [263, 281, 
56, 81, 240, 241] based on medium-induced radiative energy loss, can reproduce the magnitude of the vr'' 
suppression assuming the formation of a hot and dense partonic system characterized by different, but 
related, properties: i) large initial gluon densities dN'^ /dy ^ 1000 [263], ii) large "transport coefficients" 
(7o ~ 3.5 GeV/fm^ [56], iii) high opacities L/A 3 - 4 [81], or iv) effective parton energy losses of the 
order of dE/dx ~ 14 GeV/fm [281]. 

The px dependence of the quenching predicted by all models that include the QCD version of the 
Landau-Pomeranchuck-Migdal (LPM) interference effect (BDMPS [74] and GLV [159] approaches) is 
a slowly (logarithmic) increasing function of pT, a trend not compatible with the data over the entire 
measured pt range. Other approaches, such as constant energy loss per parton scattering, are also not 
supported as discussed in [173]. Analyses which combine LPM jet quenching together with shadowing 
and initial-state pT broadening ("pQCD-II" [263] in Fig. 56) globally reproduce the observed flat pT 
dependence of Raa, as do recent approaches that take into account detailed balance between parton 
emission and absorption ("pQCD-I" [281] in Fig. 56, left). 

At variance with parton energy loss descriptions, a gluon saturation calculation [178] is able to 
predict the magnitude of the observed suppression, although it fails to reproduce exactly the flat pT 



dependence of the quenching [14]. Similarly, semi-quantitative estimates of final-state interactions in a 
dense hadronic medium [142] yield the same amount of quenching as models based on partonic energy 
loss, however it is not yet clear whether the px evolution of the hadronic quenching factor is consistent 
with the data or not [14]. 

The amount of suppression for TfO [25] and [14, 31] is the same above w 4 - 5 GeV/c for all 
centrality classes [31]. However, below ~ 5 GeV/c, vr'^'s are more suppressed than inclusive charged 
hadrons in central collisions (as can be seen by comparing the right and left plots of Fig. 56). This is due 
to the enhanced baryon production contributing to the total charged hadron yield in the intermediate pT 
region {pT ~ 1 - 4 GeV/c) in Au+Au collisions [28, 15] (see section 4.2 below). 

High Pt suppression: ^/s^ dependence Fig. 57 shows Raa{pt) for several vr*^ measurements in 
high-energy AA collisions at different center-of-mass energies [116]. The PHENIX Raa{pt) values for 
central Au-i-Au collisions at 200 GeV (circles) and 130 GeV (triangles) are noticeably below unity in 
contrast to the enhanced production (Raa >1) observed at CERN-ISR (min. bias a + a [50], stars) and 
CERN-SPS energies (central Pb-i-Pb [34], squares) and interpreted in terms of initial-state px broadening 
("Cronin effect" [51]). 
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Fig. 57: Nuclear modification factor, -Raa(pt), for 7r° measured in central ion-ion reactions at CERN-SPS [34] (squares), 
CERN-ISR [50] (stars), and BNL-RHIC (triangles [18], circles [25]) energies. 

Fig. 58 shows -R200/13o(pt), the ratio of Au-(-Au charged hadron yields at = 130 and 200 

GeV in 4 centrality classes, compared to pQCD and gluon saturation model predictions [14]. The in- 
crease in high px yields between the two center-of-mass energies is a factor ~2 at the highest pt's, 
whereas at low px, the increase is much moderate, of the order of 15%. The large increment of the hard 
cross sections is naturally consistent with pQCD expectations due to the increased jet contributions at 
high transverse momenta. In the saturation model [178] the increase at high px is accounted for by the 
enhanced gluon densities at ^/s^ - 200 GeV compared to 130 GeV in the "anomalous dimension" xt 
region of the Au parton distribution function. 

PHENIX [31] has addressed the ^/s dependence of high pT production by testing the validity of 
"xt scaling" in Au-i-Au, i.e. verifying the parton model prediction that hard scattering cross sections can 
be factorized in 2 terms depending on ^/s and xt = Ipj-j ^ respectively: 

^% = ^ ^(-r) ^ e'^^ = 0(xrt . (147) 

In (147), F{xt) embodies all the xt dependence coming from the parton distribution (PDF) and frag- 
mentation (FF) functions (PDFs and EEs, to first order, scale as the ratio of pT at different ^/s.), while 
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Fig. 58: ^^200/130 (pt) vs. pt for [h'^ + h )/2 for four different centrality bins measured by STAR compared to pQCD and 
gluon saturation model predictions [14]. 



the exponent n, related to the underlying parton-parton scattering, is measured to be n 4 - 8 in a wide 
range of p + p,p collisions [31]. Fig. 59 compares the XT-scaled hadron yields in ys/v7V = 130 GeV 
and 200 GeV Au+Au central and peripheral collisions. According to Eq. (147), the ratio of inclusive 
cross-sections at fixed xt should equal (200/f30)". On the one hand, scaling holds (in the kine- 
matical region, xt > 0.03, where pQCD is expected to hold) in Au+Au with the same scaling power 
n = 6.3 ± 0.6 for neutral pions (in central and peripheral collisions) and charged hadrons (in peripheral 
collisions) as measured in p+p [31]. This is consistent with equal (pQCD-like) production dynamics 
in p+p and Au+Au, and disfavours final-state effects described with medium-modified FF's that violate 
XT scaling (e.g. constant parton energy losses independent of the parton px)- Equivalently, models that 
predict strong initial-state effects (e.g. gluon saturation) respect xt scaling as long as their predicted 
modified nuclear PDFs are depleted, independently of ^/s, by the same amount at a given xt (and cen- 
trality). On the other hand. Fig. 59 (right) shows that charged hadrons in central collisions (triangles) 
break xt scaling which is indicative of a non perturbative modification of particle composition spectra 
from that of p+p at intermediate p^'s (see section 4.2 below). 
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Fig. 59: XT scaled spectra for vr" (left) and (h'^ + h )/2 (right) measured in central and peripheral collisions at y^s^Tiv" - 
130 and 200 GeV by PHENIX [31]. Central (Peripheral) xt spectra are represented by triangles (squares), and solid (open) 
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130 GeV scaled by a factor of [130/200]'^-^). 



High Pt suppression: centrality dependence In each centrality bin, the value of the high pT suppres- 
sion can be quantified by the ratio of Au+Au over NcoU-scaled p+p yields integrated above a given (large 



enough) px- The centraUty dependence of the high suppression for vr*^ and charged hadrons, given 
by Raa{pt > 4.5 GeV/c), is shown in Fig. 60 (left) as a function of {Npart) for PHENIX data. The 
transition from the NcoU scaling behaviour (Raa ~ 1) apparent in the most peripheral region, {Npart) ^ 
40, to the strong suppression seen in central reactions {Raa ~ 0.2) is smooth. Whether there is an abrupt 
or gradual departure from NcoU scaling in the peripheral range cannot be ascertained within the present 
experimental uncertainties [116]. The data, however, is inconsistent with NcoU scaling (at a 2a level) for 
the 40-60% centrality corresponding to {Npart) ~ 40 - 80 [21, 116], whose estimated "Bjorken" energy 
density (eBj ~ 1 GeV/fm'^) [116] is in the ball-park of the expected "critical" QCD energy density. A 
similar centrality dependence of the high pT suppression is seen in STAR /i^ data (Fig. 60, right) 



< 
< 
oc 



1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 




• 71° (Pt > 4.5 GeV/c) 

* h± (PT>4-5GeV/c) 



PHENIX 



;TTt| 



50 100 150 200 250 300 350 



f 1.4 



1.2 
1 

0.8 
0.6 
0.4 
0.2 








50 100 150 200 250 300 350 

<Np„rt> 



Fig. 60: Evolution of the high pr vr" and suppression, Raa(pt > 4.5 GeV/c), as a function of centrality given by (Npart) 
(PHENIX, left). Same evolution shown as Raa{pt >6.0 GeV/c) for STAR h"^ data [282] (right). 



Npart (instead of Ncou) scaling at high pT is Np^^.^ = 2) expected in scenarios dominated either 
by gluon saturation [178] or by surface emission of the quenched jets [222]. "Approximate" Npart 
scaling has been claimed by PHOBOS [66]: the ratio of central to a fit to mid-central yields in the range 
Pt ^ 2. - 4. GeV/c stays flat as a function of centrality (Fig. 61, left). However, at higher pT values, 
where the suppression is seen to saturate at its maximum value, the centrality dependence of the ratio 
of A'part-scaled Au+Au over p+p yields for vr*^ and measured by PHENIX (Fig. 61, right) does not 
show a true participant scaling {R^^^ > 1 for all centralities). Nonetheless, the fact that the production 
per participant pair above 4.5 GeV/c is, within errors, approximately constant over a wide range of 
intermediate centralities, is in qualitative agreement with a gluon saturation model prediction [178]. 



High Pt suppression: particle species dependence One of the most intriguing results of the RHIC 
program so far is the different suppression pattern of baryons and mesons at moderately high pt's. 
Fig. 62 (left) compares the Ncoii scaled central to peripheral yield ratios ^ for {p + p)/2 and vr*^: Rep = 
(y^e/d(°-l°^°)/^^°J/°^")/(y^e^d(60-92%)/^60-92%^)_ ^^^^ ^ ^ 4 5 ^^y/^ (anti)protons are not 

suppressed {Rep ~ 1) at variance with the pions which are reduced by a factor of 2 - 3 in this pT range. 
If both vr'^ and p, p originate from the fragmentation of hard-scattered partons that lose energy in the 
medium, the nuclear modification factor Rep should be independent of particle species contrary to the 
experimental result. The same discussion applies for strange mesons and baryons as can be seen from 
the right plot of Fig. 62. Whereas the kaon yields in central collisions are suppressed with respect to 
"NcoU scaling" for all measured pT, the yield of A + A is close to expectations from collision scaling in 
the Pt range 1.8 — 3.5 GeV/c. Interestingly, above pT ~ 5.0 GeV/c, the K^K^, A + A, and charged 
hadron yields are suppressed from binary scaling by a similar factor. 

'since the 60-92% peripheral Au+Au (inclusive and identified) spectra scale with Ncoii when compared to the p+p 
yields [25, 14, 28], -Rep carries basically the same information as Raa. 
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Fig. 61: Left: Ratio of central to semi-central h yields normalized by Npart in 4 different pr bins, as a function of Npart 
measured by PHOBOS [66]. The dashed (solid) line shows the expectation for Ncoii (Npart) scaling. Right: Ratio of Au+Au 
over p+p TT*^ and yields above 4.5 GeV/c normalized by Npart (-RaT"^*) ^ function of centrality given by {Npart) as 
measured by PHENIX [25, 31]. The dashed line indicates the expectation for Npart scaling. 
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Fig. 62: Ratio of central over peripheral Ncoii scaled yields, Rcp{pt), as a function of pr for different species measured in 
Au-l-Au colhsions: (p + p)/2 (dots) and 7r° (circles) by PHENIX [28] (left), and A, A (circles) and K",K^ (triangles) by 
STAR [15] (right). 



Fig. 63 (left) shows the ratios of {p + p)/2 over vr*^ as a function of pT measured by PHENIX 
in central (0-10%, circles), mid-central (20-30%, squares), and peripheral (60-92%, triangles) Au+Au 
collisions [28], together with the corresponding ratios measured in p-i-p collisions at CERN-ISR ener- 
gies [47, 50] (crosses) and in gluon and quark jet fragmentation from e~^e~ collisions [9] (dashed and 
solid lines resp.). Within errors, peripheral Au-i-Au results are compatible with the p-i-p and e+e~ ra- 
tios, but central Au-i-Au collisions have a p/vr ratio ~ 4 - 5 times larger. Such a result is at odds with 
standard perturbative production mechanisms, since in this case the particle ratios p/vr and p/ir should 
be described by a universal fragmentation function independent of the colliding system, which favors 
the production of the lightest particle. Beyond p^- w 4.5 GeV/c, the identification of charged particles 
is not yet possible with the current PHENIX configuration, however the measured /i/vr'' ~ 1.6 ratio 



above ~ 5 GeV/c in central and peripheral Au+Au is consistent with that measured in p+p collisions 
(Fig. 63, right). This result together with STAR Rep result on strange hadrons (Fig. 62, right) supports 
the fact that for large pt values the properties of the baryon production mechanisms approach the (sup- 
pressed) meson scaling, thus limiting the observed baryon enhancement in central Au+Au collisions to 
the intermediate transverse momenta pt ^ 5 GeV/c. 

Several theoretical explanations (see refs. in [19, 28]) have been proposed to justify the dif- 
ferent behaviour of mesons and baryons at intermediate pr's based on: (i) quark recombination (or 
coalescence), (ii) medium-induced difference in the formation time of baryons and mesons, (iii) dif- 
ferent "Cronin enhancement" for protons and pions, or (iv) "baryon junctions". In the recombination 
picture [168, 149, 136] the partons from a thermalized system coalesce and with the addition of quark 
momenta, the soft production of baryons extends to much larger values of pT than that for mesons. In 
this scenario, the effect is limited to < 5 GeV, beyond which fragmentation becomes the dominant 
production mechanism for all species. 
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High Pt suppression: pseudorapidity dependence BRAHIVIS is, so far, the only experiment at RHIC 
that has measured high px inclusive charged hadron spectra off mid-rapidity. Fig. 64 (left) shows the 
nuclear modification factors Raa{pt) for central and semi -peripheral Au+Au measurements at mid- 
pseudorapidity (77 = 0) and at r/ = 2.2 [62]. The high pr suppression is not limited to central rapidities 
but it is clearly apparent at forward ry's too. Fig. 64 (right) shows the ratio of suppressions at the two 
pseudorapidity values, i?^ = Rcp{r] = 2.2) / Rcp{rj = 0). The high px deficit at 7? = 2.2 is similar 
to, or even larger than, at r/ = indicating that the volume causing the suppression extends also in the 
longitudinal direction. 

High Pt azimuthal correlations in Au+Au collisions There are two main sources of azimuthal corre- 
lations at high PT in heavy-ion collisions: 

• The fragmentation of hard-scattered partons results in jets of high pT hadrons correlated in both 
rapidity and azimuthal angle. Such correlations are short range (A?? < 0.7, A(/) < 0.75), involve 
comparatively large transverse momentum particles {pT > 2 GeV/c), and are unrelated (in princi- 
ple) to the orientation of the AA reaction plane. 
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Fig. 64: Le/j": Raa (pt) measured by BRAHMS at = and r/ = 2.2 for 0-10% most central and for semi-peripheral (40-60%) 
Au+Au collisions. Right: Ratio Rri of Rep distributions at ?; = 2.2 and rj = 0. Figs, from [62]. 



• Collective (elliptic) flow: The combination of (i) the geometrical asymmetry in non-central AA re- 
actions ("almond"-hke region of the overlapping nuclei), and (ii) multiple reinteractions between 
the produced particles in the overlap region; results in pressure gradients in the collision ellipsoid 
which transform the original coordinate-space asymmetry into a momentum-space anisotropy. The 
amount of elliptic flow (a true collective effect absent in pn-p collisions) is measured by the sec- 
ond harmonic coefficient, V2 = (cos (20)), of the Fourier expansion of the particles azimuthal 
distribution with respect to the reaction plane. 

Additionally, there are other second-order sources of angular correlations like resonance decays, 
final state (particularly Coulomb) interactions, momentum conservation, or other experimental effects 
like photon conversions, which have to be subtracted out in order to extract the interesting "jet-hke" or 
"flow-like" signals. 



4.3 High-pT Azimuthal Correlations: Jet Signals 

Although, standard jet reconstruction algorithms fail below pr ~ 40 GeV/c when apphed to the soft- 
background dominated environment of heavy-ion colUsions, angular correlations of pairs of high pT 
particles [24, 103] have been very successfully used to study on a statistical basis the properties of the 
produced jets. For each event with "trigger" particle(s) with p^- = 4 - 6 GeV/c and "associated" particle(s) 
with pt = 2-4 GeV/c and |?7| < 0.7, STAR [24] determines the two-particle azimuthal distribution 

D{A^) oc — ^ . (148) 

Fig. 65 shows D{A(p) for peripheral (left) and central (right) Au-i-Au collisions (dots) compared to 
D{A4>) from p-Fp colhsions (histogram), and to a superposed cos{2A4>) flow-like term (blue curve). 
On the one hand, the correlation strength at small relative angles (A(6 ^ 0) in peripheral and central 
Au-i-Au as well as at back-to-back angles (Acf) ~ vr) in peripheral Au-i-Au are very similar to the scaled 
correlations in p-np collisions. The near-side peaks in all three collision systems are characteristic of jet 
fragmentation [24] (a result also observed by PHENIX using neutral trigger particles [103]). On the 
other hand, the away-side peak (A0 ~ tt) in central colUsions is strongly suppressed. 

In order to study the evolution as a function of centrality of the the near- side, L>^"^"(A(/) < 0.75), 
and away-side, D^"^"(A(^ > 2.24), angular correlations in Au-hAu compared to p-i-p, D^^, STAR has 
constructed the quantity 

, Jit ■Hmo^-^- - 5(1 + 2-1 cos(2Aj,))l 
/^^d(A,)i,PP • 




Fig. 65: Azimuthal correlations for peripheral {left) and central (right) Au+Au collisions compared to the pedestal and flow- 
scaled correlations in p+p collisions. Fig. from [171]. 



where B accounts for overall background and V2 the azimuthal correlations due to elliptic flow. Fig. 66 
shows Iaa for the near-side (squares) and away-side (circles) correlations as a function of the number 
of participating nucleons (Npart)- On the one hand, the near-side correlation function is relatively sup- 
pressed compared to the expectation from Eq. (149) in the most peripheral region (a result not completely 
understood so far) and increases slowly with Npart- On the other hand, the back-to-back correlation 
strength above the background from elliptic flow, decreases with increasing Npart and is consistent with 
zero for the most central collisions. The disappearance of back-to-back jet-like correlations is consistent 
with large energy loss effects in a system that is opaque to the propagation of high-momentum partons 
or their fragmentation products. 
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Fig. 66: Ratio of Au+Au over p+p integrated azimuthal correlations, Eq. (149), for small-angle (squares, | A<^j < 0.75 radians) 
and back-to-back (circles, \A(f)\ > 2.24 radians) azimuthal regions versus number of participating nucleons for trigger particle 
intervals 4 < p^"" < 6 GeV/c (solid) and 3 < p^"^ < 4 GeV/c (hollow) [24]. 



High pt azimuthal correlations: Collective elliptic flow At low pT the strength of the elliptic flow 
signal is found to be large and consistent with hydrodynamics expectations. Above p^- ~ 2 GeV/c where 
the contribution from collective behaviour is negligible, V2 is found to be still a sizeable signal which 
saturates and/or slightly decreases as a function of pT [22, 27, 15]. The large value V2{pt > 2 GeV/c) ~ 
0.15 implies unrealistically large parton densities and/or cross-sections according to standard parton 
transport calculations [220]. Various interpretations have been proposed to account for such a large V2 



parameter within different physical scenarios. In jet quenching models [157] the resulting momentum 
anisotropy results from the almond-like density profile of the opaque medium (see, however, [247]). 
Calculations based on gluon saturation [190] yield a ("non-flow") azimuthal asymmetry component from 
the fragmentation of the released gluons from the initial-state saturated wave functions of the colliding 
nuclei. Finally, quark recombination effects [220] can naturally enhance the elliptic flow of the produced 
hadrons compared to that of partons. The measured V2{pt) for mesons and baryons shows a distinct 
pattern (Fig. 67): > V2 at low px, « at pj- « 2 GeV/c, and < V2 at higher pt's; which 
further constraints the proposed theoretical explanations. 




Fig. 67: as a function of transverse momentum for identified particles at RHIC: tt*, and p, p from PHENIX (left), and 
K° and A, A from STAR (right). 

Quark coalescence models [220] naturally lead to weaker baryon flow than meson flow at low pT, 
while the opposite holds above 2 GeV/c. This simple mass ordering expectation recombination models 
is confirmed by the identified particle data from PHENIX and STAR (Fig. 68). The fact that the V2 
parameters scaled by the number of constituent quarks (n = 2 for mesons, n = 3 baryons) versus pT / n, 
globally fall in a single curve, supports the scenario where hadrons at moderate pr's form by coalescence 
of co-moving quarks. 
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Fig. 68: Tfie V2 parameter scaled by the number of constituent quarks (n) versus pr/n for vr*, and p,p (PHENIX [27], 
left) and A'" and A, A (STAR [15], right). 



4.4 High-pT Hadron Production in d+Au Collisions 

Proton- (or deuteron-) nucleus collisions constitute a reference "control" experiment needed to determine 
the influence of cold nuclear matter effects in high pT hadro-production. Since final-state medium effects 
are marginal in p,d-i-Au collisions, they are basic tools to ascertain whether models based on initial- or 



final- state QCD effects can explain the distinct hard scattering behaviour observed in Au+Au collisions 
at RHIC. During the third year of RHIC operation, the 4 experiments collected data from d+Au collisions 
at ^/s^ = 200 GeV. The resulting high pr results at mid-rapidity from PHENIX [30], STAR [16], 
PHOBOS [67], and BRAHMS [62] consistentiy indicate the following: 

• High pt inclusive [30, 16, 67, 62] and vr*^ [30] spectra from d-i-Au minimum bias (MB) col- 
lisions are not suppressed but are enhanced compared to p-i-p collisions {RdAu plots in Fig. 69), 
in a way very much reminiscent of the "Cronin effect" observed in fixed-target p-i-A collisions at 
lower ^/s [51]. As a matter of fact, p-i-Au collisions (from neutron-tagged d-i-Au events [30]) show 
a similar behaviour as minimum bias d-i-Au collisions. 

• Above Pt 2.5 GeV/c the nuclear modification factor of inclusive charged hadrons in MB d-i-Au 
collisions saturates at ^ RdAu ~ 1-4. Above 6 GeV/c, STAR and PHENIX vr*^ results seem to 
indicate that RdAu decreases as a function of px, becoming consistent with 1 at around 8 GeV/c. 

• The "Cronin enhancement" for unidentified hadrons at high pT (RdAu ~ 1-35) is larger than for 
neutral pions (Rf^^ 1.1) [30]. 

• The degree of "enhancement" in d-i-Au compared to p-i-p increases with collision centrality [67, 
16], an opposite trend to Au-i-Au results. 

• The azimuthal correlations in MB and central d-i-Au collisions are very similar to that of p-i-p 
and do not show the significant suppression of the away-side peak observed in central Au-i-Au 
reactions [16]. 
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Fig. 69: Top: Nuclear modification factor RdAu for MB d+Au: and tt" (PHENIX, left), (BRAHMS, right). Bot- 
tom: RdAu for measured by PHOBOS in 4 different d+Au centralities (left), and comparison of two-particle azimuthal 
distributions for central d+Au, p+p, and central Au+Au collisions (STAR, right). 



All these results lead to the conclusion that no "cold" nuclear matter (or initial-state) effects, - 
like a strong saturation of the nuclear parton distribution functions in the relevant (x, Q^) kinematical 

''R™T"'(PT = 2-7 GeV/c) ~ 1.35, Rf^^{pT = 2-6 GeV/c) ~ 1.45, i?git™'(pT = 2-5 GeV/c) ~ 1.3. 



region probed by the current experimental setups-, can explain the high pT behaviour in central Au+Au. 
The data suggest, instead, that final-state interactions are responsible of the high pr suppression and the 
disappearance of back-to-back dijet correlations observed at mid-rapidity in central Au+Au reactions. 

In summary, these data put strong experimental constraints on the properties of the underlying 
QCD medium produced in Au+Au reactions at collider energies. Comparison of the energy spectra 
and angular correlations data to the theoretical calculations globally supports pQCD-based models of 
final-state parton energy loss in a dense medium, although non-perturbative effects (Uke e.g. quark coa- 
lescence) are needed in order to explain the baryon-meson differences in yield and f 2 in the intermediate 
Pt window {pT ~ 2 - 5 GeV/c). Theoretical predictions of a strong saturation of the nuclear wave 
function at high energies are also in agreement with most of the data but do not seem to explain consis- 
tently Au+Au and d+Au RHIC results at midrapidity. Coming ion-ion runs at RHIC and, in the mid-term, 
Pb+Pb collisions at LHC energies will help to further strengthen our current understanding of the physics 
of QCD media at high energy densities. 

5. EXPERIMENTAL CAPABILITIES AT LHC 
5.1 Jet Physics with the ALICE Detector 

A. Morsch 

ALICE is a multipurpose heavy ion experiment with excellent tracking and secondary vertex ca- 
pabilities, electron and muon detection and a high resolution 7-spectrometer [38]. In the barrel part of 
the experiment ALICE will measure the flavor content and phase-space distribution, event by event, for a 
large number of particles whose momenta and masses are of the order of the typical energy scale involved 
(temperature Aqcd ~ 200 MeV). However, tracking and particle identification capabilities at central 
rapidity reach also far into the transverse momentum region in which particle production is expected to be 
dominated by hard processes, the production and fragmentation of high transverse momentum partons, 
i.e. jets. 

Since in its present design complete measurements are restricted to charged particles, the AL- 
ICE detector has only limited capabihties to measure the jet energy. However, energy is not the only 
jet observable. On the contrary, it is likely that the most interesting observables which can reveal the 
presence and kind of interactions of partons with deconfined partonic matter and the associated radiation 
of additional gluons (jet quenching) are mainly related to the structure of the jets, i.e. the phase space 
distribution of particle around the jet axis, fragmentation function and jet shape. Recall also that histori- 
cally the measurement of particle transverse momenta relative to the jet-axis has been used to show first 
evidence for gluons radiated from quark jets produced in e+c~ collisions [96]. Similar distributions can 
be measured by ALICE down to very low particle momenta and for identified particles. This analysis 
will be performed as a function of the energy density by varying the centrality and the size of the AA 
collision system. The study of pA collisions will establish the reference for comparison with cold nu- 
clear matter. Moreover, the observables can be studied as a function of the distance of the jet or leading 
particle direction to the reaction plane. 

The STAR experiment at RHIC has shown [192] that the combination of an electromagnetic 
calorimeter with a TPC tracking system is functionally equivalent to an ideal jet detector in a heavy 

ion collision environment. An electromagnetic calorimeter for ALICE (EMCAL) [111] has been pro- 
posed by a group of collaborating US institutes. This would provide an opportunity to measure jet energy 
and the jet structure in heavy ion collisions. 

5.11 Tracking Particle Identification 

ALICE has been designed to measure single- and multi-particle distributions at the highest anticipated 
charged particle multiplicities for Pb-Pb reactions {dN/drj = 8000). The central tracking devices mea- 
sure particles in the pseudo-rapidity range [rjl < 0.9 with full azimuthal coverage. The inner tracker 



(ITS) [39] provides secondary vertex reconstruction of charm and hyperon decays, particle identification 
and tracking of low-momentum particles (pr > 100 MeV). The main tracking system is a Time Pro- 
jection Chamber (TPC) [40] providing robust tracking capability, good momentum resolution and two 
particle separation, in a high particle density environment. The Transition Radiation Detector (TRD) 
[41] for electron identification adds additional tracking information at high radii improving the momen- 
tum resolution of high-px particles. The momentum resolution Ap/p for particle momenta < 2 GeV is 
below 1%. It rises almost linearly to up to 16% for p = 100 GeV. Tracking efficiency reaches almost 
100% for momenta above 1 GeV. A geometrical inefficiency for high-px particles of 10% results from 
dead zones between the TPC read-out chambers. 

A Time of Flight System (TOF) [42] based on multi-gap resistive plate chambers (RPCs) extends 
particle identification in the barrel region to 2 GeV and 3.5 GeV for K/tt and K/p separation, respec- 
tively. Recent studies with an improved TPC tracking algorithm have shown that particle identification on 
a statistical basis via dE/dx can be extended into the relativistic rise up to « 50 GeV. The HMPID (High 
Momentum Particle IDentificationP) system [43] using a single-arm Ring Imaging Cerenkov Counter is 
devoted to the discrimination of hadrons in the hard part of the momentum spectrum. It will enhance the 
PID capability of ALICE by allowing to identify particles beyond the momentum interval covered by the 
energy loss measurements and the TOF. The useful range for the identification of vr /K and K/p of a track 
by track basis is extended up to 3 GeV and 5 GeV, respectively. The PHOS (PHOton Spectrometer) [44] 
is a single-arm high-resolution electromagnetic calorimeter of high granularity. It is optimised for mea- 
suring photons, vr^'s and r] mesons in broad energy ranges. The TRD provides excellent e/7r-separation 
in the 1 — 100 GeV momentum range. 

5.12 Event Geometry 

The combined measurement of energy in the zero-degree calorimeters (ZDC) [45] and the forward elec- 
tromagnetic calorimeter (ZEM) allows ALICE to determine the collision centraUty in five impact param- 
eter bins between and 16 fm. Non-zero impact parameters lead to pressure gradients in the reaction 
volume producing directed flow that can be observed by ALICE as an anisotropy in the azimuthal distri- 
bution of final state particles. The measurement of directed flow allows an unambiguous determination 
of the event plane. Thus jet observables can be studied as a function of the centrality and the distance of 
the jet axis to the reaction plane. 

5.13 Inclusive Transverse Momentum Spectra 

As described in section 3.43 transverse momentum spectra are a mean to study final state radiative 
energy loss via the nuclear modification factor Raa{pt)- Using a sample of 10^ Pb-Pb minimum bias 
events and in the absence of quenching ALICE will be able to measure the charged particle transverse 
momentum spectrum up to ^ 100 GeV and down to 100 McV. 0(10) events are in the highest 
momentum bins where jet quenching is expected to reduce the charged particle yield by a factor 2 — 3. 
ALICE is planning to use a High Level Trigger (HLT) to increase the statistics in the high momentum 
region. Moreover, neutral pions will be identified on an event-by-event basis in the high-momentum 
range from 30 - 100 GeV. 

5. 14 Rates and the Need for Triggering 

ALICE wants to study the whole spectrum of jet production ranging from mini-jets (£'t > 2 GeV) to 
high-£'T jets of several hundred GeV. Concerning the experimental capabilities one has to distinguish 
four energy regions, which are here discussed for central (b < 5fm) Pb-Pb collisions at ^snn = 
5.5 TeV/nucleon: 

• In the region E'x < 20 GeV several jets of these energies overlap in one event within the ALICE 
acceptance. This means jet identification in the ttaditional sense is not possible. However, their 




Fig. 70: Distance R of particles witli pt>3 GeV/c to leading particles with pt>5 GeV/c (left) and f<T> 10 GeV/c (right). The 
dashed line describing the uncorrelated background is obtained using randomized leading particle directions. 



presence reveals in inclusive studies of particle correlations. 

• For Et < 100 GeV the jet rate of > 1 Hz is high enough so that even with a read-out rate Umited 
by the TPC to 40 Hz an event sample of > 10^ jets can be collected in one effective month of 
running (10^ s). 

• For Et > 100 GeV triggering will be necessary to collect enough data. 

• Assuming that for a fragmentation function analysis one needs about 10^ — 10^ events the statistics 
limit is reached at about 250 GeV. 

In the absence of calorimetry, triggering on jets is only possible using a High Level Trigger (HLT). 
Presently ALICE us studying the possibility to trigger on event topologies where two or more high-pT 
tracks are found in a small area of the rj — (f) plane. The search has to be performed on track candidates 
which are themselves the result of a HLT fast clustering and tracking procedure. 

5. 15 Inclusive Jets at Low Et 

In proton-antiproton collisions at 1.8 TeV evidence for low energy charged particle clusters has been 

seen [33]. These charged particle jets become apparent event by event at a charged jet energy of about 
2 GeV with, on the average, two charged particles with > 0.5 GeV and grow to, on the average, 
about 10 charged particles at a charged jet energy of 50 GeV. 

In central Pb-Pb collisions at 5.5 TeV with an anticipated charged particle multipUcity of up to 
8000 charged particles per unit of rapidity the situation will be completely different. It is expected, 
that in one central collision about 100 jets with > 5 GeV, are produced within the acceptance of 
the ALICE central tracking system. The jet multiplicity decreases to an average of one, for pT > 
20 GeV. The individual structures of these low-px jets dissolve in the overall event structure and are 
not distinguishable event by event. Nevertheless, their properties can be studied using inclusive particle 
correlation distributions as will be shown in the following. 

This effort has two main objectives. First the study of the underlying event properties is important 
for the understanding of the limits of the energy resolution for the reconstruction of high-^^T jets in HI 
colhsions, i.e. the underlying event fluctuations. Second, it is expected that in-medium modifications 
of the jet structure will be stronger for low jet energies. ALICE wants to study changes of the particle 
momenta parallel to the jet axis (jet quenching) and in the transverse direction (transverse heating). 

In order to study inclusive particle correlations, three principle methods can be used: 

• Event by event and region by region fluctuations of number of particles or energy; 




Fig. 71: Distribution of z = pT /p!^^'*'"^ for particles within 7?<0.1 of the leading particle direction (solid line). The contribu- 
tion of uncorrelated particles obtained from a guard band l<ii<2 has been subtracted and is shown as a dashed line. 




Fig. 72: Charged particle transverse momentum distribution for particles within R<0.1 of the leading particle direction (solid 
line). The contribution of imcorrelated particles obtained from a guard band l<i?<2 has been subtracted and is shown as a 
dashed line. 

• Correlations with leading particles; 

• Spatial spectrum analysis (autocorrelation). 

These three methods will be outUned in the following and some feasibiUty studies will be pre- 
sented. 

Underlying event fluctuations Here, we will consider fluctuations of the energy contained in a cone 
of size Rq. These fluctuations Umit the energy resolution of high-energy jets obtained with the cone 
algorithm. For uncorrelated particle production, the relation between the rms energy variation AE and 
the number of particles, N, the mean transverse momentum, < >, and the rms of the pT-distribution, 

ApT, is AE = y/N \J <pt +Ap\. As an example, this value increases by a factor \/2, if the 
multiplicity N results from N/2 clusters of multiplicity 2. In fact, in central Pb-Pb collisions simulated 
with HIJING [270, 155], we observe for Rq = 0.3 fluctuations that are by factors of 1.3 — 1.5 higher 
than expected for uncorrelated particles, the exact value depending on the pT-cut. 



Correlations with leading particles In order to study correlations with leading particles we apply an 
algorithm similar to the one used for the CDF charged jet analysis. All particles with a pt > Pt^*^ 
leading particle candidates and are ordered according to their p^. We start with the highest pT candidate 
and record the distances R in the t] — 0-plane between all charged particles and the leading particle. If 
another candidate is found within a distance R < i?sep it is eliminated from the list of leading particles. 
The procedure continues with the next leading particle until none is left. 

This algorithm is a natural extension of the cone algorithm used for jet reconstruction to inclusive 

studies in the low-px jet region for heavy ions collisions. To see possible angular correlations we plot 
the particle density (27rii)^^dA^/di?. Uncorrelated particle production should result in a flat distribution 
close to the leading particle direction R = 0. 

HUING simulations of central Pb-Pb collisions at the LHC centre of mass energy have been 
performed to study these leading particle correlations. Fig. 70 shows such distributions for p^^'^ = 5 
and lOGeV with a cut on all other particles of > 3GcV and i?scp = 0.5. The dashed lines are 
the corresponding distributions for randomized leading particle directions. Clear correlation signals are 
visible for R < 0.3. As expected, the significance of the signal increases with the transverse momentum 
cut. 



Fragmentation function The distributions of the correlated particle transverse momentum normalized 
by the leading particle transverse momentum z = p^ /pl^;^*^™^ is related to the jet fragmentation function. 
In Fig. 71 we show this distribution for particles with R < 0.1. The background distribution obtained 
from particles in a guard band 1 < i? < 2 has been subtracted and is shown as a dashed Une. The signal 
dominates at high z values. The corresponding "raw" pt distributions are shown in Fig. 72. 

pleading ^j^jy ^ ^g^y p^^^ estimator of the jet energy and the dN/dz distribution can only be in- 
terpreted as a smeared out "pseudo"-fragmentation function. This explains the rather modest decrease as 
z approaches 1 . Nevertheless, since in-medium modifications of the fragmentation function are expected 
to be strong for low energy jets, we expect that the effect can be observed in the measured distribution 
by comparing pPb to Pb-Pb data. 

Momentum transverse to the jet axis The distributions of the correlated particle momentum perpen- 
dicular to the jet axis Jt = P ■ sin(0 (particle, leading particle)) is shown in Fig. 73. The background 
shown as dashed Unes is obtained using randomized leading particle direction and has been subtracted. 
Signal and background have a similar distribution. However, this has to be interpreted in the light of the 
results of the previous paragraph, where we saw that the signal has a harder px-spectrum. Ordering of 
Pt in the jet fragmentation leads to limited jt with mean value of about 250 MeV. Higher px particles 
are, on the average, closer to the jet axis. It is expected that in-medium effects can significantly broaden 
the Jt distribution of the particles inside the jet. For the expected S/B > 1 such changes can be easily 
measured. 



Forward-backward correlations The backward jets observed through forward-backward (;^-correlations 
are more difficult. This is due to the Umited r/-acceptance of the experiment which reduces the number 
of two-jet events within the acceptance, especially for low-px jets. Higher statistics or harder cuts are 
needed than for the /^-correlation to obtain significant signals. In Fig. 74 we show the Acf) distribution 
for particles with px > 5 GeV with respect to the leading particles direction with pJ^;^*^"^^ > 10 GeV. 
Already for 10^ HUING events a significant backward peak is observed indicating that such an analysis 
is possible at the LHC. 

Spatial spectrum analysis Jets reveal their presence through repeating length scales. It is well known 
from signal processing that hidden periodicities can be detected using Fourier analysis techniques e.g. 




Fig. 73: Distributions of tlie correlated particle momentum perpendicular to the jet axis Jt = p ■ 
sin(6(particle, leading particle)) obtained from particles within i?<0.1 of the leading particle direction (solid line). The 
contribution of uncorrelated particles obtained by randomizing the leading particle direction has been subtracted and is shown 
as a dashed line. 



the spectral power density (SPD) or its back-transformation the autocorrelation function (ACF). For 
jet analysis the discretized energy density distribution in the ry — plane is used as an input. The 
transformations are performed using a two-dimensional discrete fast Fourier transformation algorithm 
(DFFT). In heavy ion collisions the length scales not only repeat inside the same event but also from 
event to event. In order to profit from this fact, the spectral power densities are averaged over many events 
before the back-transformation is performed to obtain the event averaged auto-correlation function. 

5.16 Reconstructed Jets 

Charged Jets Leading particle analysis provides only a very poor energy measurement. It is known 

from pp data that leading charged particles carry on the average 24% of the charged jet energy. The distri- 
bution of this energy fraction has a rms of 50% hmiting the jet energy measurement to approximately 
the same value. Nevertheless, there exists some limited selection power for jet energy classes, since the 
jet energy can not be lower than the leading particle energy. Due to the steeply falling jet production 
spectrum a leading particle cut of for example > 5 GeV will select low-energy jets mainly in the 
range 5 < Et < 25 GeV. In the low-energy region underlying event fluctuation are of the order of the 
reconstructible jet energy, i.e. the energy of reconstructed charged particles from the jet in a fixed cone 
after pT-cuts. Hence, true jet reconstruction is only possible in the high energy region, > 50 GeV 
being our current estimate. 

Since the energy resolution will be limited due to the incomplete reconstruction, a simple recon- 
struction algorithm is sufficient. The cone-algorithm similar to the one developed by the UAl collabora- 
tion [61] has been used to evaluate the jet energy resolution in Pb-Pb collisions using charged particles 
only. In this algorithm, charged particles above a threshold > p^^'^ are considered as seeds. A jet is 
defined by all charged particles within a cone R < iimax with respect to the jet axis, which for the first 
iteration is given by the leading particle direction. In the next iteration a new jet direction is obtained 
from the sum of the momentum vectors of these jet particles. The procedure is repeated until conver- 
gence is reached, i.e. the difference between new and old direction falls below a minimum value. The 
energy from the underlying event is determined from the particle outside the jet cone and is subtracted 
from the cone energy. In central Pb-Pb colhsions, for a jet E^ of 50 GeV the energy resolution is similar 
to that obtained from the leading particle analysis (AEt/E^ « 50%). For jet energies above 100 GeV 
the resolution improves to 40 %. The resolution is mainly limited by the fluctuation in the part of the jet 
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Fig. 74: /^(f) distribution for particles with pt>5 GeV/c with respect to the leading particle direction with prp'^'iding ^ jq Qgy/^ 
(solid line). The dashed lines represent the uncorrelated background and has been obtained by randomizing the leading particle 
direction. The dotted line indicates the uncorrelated background level below the signal and has been obtained from fit to the 
region between the peaks. 

energy cairied by charged particles (30%) and the small cone-size (i? < 0.3) that is needed to reduce the 
energy from the underlying event. The latter reduces the reconstructed jet energy and thus increases the 
ratio /S.E/E. 

Energy resolution with EMCAL The proposed EMCAL adds additional energy information from 
neutral particles. A linear combination of the energies of charged particles and calorimeter cells that 

minimizes the energy dispersion has been found from MC studies and is used to reconstruct the jet 
energy. In this case jets with E-y > 50 GeV can be reconstructed. The expected resolution is 34% for 
= 50 GeV and 29% for Et: = 100 GeV (preliminary results). 

Fragmentation function with reconstructed jets The fragmentation function represents the average 
distribution of jet energy among the particles produced in the fragmentation of the initial parton. The ef- 
fect of additional gluon radiation when the particle is propagating through a deconfined partonic medium 
will change the fragmentation function in two ways: the relatively low-energetic radiated gluons will 
increase the number of particle carrying a small fraction z of the jet energy. On the contrary, the high-z 
part will be depleted. It is important to observe both effects experimentally. The challenge is to obtain a 
reasonable S/B in the low-z region and enough statistics in the high-z region. 

As an example we show in Fig. 75 the fragmentation function for jets that have a reconstructed jet 
energy in the range 95 < £'t < 105 GeV. The solid line is without jet quenching. The dotted line is for 
a quenching scenario which has been simulated using the superposition of two jets with £'t = 80 GeV 
and Et: = 20 GeV, respectively. The dashed line is the underlying event. 

Fragmentation function with identified leading particles Reconstructed D-mesons and high-px elec- 
trons from semi-leptonic heavy-quark decays with large impact parameter are leading particles per se. 
They tag the production and fragmentation of charm and beauty quarks, respectively. In the past, exper- 
iments like UAl without secondary vertex detection capability for charm and beauty hadrons have used 
non-isolation cuts as an additional criterion to select these heavy particles, where non-isolation is defined 
as the energy in a cone around the particle. 

In ALICE, the secondary vertex reconstruction capabiUty using the ITS allows to select these par- 



xO.5 



Fig. 75: Fragmentation function for jets that have a reconstructed jet energy in the range 95 < Et < 105 GeV. The solid line 
is without jet quenching. The dashed Une is for a quenching scenario which has been simulated using the superposition of two 
jets with Et = 80 GeV and Et = 20 GeV, respectively. The dotted line is the ^-distribution for charged particles from the 
underlying event. 

tides using impact parameter cuts. mesons decaying into ttK can be reconstructed in the transverse 
momentum range < px < 15 GcV. Electrons from b-hadron decays cover the transverse momentum 
range pt > 2 GeV. Additional particles close to the leading particle come from the fragmentation and 
decay of the heavy quark. The jt and pt distribution of these particles can be studied in pp, pPb and 
Pb-Pb collisions. Again, differences between these systems are the result of interactions between partons 
and the medium. They are expected to be weaker for heavy quarks than for light quarks, which makes 
this measurement especially interesting. 

Moreover, for identified charmed mesons the inclusive pT-spectra can be studied in the same 
way as for example charged light hadrons. Comparing Raa{pt) of light and heavy mesons will show 
whether energy loss is similar for both particle types. Theory expects a smaller energy loss for heavy 
quarks (dead-cone effect). 

Photon-tagged jets Photons produced at the earliest stage of the collisions, preserve almost all their 
energy after traversing through the dense medium. Hence, the attenuation of the jet energy can be 
measured via comparison of the prompt photon and jet kinematics. Prompt photons accompanied by a 
jet produced in the opposite direction at high pT can be studied as a probe of the dense medium formed in 
heavy-ion collisions. ALICE will be able to detect and identify prompt photons using the PHOS detector, 
while hadrons from jets will be detected by the TPC and, optionally, by the EMCAL. 

5.2 Jet Physics with the CMS Detector 

O.L. Kodolova, LP. Lokhtin, S.V. Petrushanko, C. Roland, L.I. Sarycheva, S.V. Shmatov, I.N. Vardanian 
5.21 CMS Detector 

The Compact Muon Solenoid (CMS) is a general purpose detector designed primarily to search for the 
Higgs boson in proton-proton colhsions at LHC [104]. The detector is optimized for accurate mea- 
surements of the characteristics of high-energy leptons and photons as well as hadronic jets in a large 
acceptance, providing unique capabilities for "hard probes" in both pp and AA colhsions. In particular, 
jet quenching studies of hard jets, high-mass dimuons and high-py hadrons are primary goals of the CMS 
Heavy Ion Programme [82]. 
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Fig. 76: The longitudinal view of CMS detector. 
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A detailed description of the detector elements can be found in the corresponding Technical Design 
Reports [105, 106, 107, 108]. The longitudinal view of the CMS detector is presented in Fig. 76. The 
central element of CMS is the magnet, a 13 m long solenoid with an internal radius Ri 3 m, which will 
provide a strong 4 T uniform magnetic field. The An detector consists of a 6 m long and 1.3 m radius 
central tracker, electromagnetic (ECAL) and hadronic (HCAL) calorimeters inside the magnet and muon 
stations outside. The tracker and muon chambers cover the pseudorapidity region |ry| < 2.4, while the 
ECAL and HCAL calorimeters reach \r)\ = 3. A pair of quartz-fiber very forward (HE) calorimeters, 
located ±11 m from the interaction point, cover the region 3 < |r/| < 5 and complement the energy 
measurement. The tracker is composed of pixel layers and sihcon strip counters. The CMS muon stations 
consist of drift tubes in the barrel region (MB), cathode strip chambers in the end-cap regions (ME), 
and resistive plate chambers in both barrel and endcap dedicated to triggering. The electromagnetic 
calorimeter is made of almost 76000 scintillating PbW04 crystals. The hadronic calorimeter consists of 
scintillator sandwiched between brass absorber plates. The main characteristics of the calorimeters, such 
as energy resolution and granularity are presented in Table 5. 

5.22 Dijet, "y+jet and Z+jet Production at CMS 

The following signals of jet quenching by medium-induced parton energy loss have been identified as 
being measurable in heavy ion collisions with CMS [82]. 

Jet pairs are produced in the initial scattering processes in pp and AA colhsions through reactions 



Table 5: Energy resolution, a/E, and granularity of the CMS calorimeters in the barrel (HB, EB), endcap (HE, 
EE) and very forward (HF) regions. The energy resolution is shown for the total energy of electrons and photons 
(ECAL) and transverse energy of hadronic jets (HCAL, HF). 



Rapidity 
coverage 


<| 7/ |< 1.5 
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3.0 <| r] \< 5.0 


Subdetector 


HCAL (HB) 


ECAL (EB) 


HCAL (HE) 


ECAL (EE) 


HF 


a 
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1.16 
0.05 


0.027 
0.0055 


0.91 
0.05 


0.057 
0.0055 


0.77 
0.05 


granularity 
A?7 X A<^ 


0.087 X 0.087 


0.0174 X 0.0174 


0.087 X 0.087 


0.0174 X 0.0174 
to 0.05 X 0.05 


0.175 X 0.175 



such as 

99^ 99 , Q9^Q9 , ^ 99 , 99^09, 

where the gg gg process is dominant. High px jet pairs produced in nucleus-nucleus collisions can be 
suppressed relative to their production in independent nucleon-nucleon interactions [153]. The dijet rates 
depend on impact parameter [204] and may exhibit azimuthal anisotropy in non-central collisions [205]. 

Single jets may be produced opposite a gauge boson in 7-l-jet [273, 274] and Z+jet [177] final 
states, dominantly through processes such as 

Q9 -^qi , qg qZ ■ 

In heavy ion collisions, the relative px between the jet and the boson becomes imbalanced due to inter- 
actions of the jet within the medium. The Z is detected in the dimuon channel, Z 

In this section, we will consider mainly production and measurement of hard jets, the two points 
described above. Although not discussed in detail here, CMS can also measure parton energy loss in two 
other channels, leading particles and heavy quarks. Leading particles in a jet may have their momentum 
suppressed due to medium modifications of the jet fragmentation functions [271]. The capabihty of the 
CMS tracker to measure the momenta of charged particles in heavy ion collisions is discussed later. 
Heavy quark energy loss, particularly b quark loss, can be measured in two channels. Semileptonic B 
and D decays contribute to the high mass dimuon spectra and hadronic B decays to J/ tp are a substantial 
contribution to secondary charmonium production [202, 206]. 

By considering "jet energy loss" here, we concentrate on the energy that falls outside the jet 
cone and is truly lost to the jet, see Refs. [203, 71]. In fact, since coherent radiation induces a strong 
dependence of the radiative energy loss of a jet on the angular cone size, it will soften particle energy 
distributions inside the jet, increase the multiplicity of secondary particles, and, to a lesser degree, affect 
the total jet energy. On the other hand, coUisional energy loss turns out to be practically independent 
of jet cone size and causes the loss of total jet energy. Moreover, the total energy loss of a jet will be 
sensitive to the experimental capabiUties for low-p^ particles, products of soft gluon fragmentation. In 
CMS, most of these low-pr hadrons may be cleared out of the central calorimeters by the strong magnetic 
field. 

Table 6 presents the event rates for various channels, including hard jets, in a one month Pb— Pb 
run (assuming two weeks of data taking), T = 1.2 x 10^ s, with luminosity L = 5 x 10^^ cm~^s~^ so 
that 

iVev = Ta'Xj^L . 



The production cross sections in minimum bias nucleus-nucleus collisions were obtained from those in 
pp interactions at the same energy, i/s = 5.5 TeV, using the simple parameterization = cr^A^. 
The pp cross sections were evaluated using PYTHIA 6.1 [248] with the CTEQ5L parton distributions. 
The jet production cross sections in the CMS acceptance will be large enough to carefully study the dijet 
rate as a function of impact parameter as well as the azimuthal angle and rapidity distributions of jet 
pairs. The estimated statistics for 7-l-jet production are satisfactory for studying the ^^T-imbalance of 
the process but the large background from jet +jet(^ vr'') is still under investigation. The corresponding 
statistics for Z{—i- /x"'"/x~)+jet are rather low, but the background is less than 10% in this case. 



Tabic 6: Expected rales for jet production channels in a one month Pb— Pb run. 



Channel 


Barrel 


Barrel+Endcap 


jet+jet, E}^^ > 100 GeV 


2.1x10^ 


4.3 xlO^ 


7+jet, E^J.^*'^ > 100 GeV 


1.6x10^ 


3.0x10^ 


fi+fi-)+}et, E}^\p^ > 100 GeV 
/x+/x-)+jet, > 50 GeV 


30 
180 


45 
300 



Of course, there are some theoretical uncertainties in the absolute jet rates in pp collisions due 
to the choice of parton distribution functions, the importance of next-to-leading order of ag corrections, 
etc. Thus jet measurements in pp or dd collisions at the same or similar energies per nucleon as in the 
heavy ion runs are strongly desirable to determine the basehne rate precisely. One complementary way 
to reduce uncertainties in the analysis of jet quenching is the introduction of a reference process unaf- 
fected by medium-induced energy loss and with a production cross section proportional to the number 
of nucleon-nucleon collisions, such as Z{-^ production [82, 177]. 



5.23 Jet Reconstruction 

The main difficulty of jet recognition in heavy ion collisions arises from the "false" jet background - 
transverse energy fluctuations coming from the high multipUcity of "thermal" secondary particles in the 
event [193]. Predictions give between 3000 to 8000 charged particles per unit rapidity in central Pb— Pb 
collisions at the LHC. In these circumstances, reconstruction of "true" QCD jets resulting from hard 
parton-parton scatterings is very important. The definition of an object like a jet is quite non-trivial even 
in pp collisions. In particular, the energy and spatial resolution of hard jets are sensitive to the parameters 
of the jet finding algorithm, see Ref. [130]. 

In the CMS heavy ion physics programme, the modified sliding window-type jet finding algorithm 
has been developed to search for "jet-like" clusters above the average energy and to subtract the back- 
ground from the underlying event [82, 184], see Fig. 77 for an illustration. The algorithm is described 
step-by-step here. 

• As a function of pseudorapidity rj, one calculates for each event the average transverse energy, E^^^{r]), 

2 



and the dispersion, "(ry) = y (Efp^^ (r]))^ — \^E^^^{r])j , in all the calorimeter cells. The superscript 
"cell" means averaging over the calorimeter cells in the given event. 

• All possible rectangular windows, with overlaps, in the calorimeter map of rj — ip space are constructed. 
Each window consists of an integer number of calorimeter cells. The numbers of cells per window in rj, 
A^^™*^, and </?, A^™™*^, are calculated separately, 

where A^*°* and are the total number of calorimeter cells in 77 and ip, r/max is the maximum pseudo- 
rapidity and R is the jet cone radius, an external parameter of the algorithm. 




Fig. 77: Sliding window-type jet finding algorithm. 



• The window energy is calculated as a sum of the call energies exceeding background, which is la, 
D^^^{ri), above the averaged energy, E'^^". If the transverse energy of the cell is negative after back- 
ground subtraction, it is set to zero. 

• The search for jets and the evaluation of their energies is started from the window with the maximum 
transverse energy. 

• Non-overlapping windows with energy E^^^'^ > E™^ are considered to be jet candidates. 

• The center of gravity of the window is considered as a center of the jet. 

• For correction of the jet axis a cell with maximum transverse energy in cone is found and considered 
as a new geometrical center of this jet. Cells within radius R around the new geometrical center are 
collected and center of gravity of jet is recalculated. 

• Cells in a cone with radius R around jet center are collected. 

• The values of E^^^\r]) and D^^\r]) are recalculated using cells which were not included in the jets. 

• The jet energy is then the difference between energies in collected cells, E^^^, and the background 
energy per cell, 

4^* = j2{Ef' - \E^) + Df'm. 

Tie 



Table 7: Purity, noise (contamination levels, false jets / generated jets) and jet transverse energy resolution in 
central Pb-Pb collisions with dN^/dy{y = 0) = 3000 and 8000. 



Et min 


Purity 


Noise 


a{ET)/ET{%) 


(GeV) 


3000 


8000 


3000 


8000 


3000 


8000 


75 


0.96 ± 0.03 


0.88 ± 0.03 


0.021 ±0.006 


0.083 ± 0.009 


19.0 


17.8 


100 


0.99 ± 0.03 


0.97 ± 0.03 


0.002 ± 0.001 


0.011 ±0.003 


16.4 


18.4 


125 


1.00 ± 0.03 


0.99 ± 0.03 


0.000 ± 0.000 


0.004 ± 0.002 


15.1 


16.8 


200 


1.00 ± 0.03 


0.99 ± 0.03 


0.000 ± 0.000 


0.001 ± 0.001 


10.7 


12.7 



Jet reconstruction was studied in the barrel calorimeters, |r/| < 1.5, with the GEANT-based pack- 
age CMSIM_123 (CMS simulation package, version 123) adapted to heavy ion collisions. The initial 
jet distributions in the nucleon-nucleon sub-collisions at ^/s = 5.5 TeV were generated using PYTHIA 
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Fig. 78: Correlation between reconstructed and gener- 
ated jet transverse energies in Pb— Pb and pp events. 



Fig. 79: Jet energy resolution in Pb— Pb and pp events. 
(See text for explanation). 



6.1 [248], as described before. This dijet event is then superimposed on the Pb-Pb event, obtained using 
a hydrodynamical model of the hadron spectrum [193] as a superposition of thermal hadron distributions 
and collective flow. The average hadron transverse momentum, (p^) = 0.55 GeV/c, was fixed in the 
model. The analysis was done for two estimates of the charged particle multiplicity at y = in central 
collisions, dN^/dy{y = 0) = 3000 and 8000. The calorimeter occupancy is high enough in these 
cases: mean energy per tower is about 4.4 (1.7) GeV for dN^/dy{y = 0) = 8000 (3000) in the barrel 
part, and it increases by a factor ~ 2 in endcaps. Fig. 78 shows the correlation between reconstructed 
and generated transverse energies of jets in Pb— Pb and pp events. The generated jet cone has a radius 
R = 0.5 while the radius of the reconstructed cone is larger, R = 0.6. Since the average measured jet 
energy in Pb— Pb collisions is the same as in pp, the pp interactions are a baseline for jet physics in heavy 
ion collisions. However, the jet transverse energy resolution is degraded by a factor ~ 2 in the high 
multiplicity central Pb— Pb collisions compared to pp interactions, as shown in Fig. 79. The jet energy 
resolution is defined here as a{E^'^° /E^^)/ < E^'^°/E^'^ >, where S^^™ is the reconstructed trans- 
verse energy, and E^'^ is the transverse energy of all generated particles inside the given cone radius R. 
The resolution of 75 GeV jets with dN^ /dy{y = 0) = 8000 is smaller than for 100 GeV jets (as well 
as the resolution of 50 GeV jets with dN^ /dy{y = 0) = 3000 is smaller than for 75 GeV jets) since the 
background is not fully subtracted and, as a result, the reconstructed jet energy is larger in Pb— Pb than 
pp. We can define the purity of jet reconstruction as the number of events with a true QCD jet divided 
by the number of events with reconstructed jets. Then, for example, the purity is ~ 50% for 50 GeV jets 
for dN^/dy{y = 0) = 8000, because the average energy of the false jets in the background events is 
also ~ 50 GeV. The purity increases rapidly with increasing and becomes ~ 100% at 100 GeV. 

Table 7 summarizes the purity, contamination levels (false jets /generated jets) and jet transverse energy 
resolution in central Pb— Pb collisions with dN^ /dy{y = 0) = 3000 and 8000. 

Since the azimuthal angle and rapidity distribution of jets are of particular interest for jet quench- 
ing observables, the angular resolution is important. Figures 80 and 81 show the differences in azimuthal 
angle tp between generated and reconstructed 100 GeV jets in events without and with Pb— Pb back- 
ground. Even in the most pessimistic case, dN^ /dy{y = 0) = 8000, the if resolution is degraded only 
by a factor of 1.6 in Pb— Pb compared to pp collisions. The resolution is still less than the azimuthal 
size of a calorimeter tower, A^p = 0.087. A similar result is also found for the -q resolution (shown in 
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Fig. 80: Distribution of differences in azimuthal angle 
ip between generated and reconstructed jets with Et = 
100 GeV in pp events. 
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Fig. 82: Distribution of differences in pseudorapidity 77 
between generated and reconstructed jets with Et = 
100 GeV in pp events. 




Fig. 81: The same as in Fig. 80 but for Pb— Pb events at 

dN^/dy{y = 0) = 8000. 
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Fig. 83: The same as in Fig. 82 but for Pb— Pb events at 

dN^/dy{y = 0) = 8000. 
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Fig. 84: Energy dependence of jet azimuthal angle reso- Fig. 85: Energy dependence of jet pseudorapidity reso- 
lution in Pb— Pb and pp events. lution in Pb— Pb and pp events. 

Figs. 82 and 83), which is however somewhat worse than the tp resolution. The reason for that is the 
vertex was not fixed here and there was not correction of rj due to fluctuation of Z coordinate in "pile 
up" subtraction algorithm (the latest version of the algorithm included this facility). 

Figures 84 and 85 shows the energy dependence of the spatial resolution for pp and Pb— Pb events 
with dN^/dy{y = 0) = 3000 and 8000. Thus the spatial position of a hard jet can be reconstructed in 
heavy ion collisions at CMS with high enough accuracy for analysis of jet production as a function of 
azimuthal angle and pseudorapidity. 

The Level 1 single jet and electron/photon trigger rates in Pb— Pb collisions have been estimated 
using the trigger algorithms developed for pp collisions with a parameterization of HIJING results for the 
background [82]. The dominant contribution to the trigger rate comes from the single jet trigger which 
uses the transverse energy sums (electromagnetic + hadronic) computed in the calorimeter regions (4x4 
trigger cells) Ar] x Acj) = 0.348 x 0.348. For a threshold of 40—50 GeV it gives an acceptable output 
rate of about 400—200 Hz and is fully efficient for most central collisions. Assuming that with the high 
level trigger full jet reconstruction is possible, the rate can be further reduced to a level lower than 10 Hz 
for jets with reconstructed transverse energies greater than 100 GeV. The rate of the single photon trigger 
is less than 1 Hz for a threshold of 50 GeV. With such a threshold, the trigger efficiency is close to 100% 
for 7+jet events useful for off-line analysis. 

5.24 Tracking 

Track finding in heavy ion collisions is difficult due to the large number of tracks in an event. We 
consider heavy ion multiplicities up to the worst case, 8000 charged particles per unit rapidity in a central 
Pb— Pb event. In addition to the primary tracks, the CMS tracker is occupied by secondaries produced by 
interactions with the detector material. The CMS track reconstruction algorithm, originally developed 
for pp collisions, is based on Kalman Filtering and includes seed generation, track propagation, trajectory 
updating and smoothing. However, this track finder, optimized for the highest efficiency at low density, 
fails when dealing with central Pb— Pb events due to the large hit combinatorics. Therefore, modification 
of the algorithm is necessary [239]. 

The essential difference for heavy ion collisions is that the primary vertex can be determined with 



dispersion of 200 //m using only two barrel layers before any tracking. Using this constraint reduces 
the combinatorial background during track seeding. Thus these two features, primary vertex finding and 
restriction of the vertex region, were added to the standard package. 

In order to achieve good rejection power to fake tracks from random hit combinations we need to 
require the track to have as many hits as possible. Only tracks that leave the tracker through the outermost 
layer are considered. This requirement leads to a minimum transverse momentum cutoff ofpr > I GeV 
for the track to be considered reconstructable. 

Given this constraint the modified pp reconstruction algorithm gives about 80% geometrical 

acceptance and close to « 100% algorithmic reconstruction efficiency in a low multiplicity environment. 
The acceptance varies slightly ±10% over the ?7-coverage of the tracker and is independent of px 

We have tested track propagation algorithms for the case of high occupancy using Monte Carlo 
tracks seeds. The estimated efficiency under these conditions appears high enough, 80%. The momen- 
tum resolution is less than 1% at px < 100 GeV. Thus track propagation at high density is quite effective. 
The reconstruction time is around 1500 s/per event. We have also investigated using a regional track re- 
construction in a limited rj — ip area. This option is essentially needed for jet finding and correcting. The 
efficiency is very high, > 95%, but we see a large number of fake tracks. 

Thus we believe that pattern recognition is possible in heavy ion colUsions with the CMS tracking 
system. The existing pp track reconstruction package has been shown to be a valuable framework, but 
needs major restructuring to reduce combinatorics and computation time. 



5.25 Impact Parameter Determination 

It is important to study hard jet and high-py hadron production in heavy ion colfisions as a function of 
centrality. In CMS, the best way to determine the impact parameter, b, event-by-event is the transverse 

energy deposited in the calorimeters, E^^, which strongly decreases from central to peripheral colli- 
sions [82]. The jet production rates, N^'^^, can be measured in bins of E'^?*. Then the b— and E^^— 
dependencies of ATJ^* can be related by the E^^ — b correlation functions Caa, 



E^^-E^^\b) 



Ni'\b) = J dEl^Ni^{E!f )CAA(h,E!f ), C^a(6,£¥') = ^=-i^exp 

Since very forward rapidity region, |?7| > 3, is almost free of final state interactions, the (trans- 
verse) energy deposition in HF is determined mainly by the initial nuclear geometry of a collision rather 
than by final state dynamics. Determining the impact parameter via (transverse) energy deposition in 
the HF calorimeter thus avoids some possible uncertainties [82, 113]. The impact parameter dependence 
of the total transverse energy produced in the pseudorapidity interval 3 < I77I < 5 obtained with HI- 
JING [270, 155] is presented in Fig. 86 from [113] for Pb— Pb interactions. The Et — b correlation is 
diffuse due to fluctuations in the nucleus-nucleus collision dynamics, including fluctuations in the num- 
ber of nucleon-nucleon sub-collisions at a given b and fluctuations of transverse energy flow in each 
nucleon-nucleon interaction. The correlation of the total energy flow is of similar shape. Most of the 
energy produced in the very forward direction is between 10 and 100 TeV. It is then possible to mea- 
sure the total energy with high accuracy, reducing uncertainties in b. The impact parameter distribution 
at fixed values of Et is Gaussian-like with a width, a^,, dependent on impact parameter, see Fig. 87 
from Ref. [113]. The absolute accuracy is defined here as =b2crb, about ~ 1 fm for Pb— Pb colfisions 
with 1 < ^ < 13 fm. It is degraded by a factor of ~ 2 for very peripheral events, 5 > 13 fm, due to 
the diminution of the produced energy in the pseudorapidity region. At the same time, we see that the 
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Fig. 86: Correlation between the transverse energy flow, Et, and impact parameter b at very forward rapidities, 
3 < < 5, for 10000 minimum bias Pb-Pb coUisions [113]. 
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Fig. 87: Impact parameter dependence of (a) the Gaussian width ai, and (b) the relative error 6 = 2ai,/b for Pb— Pb 
collisions [113]. 



relative error is minimal for peripheral collisions since the statistics are increased. 

To summarize, we note that although these results were obtained on the particle level, similar 
conclusions are expected to be valid when detector effects are taken into account. It has been shown that 
the finite energy and spatial resolution of the HF calorimeter do not substantially degrade the accuracy 
of the impact parameter determination in heavy ion collisions [113]. 



5.26 Reconstruction of Nuclear Reaction Plane 

The azimuthal anisotropy of jet and high-py particle production in semi-central heavy ion collisions is 
predicted to be a signal of partonic energy loss in an azimuthally asymmetric volume of quark-gluon 
plasma [205, 277, 159, 158]. The advantage of azimuthal jet observables is that one needs to reconstruct 
only the azimuthal position of jet, not its total energy. It can be done easily and with high accuracy while 
reconstruction of the jet energy is more ambiguous. The methods summarized in Ref. [267, 234] present 
ways to determine the reaction plane. They are applicable to the study of anisotropic flow of soft and 
semi-hard particles in the current dedicated heavy ion experiments at the SPS [54] and RHIC [13] and 
may also be used at the LHC [205]. When the azimuthal distribution of particles is described by the 
elliptic form, 

dN No . ^ . M ? J dN 

= — [1 + 2i;2 cos 2((/J - (^Creac)] , ^ j '^^ ' 

— TT 

the nuclear reaction plane angle, </?reac5 is 

LVi sin 2ipi 



tan (2(/3reac) 



where the weights, ui, are selected to optimize the resolution. The coefficient V2 of the azimuthal 
anisotropy of particle flow is an average over cos(2<^). In CMS, the weights can be introduced [209] 
as energy deposition in calorimeter sector i of position ipi, toi = Ei{(pi). Fig. 88 illustrates the energy 
deposition in sectors of the barrel and endcap regions of the CMS hadronic and electromagnetic calorime- 
ters, \rj\ < 3, for generated with hydro-code [193, 207] Pb— Pb events at 6 = 6 fm. The detector response 
is obtained with CMSIM_125. The estimated resolution of the (^reac determination, a^p reac = 0.15 rad, 
see Fig. 89, allows measurement of the coefficient of the jet azimuthal anisotropy with ~ 90% accuracy, 
defined as the ratio of the average of (cos 2(/?jet) over all events "measured" to its "true" value. The 
estimates obtained are quite optimistic because hydrodynamic models give rather large values of elliptic 
flow at high-pT- On the other hand, the majority of microscopic Monte Carlo models underestimate 
flow. For example, under the same conditions HUING [270, 155] with jet quenching predicts much less 
anisotropic flow and yields aip reac = 0.8 with only ~ 20% accuracy. 

Recently a method was also suggested [207] for measuring the jet azimuthal anisotropy coeffi- 
cients without event-by-event reconstruction of the reaction plane. This technique is based on the corre- 
lations between the azimuthal position of the jet axis and the angles of particles not incorporated in the 
jet. Then 

V'^ = (cos 2(^jet)event = ( , ^ / ^ ) 

\ V(C0S2((^1 - ^2) Wl(^l) '^2m)) I ^^^^^ 

where the weights Ui are defined as before. In some sense, this represents the development and gener- 
alization of the well-known method for measuring the azimuthal anisotropy of particle flow originally 
considered in a number of works, see for example Refs. [267, 234, 269, 228, 229]. The accuracy of the 
method improves with increasing multiplicity and particle (energy) azimuthal anisotropy and is practi- 
cally independent of the absolute value of the azimuthal anisotropy of the jet itself [207]. The accuracy 
of achieved by such a method are estimated to be 94% for hydrodynamics and 30% for HUING, 
somewhat better than those obtained from direct reconstruction of the reaction plane angle. 




^(gen)-^(rec), rod 

Fig. 88: Energy deposition in the barrel and endcap re- Fig. 89: Distribution of differences in the generated and 
gions of the CMS hadronic (HCAL) and electromag- reconstructed azimuthal angles of the Pb— Pb reaction 
netic (ECAL) calorimeters for Pb— Pb collisions at 6 = plane for 6 = 6 fm (hydrodynamics with CMSIM_125). 
6 fm (hydrodynamics with CMSIM_125). 

5,3 Jet Physics with the ATLAS Detector 

S. Aronson, K. Assamagan, B.Cole, M. Dobbs, J. Dolesji, H. Gordon, F. Gianotti, S.Kabana, M.Levine, F. 
Marroquim, J.Nagle, P. Nevski, A. Olszewski, LRosselet, P. Sawicki, H.Takai, S. Tapprogge, A. Trzupek, 
M.A. B. Vale, S.White, R. Witt, B. Wosiek and K. Wozniak 

The ATLAS detector is designed to study high pT physics in proton-proton collisions at high 
LHC machine luminosity. Most of the detector subsystems will be available for the study of heavy ion 
collisions. One of the highlights of the ATLAS detector is its calorimeter subsystem. Both the elec- 
tromagnetic and hadronic compartments are finely segmented and well suited for jet quenching studies. 
RHIC results suggest that partons may radiate gluons in the dense matter formed in heavy ion collisions. 
This phenomena can be certainly be well explored in the ATLAS detector. We report on early assessment 
of the detector capabilities in the heavy ion environment. 

5.31 The ATLAS Detector 

The ATLAS detector is designed to study proton-proton collisions at the LHC design energy of 14 TeV 

in the center of mass. The physics pursued by the collaboration is vast and includes: Higgs boson search, 
searches for SUSY, and other scenarios beyond the Standard Model, as well as precision measurements 
of process within (and possibly beyond) the Standard Model. To achieve these goals at a full machine 
luminosity of 10^^ cm~^s~^. ATLAS will have a precise tracking system (Inner Detector) for charged 

particle measurements, an as hermetic as possible calorimeter system, which has an extremely fine grain 
segmentation, and a stand-alone muon system. An overview of the detector is shown in Fig. 90. 

The Inner Detector is composed of (1) a finely segmented silicon pixel detector, (2) silicon strip 
detectors (Semiconductor Tracker (SCT)) and (3) the Transition Radiation Tracker (TRT). The segmen- 
tation is optimized for proton-proton collisions at design machine luminosity. The inner detector system 
is designed to cover a pseudo-rapidity of | r/ |< 2.5 and is located inside a 2 T solenoid magnet. 

The calorimeter system in the ATLAS detector surrounds the solenoid magnet is divided into 
electromagnetic and hadronic sections and covers pseudo-rapidity | r] |< 4.9. The EM calorimeter is 




Fig. 90: The overall layout of the ATLAS detector 



an accordion liquid argon device and is finely segmented longitudinally and transversely for | r/ |< 3.1. 
The first longitudinal segmentation has a granularity of 0.003 x 0.1 (A?7 x Ac/)) in the barrel and slightly 
coarser in the endcaps. The second longitudinal segmentation is composed of Aij x Ac/) = 0.025 x 0.025 
cells and the last segment A?; x = 0.05 x 0.05 cells. In addition a finely segmented (0.025 x 0.1) pre- 
sampler system is present in front of the electromagnetic (EM) calorimeter. The overall energy resolution 
of the EM calorimeter determined experimentally is 10% /^/E © 0.5%. The calorimeter also has good 
pointing resolution, dOmrad/V^ for photons and timing resolution better than 200 ps for showers of 
energy larger than 20 GeV. 

The hadronic calorimeter is also segmented longitudinally and transversely. Except for the end- 
caps and the forward calorimeters, the technology utilized for the calorimeter is a lead-scintillator tile 
structure with a granularity of At] x Ac/> = 0.1 x 0.1. In the endcaps the hadronic calorimeter is imple- 
mented in liquid argon technology for radiation hardness with the same granularity as the barrel hadronic 
calorimeter. The energy resolution for the hadronic calorimeters is 50% /^/E 2% for pions. The very 
forward region, up to ?? = 4.9 is covered by the Forward Calorimeter implemented as an axial drift liquid 
argon calorimeter. The overall performance of the calorimeter system is described in [1]. 

The muon spectrometer in ATLAS is located behind the calorimeters, thus shielded from hadronic 
showers. The spectrometer is implemented using several technologies for tracking devices and a toroidal 
magnet system, which provides a field of 4 T strength to have an independent momentum measurement 
outside the calorimeter volume. Most of the volume is covered by MDTs, (Monitored Drift tubes). The 
forward region where the rate is high. Cathode Strip Chamber technology is chosen. The stand-alone 
muon spectrometer momentum resolution is of the order of 2% for muons with px in the range 10 - 100 
GeV. The muon spectrometer coverage is | 77 |< 2.7. 



The trigger and data acquisition system of ATLAS is a multi-level system, which has to reduce the 
beam crossing rate of 40 MHz to an output rate to mass storage of 0(100) Hz. The first stage (LVLl) 
is a hardware based trigger, which makes use of coarse granularity calorimeter data and dedicated muon 
trigger chambers only, to reduce the output rate to about 75 kHz, within a maximum latency of 2.5 /is. 

The performance results mentioned have been obtained using a detailed full simulation of the 
ATLAS detector response with GEANT and have been validated by an extensive program of testbeam 
measurements of all components. 

5.32 Jet Physics and ATLAS 

The ATLAS calorimer coverage and its fine segmentation will be an asset for jet studies in the heavy 
ion environment. Signatures of jet quenching in central heavy ion collisions could manifest as a larger 
jet cone (as compared to proton-proton collision) and/or modifications in the jet fragmentation func- 
tion. The finely segmented (longitudinally and transversely) electromagnetic calorimeter will allow us 
to reconstruct EM clusters in the jet environment and in particular measure the tt" containt in the jet. 

There is excellent opportunity in ATLAS to measure 7 —jet, jet-jet and Z-jet events where one can 
more fully characterize the modified fragmentation functions. In particular, the 7(or Z) in j{orZ) — jet 
processes provides a "control" over the away-side jet energy and direction that will allow the physics of 
quenching to be studied quantitatively and in great detail [274]. The effects of hard gluon radiation on 
the photon/jet energy imbalance and angular distribution can be studied in great detail using the high- 
statistics p-p data set. The 7 — jet channel requires the identification of a photon. In proton-proton 
collisions the rejection of 7/7r° is about a factor of three up to a pt of 50 GeV. However, the heavy ion 
environment presents considerable more challenge. Zq production rates have been estimated by Wang 
and Huang [274]. For px larger than 40 GeV, we expect of the order of ~500 Zq /i^/i^ events for one 
month-run. Therefore multiple runs may be required to extract relevant information on jet fragmentation. 

5.33 Expected Detector Performance for Jet Studies 

The ATLAS calorimeter granularity is shown in Table 1. The calorimeter system is not only segmented 
in r] and (j), but also longitudinally. The calorimeter is fully segmented to 77 = 5. The electromagnetic 
calorimeter is in most places 25Xo deep and designed to fully contain a 1 TeV electron or photon. 
The hadronic calorimeter is more than IQXint deep and contains all of the hadrons in typical pp high 
luminosity runs. 



Table 1: ATLAS calorimeter system segmentation. Listed are the number of longitudinal and the size of 
the transverse segmentation in the different calorimeters. 



Calorimeter System 


r] coverage 


Long. 


Transversal segmentation 


LAr Electromagnetic 


0.0 < ry < 3.2 


3 


0.003 X 0.1, 0.025 X 0.025, and 0.05 x 0.05 


LAr Hadronic 


1.5 < r? < 3.2 


4 


0.1.1 for 1.5 < 77 < 2.5 and 0.2 x 0.2 otherwise. 


Hadronic Tile 


0.0 < ?7 < 1.7 


3 


0.1 x 0.1, 0.1 X 0.1, and 0.2 x 0.1 


Forward Calorimeter 


3.K ?7 < 4.9 


4 


0.2 X 0.2 



The jet energy resolution for high luminosity proton-proton run is hQ%/\fE 2% and the jet 
reconstruction threshold approximately Et ~ 2{)GeV for proton-proton run. These numbers are ex- 
pect to be different for the heavy ion environment. Preliminary simulations indicate that the detection 
thresholds should be at around 40~50 GeV. We expect a worsening of the energy resolution because of 
the soft background but the energy scale should remain untouched, unless the ratio of EM to Hadronic 
components in the jet changes substantially. 

Preliminary figures for energy deposition in the EM calorimeter from full HIJING central events 
indicate that approximately 4 GeV of transverse energy is deposited in a tower r/X(/> = 0.1x0.1. 



This number is consistent with reported by CMS. However, due to the longitudinal segmentation and 
predominantly low pr (< I Gey)nature of the particles in the background, we do expect that most of 
the energy to be deposited in the first compartament of the calorimeter. Thus jets could be reconstructed 

on the basis of the remaining compartaments. Detailed studies are under way. 

To study jet quenching in a direct way is to measure its fragmentation function and possible 
changes in the jet cone radius. The fragmentation function can only be measured if particles are identified 
within the jet. The segmentation of the ATLAS calorimeter is such that allows, in principle, for the identi- 
fication of TTo's. In spite of the soft background preliminary studies shows encouraging results. However, 
we have observed a significant worsening of the EM cluster energy resolution. Studies performed in 
proton-proton collisions do have not addressed the issue of EM cluster reconstruction, specially for vr'^s 
at low energy. 
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